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Preface 


Strategic  airlift  air-refueling  planning  is 
becoming  more  important  as  the  MAC  airlift  fleet  gains  an 
operational  in-flight  refueling  capability.  As  fuel  sup¬ 
plies  become  more  expensive  and  scarce,  fuel  conservation 
becomes  increasingly  more  important.  We  sincerely  hope 
that  airlift  and  tanker  operational  planners  will  use  the 
results  of  this  research  when  selecting  air-refueling 
points  and  takeoff  fuel  loads  to  reduce  total  fuel  con¬ 
sumption. 
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Abstract 


During  the  1980s  increases  in  the  potential  use  of 
Strategic  Airlift  to  transport  equipment  and  personnel  is 
anticipated.  The  capabilities  of  Strategic  Airlift  air¬ 
craft  are  extended  through  the  inclusion  of  efficient 
in-flight  refueling.  The  primary  objective  of  this  research 
was  to  develop  a  method  which  determines  the  combination  of 
in-flight  refueling  rendezvous  point,  takeoff  fuel  loads 
and  tanker  base  which  results  in  the  minimum  total  fuel  con¬ 
sumption  for  an  airlifter  and  tanker  aircraft. 

The  experimental  design  included  the  creation  of 
two  models.  An  analytic  flight  planning  model  determined 
the  optimal  rendezvous  point  and  the  takeoff  fuel  loads 
for  the  aircraft  in  a  specific  mission  scenario.  A  SLAM 
•simulation  model  verified  the  operational  feasibility  of 
the  results  of  the  analytic  model  by  simulating  the  flights 
of  aircraft. 

The  optimal  rendezvous  can  only  be  determined  by 
analyzing  the  interaction  of  the  airlifter  route  distance, 
the  cargo  load,  and  the  location  of  the  tanker  base.  As 
the  total  distance  or  cargo  load  decreases  or  the  tanker 
base  is  located  farther  along  the  airlifter 's  route  of 
flight,  the  optimal  rendezvous  point  is  located  incre¬ 
mentally  farther  from  the  boundary  established  by  the 
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maximum  feasible  range  of  the  airlifter  from  its  destina¬ 
tion.  The  optimal  takeoff  fuel  loads  are  dependent  on  the 
aircraft  combination  and  will  result  in  the  smallest  sum 
of  the  total  fuel-carrying  capacity  used.  By  using  the 
optimal  rendezvous  point  for  an  in-flight  refueling  mission, 
significant  fuel  savings  are  realized. 
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OPTIMIZATION  OF  STRATEGIC  AIRLIFT  IN-FLIGHT  REFUELING 


I .  Introduction 

The  Problem 

Efficient  in-flight  refueling  of  strategic  airlift 
aircraft  is  being  recognized  as  vital  to  the  Military  Air¬ 
lift  Command  (MAC)  and  the  Department  of  Defense.  Cur¬ 
rently,  all  C-5A  aircraft  are  air-refuelable  and  by  the  end 
of  1982  all  MAC  C-141  aircraft  will  be  air-refuelable 
(Ref  21:48).  As  a  result  of  this  increased  capability, 

MAC  is  constantly  striving  through  innovation  and  research, 
to  make  air-refueling  missions  more  productive  and  effi¬ 
cient.  One  way  to  increase  the  efficiency  of  airlift  air¬ 
refueling  missions  is  to  determine  the  combination  of  plan¬ 
ning  factors;  the  takeoff  fuel,  tanker  basing,  and  air¬ 
refueling  rendezvous  point,  which  would  minimize  the  total 
fuel  consumed  by  an  airlift  and  tanker  aircraft.  No  current 
method  determines  an  optimal  combination  of  these  factors 
for  in-flight  refueling  between  strategic  airlift  and  tanker 
aircraft . 

Objectives 

The  objective  of  this  research  effort  was  to 
develop  a  method  which  determines  the  optimal  combination 
of  in-flight  refueling  rendezvous  point,  takeoff  fuel  loads 


and  tanker  base  for  a  refueled  aircraft  lAission.  By 
determining  the  optimal  combination  three  questions  are 
answered.  First,  should  the  rendezvous  point  be  close  to  or 
far  from  the  points  of  departure?  Second,  how  much  fuel 
should  be  loaded  on  the  airlifter  and  tanker  prior  to  take¬ 
off?  Third,  from  what  base  should  the  tanker  depart? 

When  the  research  was  begun  it  was  believed  that 
the  following  two  hypotheses  would  prove  true. 

1.  The  minimum  total  fuel  consumed  by  the  airlift 
and  tanker  aircraft  for  their  combined  flight  will  result 
from  a  rendezvous  point  located  at  the  maximum  flight  range 
of  the  airlifter  from  its  destination  base.  This  point 

is  always  located  on  the  boundary  of  the  region  of  feas¬ 
ible  rendezvous  points  closest  to  the  airlifter  takeoff 
base . 

2.  Airlifter  aircraft  departures  with  the  maximum 
allowable  fuel  load  will  always  result  in  the  minimum  total 
fuel  consumprion  for  both  aircraft.  This  implies  that  the 
fuel  transferred  is  the  minimum  required  to  complete  the 
flight. 

By  demonstrating  the  accuracy  of  these  hypotheses 
they  could  then  be  used  to  help  define  a  general  policy 
for  efficient  in-flight  refueling. 

Two  models  were  developed  in  the  course  of  this 
research.  The  first,  FLTPLN,  computed  the  optimal  plan¬ 
ning  factors  for  an  air-refueling  airlift  mission.  The 
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FLTPLN  results  were  then  validated  for  operational  suit¬ 
ability  by  a  SLAM^  (Simulation  Language  for  Alternative 
Modeling)  simulation  model. 

These  two  models  were  used  to  explore  different 
scenarios  and  many  fuel  cargo  combinations  for  the  airlift 
aircraft.  Aircraft  combinations  investigated  are  listed 
in  Table  1.  Sensitivity  analysis  was  used  to  determine 
the  effect  on  the  response  variable,  total  fuel  consumption, 
to  changes  in  the  input  factors  of  takeoff  fuel  onload, 
cargo,  wind,  takeoff  delays,  and  total  distance. 


TABLE  1 

AIRCRAFT  COMBINATIONS 


Airlif ters 

Tankers 

C-5A 

KC-135 

C-5A 

KC-10 

C-141B 

KC-135 

C-141B 

KC-10 

Benefits 

Three  benefits  are  apparent  from  deteirmining  an 
optimal  combination  of  planning  factors  which  minimizes 
total  fuel  consumed.  The  first  benefit  is  obviously  fuel 
conservation.  To  emphasize  the  importance  of  conserving 

^SLAM  is  a  simulation  language  developed  by 
Pritsker  and  Associates,  Inc.,  West  Lafayette,  Indiana. 
SLAM  has  been  designed  to  support  engineers,  managers, 
and  researchers. 
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fuel,  the  Chairman  of  the  Joint  Chiefs  of  Staff  stated  that, 
the  supply,  access  and  price  of  oil  will  dominate  future 
national  security  issues  (Ref  13:i).  By  minimizing  total 
fuel  consumed,  operating  costs  go  down,  and  funds  used  in 
procuring  the  saved  fuel  could  be  used  for  other  equally 
essential  programs.  As  the  cost  of  fuel  goes  up,  the 
savings  from  an  optimal  refueling  policy  will  also  go  up. 

For  example,  the  price  of  fuel  (JP-4)  increased  from  $.42 
per  gallon  in  1979  -to  $1.16  per  gallon  in  1980  (Ref  19). 
Additional  increases  can  be  expected  in  the  future.  The 
second  benefit  which  results  from  an  optimal  refueling 
point  is  that  individual  mission  flight  times  are  reduced. 
While  the  flight  time  of  modern  jet  aircraft  results  in  a 
relatively  short  deployment  time,  the  savings  from  the 
multiple  sorties  of  a  large-scale  deployment  would  be  sig¬ 
nificant.  This  means  closure  times  can  be  reduced  due  to 
decreased  flight  times.  Closure  time  occurs  when  the  lat»t 
aircraft  of  a  deployment  arrives  at  the  destination.  More 
importantly,  the  airframe  time,  which  cannot  be  renewed, 
is  also  reduced.  The  third  benefit  from  an  optimal  refuel¬ 
ing  point  is  that  operational  plans  and  contingency  plans 
can  be  optimally  formulated  for  a  wide  range  of  mission 
scenarios.  By  combining  all  three  benefits,  the  advan¬ 
tages  of  determining  an  optimal  refueling  point,  which 
ensures  more  efficient  airlift  operations  are  apparent. 

During  the  1980s,  significant  increases  in  the 
potential  use  of  airlift  to  transport  equipment  and 
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supplies  worldwide  is  anticipated.  Due  to  political  and 
natural  constraints  the  importance  of  air-refueling  during 
a  rapid  deployment  will  also  increase.  Therefore,  the 
efficient  use  of  air-refueling  of  airlift  aircraft  will 
result  in  significant  economies  to  the  Air  Force. 


II.  Background 


Military  Airlift  Command  policies  regarding  air¬ 
refueling  operations  for  airlift  aircraft  are  relatively 
new  and  untested,  compared  to  those  of  SAC  and  TAC. 
Employment  plans  utilizing  air-refueling  are  founded  on 
only  a  few  years  of  experience,  and  have  not,  as  yet,  been 
tested  with  large-scale  deployments  or  war-scenario  exer¬ 
cises.  Employment  policy  and  operations  concepts  for  air¬ 
refueling  are  constantly  being  formulated  and  updated  as 
more  research  is  conducted  in  this  area.  Chapter  II  pro¬ 
vides  a  background  in  air-refueling  history,  employment 
principles  and  policy  currently  in  use  at  HQ  MAC.  Included 
is  a  description  of  the  extent  of  MAC  air-refueling 
experience,  active  air-refueling  policy  research,  current 
planning  models  in  use  at  MAC  and  HQ  USAF,  and  a  descrip¬ 
tion  of  mission  flight  planning,  using  a  typical  airlift 
flight  profile. 

History  of  Airlift  Air-Refueling 

In  the  last  two  decades,  strategic  airlift  has 
become  increasingly  important  to  the  Department  of  Defense 
to  support  the  rapid  deployment  and  resupply  of  Allied  or 
United  States  forces  overseas.  To  increase  this  resupply 
capability,  C-5A  aircraft  were  acquired  and  C-141  aircraft 
are  being  modified  to  utilize  in-flight  refueling. 


In-flight  refueling  became  availrh to  the  Mili¬ 
tary  Airlift  Command  in  December,  1969  (Ref  12) ,  with  the 
introduction  of  the  C-5A  aircraft  into  the  Air  Force  inven¬ 
tory.  Although  the  C-5  was  originally  built  to  be  air- 
refuelable,  this  capability  was  not  tested  and  evaluated 
until  May,  1974,  with  the  final  evaluation  report  submitted 
in  March,  1975  (Ref  5:ii,2).  Presently,  MAC  is  modifying 
C-141A  aircraft  for  air-refueling  operations  and  will 
change  the  designation  to  C-141B. 

The  need  for  quick  and  efficient  airlift  and  the 
potential  need  for  air-refueling  airlift  aircraft  was 
graphically  portrayed  during  the  Arab-Israeli  Conflict  of 
1973.  In  October,  1973,  Egypt  and  her  allies  confronted 
Israel  with  a  pre-dawn  surprise  attack  (Ref  9:27-28). 

The  initial  impact  on  Israel  was  a  Ic.rge-scale  loss  of 
equipment  (aircraft  and  tanks)  as  well  as  a  large  expendi¬ 
ture  of  munitions.  During  the  first  few  days  of  the  war, 
the  United  States  realized  that  Israel  would  require  mas¬ 
sive  and  rapid  resupply  to  avert  a  major  catastrophe.  Due 
to  the  rapid  delivery  requirement,  the  only  satisfactory 
delivery  method  was  airlift.  Lack  of  support  from  our 
European  Allies  compounded  the  delivery  problem,  as  they 
withheld  overflight  and  landing  rights  to  Israel-bound  air- 
lifters.  As  a  result,  equipment  and  munitions  stockpiled 
in  Europe  were  not  available  for  Israeli  use  and  all  sup¬ 
plies  and  munitions  destined  for  Israel  had  to  come  from 
the  Continental  United  States  (CONUS)  (Ref  15:14-15). 
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The  only  European  ally  to  offer  support  to  the  United 
States  was  Portugal,  which  controlled  the  Azores  in  the  mid 
Atlantic  Ocean.  Portugal  approved  the  Azores  as  a  MAC 
staging  base  enabling  Israel  to  be  resupplied  by  air  with¬ 
out  air-refueling  (Ref  15:14-15).  In  33  days,  MAC  crews 
airlifted  22,315  tons  of  equipment  and  supplies,  flying 
421  C-141  sorties  and  145  C-5A  sorties  (Ref  8:26-32). 

The  quick  and  effective  resupply  of  Israel's  military 
machine  enabled  Israel  to  regain  the  offensive  on  all 
borders  and  eventually  dominate  the  war  (Ref  15:14-15). 

The  Air  Force  learned  a  very  important  lesson  from 
this  resupply  effort.  Our  strategic  airlift  fleet  is  very 
dependent  on  friendly  foreign  nations  to  provide  staging 
bases.  Without  Portugal's  cooperation,  the  airlift  would 
not  have  been  possible  since  the  C-5A  did  not  have  an  opera 
tional  air-refueling  capability  at  that  time.  To  meet 
future  challenges  to  the  United  States  or  her  Allies,  the 
dependence  on  staging  bases  had  to  be  reduced  through  the 
expansion  of  air-refueling  capabilities. 

Since  1973,  the  Air  Force  has  undertaken  several 
new  programs  to  increase  the  in-flight  refueling  capabil¬ 
ity  of  the  strategic  airlift  fleet.  First,  the  C-5A  was 
made  operationally  air-refuelable  and  aircrews  were 
trained  (Ref  5) .  Second,  234  C-141A  aircraft  are  being 
modified  to  air-refuelable  C-141B  aircraft  (Ref  13:199). 
Third,  a  new  tanker  fleet  of  DC-lOs  (designated  KC-10) 
has  been  purchased.  With  the  addition  of  the  KC-10, 
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multiple  tankers  will  no  longer  be  require-'  to  refuel  the 

2 

C-5  on  a  routine  basis.  Fourth,  a  proposal  to  install 
more  efficient  engines  in  the  KC-135  tanker  is  being 
studied.  The  addition  of  these  new  engines  and  their 
reduced  fuel  consumption  will  allow  more  fuel  to  be  trans¬ 
ferred  to  the  airlifter  or  increase  the  rendezvous  range 
of  the  KC-135  for  a  fixed  offload  of  fuel.  This  results 
in  increased  range  for  the  airlifters. 

With  all  of  these  changes,  air-refueling  of  air¬ 
lift  aircraft  will  make  a  greater  contribution  to  the  stra¬ 
tegic  and  tactical  posture  of  the  United  States.  However, 
increasing  airlift  capability  by  adding  new  refueling 
assets  (KC-10  and  C-141B)  also  requires  expertise  and 
experience  in  creating  plans  and  policy  to  employ  these 
forces . 

In  the  future,  contingency  plans  which  call  for  air¬ 
lift  resources  will  be  re-evaluated  using  air-refueling  as 
an  option.  General  Huyser,  CINCMAC,  has  stated  that  future 
airlift  plans  created  to  support  war  plans  will  address 
in-flight  refueling  as  an  option  in  order  to  take  advan¬ 
tage  of  this  increased  capability  and  efficiency  (Ref  16) . 

In  addition,  studies  to  determine  how  to  optimally  use 
strategic  airlift  refueling  must  also  be  undertaken. 

These  studies  should  research  operational  problems  as  well 

2 

The  C-5A  fuel  capacity  is  315,100  pounds  and  the 
KC-10  fuel  capacity  is  349,153  pounds  of  which  any  portion 
can  be  transferred  in-flight  to  the  airlifters.  The  KC- 
135A  fuel  capacity  is  only  165,000  pounds. 
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as  economic  and  logistic  problems.  Includi_d  in  the  opera¬ 
tional  problem  area  is  the  problem  of  obtaining  an  optimal 
air-refueling  rendezvous  point. 

Current  Concepts  and  Models 

MAC  planners  are  investigating  the  possible  incor¬ 
poration  of  air-refueling  models  in  various  war  and  opera¬ 
tions  plans. ^  This  responsibility  falls  to  the  Opera¬ 
tional  Plans  Office,  HQ  MAC/XP  (Ref  16) .  As  plans  are 
periodically  revised,  air-refueling  is  being  considered  as 
a  possible  method  of  increasing  the  efficiency  and  effec¬ 
tiveness  of  the  plans.  MAC  considers  three  areas  as  poten¬ 
tially  beneficial  in  increasing  the  efficiency  of  airlift 
operations: 

1.  When  air-refueling  is  the  only  possible  method 
available  for  a  particular  mission.  An  example  would  be  a 
flight  from  Moody  AFB,  Georgia,  to  West  Cairo  AB,  Egypt, 
in  a  C-141B.  A  C-141B  cannot  fly  from  Moody  AFB  to  Cairo 
AB  without  either  stopping  enroute  for  fuel  or  air¬ 
refueling.  If,  for  political  reasons,  a  staging  base  is 
not  available,  then  the  C-141B  would  be  unable  to  accom¬ 
plish  the  deployment.  Therefore,  air-refueling  is  the 
only  way  a  C-141B  can  be  used  on  this  mission. 

2.  When  a  technical  constraint  is  imposed  upon  an 
airlift  operation  and  air-refueling  is  required.  For 

^Specific  contingency  plans  which  address  air¬ 
refueling  are  classified  and  will  not  be  covered  here. 
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example,  if  a  C-5A  were  to  airlift  a  full  load  of  the 
Army's  battle  tanks  from  Peterson  Field,  Colorado  (Colorado 
Springs) ,  the  C-5A  would  be  unable  to  take  off  with  a  full 
fuel  load  due  to  inadequate  runway  length  (11,000  FT)  at 
the  airfield's  pressure  altitude  (6,172  FT)  (Ref  10:B96). 
Air-refueling  v/ould  add  the  necessary  fuel,  after  takeoff, 
to  continue  the  flight  to  destination. 

3.  Air-refueling  can  be  used  to  improve  closure 
time  of  a  mission  or  a  planned  deployment.  An  example 
would  be  the  deployment  of  a  tactical  fighter  squadron 
from  the  CONUS  to  Korea.  By  employing  air-refueling,  the 
improved  closure  time  may  be  decisive  in  political  negotia¬ 
tions. 

Air-refueling  employment  in  a  war  plan  or  a  large- 
scale  operations  plan  has  never  been  tested  by  MAC  and 
therefore  actual  experience  is  limited.  To  determine  the 
feasibility  and  efficiency  of  using  strategic  airlift  air¬ 
refueling  in  a  MAC  plan,  the  planning  staff  must  actually 
prepare  the  plan  twice;  once  with  the  air-refueling 
included  and  once  without,  so  that  an  intelligent  choice 
between  the  two  can  be  made.  Actual  mission  profiles  must 
be  planned  so  that  logistic  and  aircrew  requirements  can 
be  determined  and  analyzed. 

In  planning  the  air-refueling  sorties  for  a  par¬ 
ticular  mission,  coordination  between  the  MAC  staff  and 
the  SAC  staff  is  essential.  This  coordination  results  in 


a  tanker  commitment  from  SAC  for  the  particular  plan. 


Agreement  on  a  rendezvous  point  betv/een  the  two  W^JCOMs 
is  accomplished  by  a  series  of  negotiations  between  the  two 
planning  staffs  (Ref  3) .  During  the  negotiating,  each 
staff  tries  to  obtain  a  rendezvous  point  most  advantageous 
for  their  MAJCOM.  However,  during  this  process,  analyti¬ 
cal  procedures  or  models  which  can  determine  an  optimal 
route  of  flight  for  each  aircraft  are  not  used  (Ref  23) . 
Therefore,  during  the  planning  process,  neither  staff 
knows  if  the  agreed-upon  rendezvous  point  is  picked  so  as 
to  minimize  total  fuel  consumed  or  total  mission  flying 
time  for  both  aircraft. 

As  HQ  MAC/XP  is  planning  routes  of  flight  for  war 
plans,  other  routine  refueling  missions  are  also  being 
planned.  Daily  air-refueling  training  missions  that  MAC 
and  SAC  fly  ere  planned  in  the  same  iterative  manner  by 
HQ  MAC/DOOMF  staff  officers.  Coordination  and  negotiations 
with  SAC  for  an  air-refueling  rendezvous  point  is  done  in 
the  same  manner,  again  without  the  benefit  of  any  opti¬ 
mizing  procedures  (Ref  3) .  Although  these  methods  are  ade¬ 
quate  and  are  operationally  feasible,  the  planning  could  be 
done  quickly  and  more  economically  using  analytical  models. 

In  an  effort  to  obtain  knowledge  in  the  area  of 
air-refueling,  HQ  MAC  has  created  air-refueling  simulation 
models.  Air-refueling  deployments  and  operations  are 
being  analyzed  in  order  to  gain  insight  into  the  world  of 
strategic  airlift  air-refueling.  At  HQ  MAC,  in  the  recent 
past,  one  primary  model  has  been  used  to  study  airlift 


air-refueling.  The  model  is  presently  ir  se  and  is  known 
as  the  M-14  model,  commonly  referred  to  as  "Collosus" 

(Ref  14) .  The  name  M-14  identifies  it  as  the  fourteenth 
macro-model  created  at  HQ  MAC.  This  model  was  originally 
created  without  air-refueling,  but  this  feature  has  been 
recently  added.  This  Monte  Carlo  simulation  model  simu¬ 
lates  the  entire  Military  Airlift  Command  operations  and 
incorporates  all  strategic  airlift  aircraft  presently 
owned  by  MAC.  M-14  includes  factors  such  as  crews,  air 
bases,  material  depots,  and  logistic  material  used  in 
daily  MAC  operations  (Ref  14) .  The  purpose  of  M-14  is  to 
study  the  entire  airlift  capacity  of  MAC  under  almost  any 
circumstances.  Other  uses  for  the  M-14  include  study  and 
analysis  of  time-phased  deployments,  such  as  JCS  exercises 
(REFORGER,  JACKFROST)  and  to  study  the  interaction  effects 
of  factors  on  airlift  operations  (crews,  supplies,  etc.). 
M-14  can  also  be  used  to  determine  or  test  optimal  global 
airlift  strategies  in  the  event  of  a  global  conflict 
requiring  massive  use  of  airlift.  Evaluation  and  valida¬ 
tion  of  war  and  operations  plans  can  also  be  accomplished 
with  this  model.  M-14  is  an  extremely  large  model  and 
cannot  be  run  at  HQ  MAC  due  to  a  lack  of  computer  capa¬ 
bility  (Ref  14) .  As  a  consequence,  a  larger  computer  at 
Kirtland  AFB,  NM,  must  be  used  whenever  the  M-14  is  used. 

Major  initial  inputs  to  the  model  include  departure 
base,  destination,  route  of  flight,  and  fuel  and  cargo 
loads.  Another  important  variable  is  tanker-air lifter 


air-refueling  rendezvous  point  for  air-rc fuel ing  missions. 
The  M-14  model  allows  the  manual  input  of  a  rendezvous 
point,  or  it  will  calculate  a  rendezvous  point  for  a  par¬ 
ticular  mission.  In  either  case,  the  rendezvous  point  is 
determined  without  the  benefit  of  analytical  optimization. 

At  a  higher  command  level,  HQ  USAF  Studies  and 
Analysis  (SAGM)  has  created  a  refueling  model.  The  pur¬ 
pose  of  this  model  is  to  study  tanker  force  sizing  (Ref 
22) .  Given  a  particular  deployment  package,  the  number  of 
tankers  required  can  be  determined,  or  given  a  particular 
number  of  tankers,  the  most  efficient  employment  methods 
can  be  determined  and  evaluated.  This  model  is  detailed  in 
its  fuel  calculations,  and  results  have  shown  that  the  fuel 
figures  from  the  model  are  within  1  to  2  percent  of  figures 
obtained  from  the  aircraft  performance  manual  and  fuel  plan¬ 
ning  publications  (Ref  22) .  The  model  is  deterministic  and 
does  not  use  simulation  methods.  Variables  such  as  wind 
factors,  which  can  change  dynamically,  are  fixed  for  a  par¬ 
ticular  run  of  the  model.  Rendezvous  points  are  either 
manually  inserted  or  can  be  calculated  inside  the  model; 
however,  no  optimization  routine  is  used  to  try  to  reduce 
fuel  consumption  by  picking  the  best  of  all  feasible 
rendezvous  points. 

Flight  Planning 


This  section  describes  and  explains  flight  planning 
and  fuel  planning  to  include  a  discussion  of  the  terminology 


and  concepts  used  later  in  describing  th_s  research. 

Flight  planning  is  an  extremely  important  part  of  aviation, 
without  which,  intercontinental  flights  would  be  hazardous 
and  chaotic.  Since  this  study  is  concerned  with  detailed 
flight  planning  calculations,  a  brief  description  of 
flight  planning  methods  will  be  presented.  Appendix  A 
will  present  a  sample  fuel  planning  calculation  for  the 
C-141B.  These  calculations  will  be  linked  to  the  SLAM 
Simulation  Model  in  Appendix  K.  For  the  purpose  of  this 
research,  flight  planning  is  defined  as  the  art  and  science 
of  determining  a  route  of  flight  and  fuel  required  to  fly 
between  a  departure  and  destination  location.  This  plan¬ 
ning  is  done  on  the  ground  prior  to  flight.  AFR  60-16, 
General  Flight  Rules,  states  that  prior  to  each  mission 
the  pilot  in  command  will  ensure  that  the  flight  path  and 
fuel  planning  will  be  performed  in  sufficient  detail  for  a 
safe  flight  (Ref  2:p.2-l).  The  initial  phase  of  flight 
planning  determines  the  route  of  flight.  This  is  usually 
listed  on  a  computer  flight  plan  obtained  from  AF  Global 
Weather  Central  (AFGWC) ,  Offutt  AFB,  NE.  The  computer 
flight  plan  provides  flight  altitude,  route  of  flight,  dis¬ 
tances,  and  flying  times  between  reporting  points  as  well 
as  other  information.  The  flight  altitude  must  conform  to 
the  hemispherical  altitude  structure  i lir  traffic  con¬ 
trolled  airspace. 

The  next  and  most  important  part  of  flight  planning 
is  fuel  planning.  Since  fuel  planning  and  fuel  consumption 
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are  of  primary  importance  to  this  research,  fuel  planning 
will  be  decomposed  into  seven  sections  for  an  air-refueling 
mission  (see  Figure  1) . 

1.  Start,  taxi  and  takeoff. 

2.  Climb  to  an  initial  cruise  altitude. 

3.  Cruise  to  an  air-refueling  rendezvous  point. 

4.  The  refueling  maneuver. 

5.  Cruise  to  destination. 

6.  Approach  and  landing. 

7.  Holding  or  cruise  to  an  alternate  airport 
(if  applicable),  and  appropriate  fuel  reserves. 

Each  of  these  seven  areas  will  be  explained  so 
that  the  terms  will  be  familiar  when  the  air-refueling 
models  are  described.  The  following  explanation  will  be 
for  a  typical  MAC  air-refueling  mission,  and  is  not 
intended  to  be  used  as  a  guide  for  planning  an  actual 
flight . 

The  first  part  of  fuel  planning  is  start,  taxi,  and 
takeoff  (STTO) .  This  encompasses  the  fuel  required  for 
engine  start,  taxi  from  the  parking  spot,  and  acceleration 
fuel  during  the  takeoff  roll.  If  ground  delays  or  Air 
Traffic  Control  (ATC)  delays  are  known  or  anticipated,  the 
fuel  consumed  for  these  is  also  computed  and  added  to  STTO 
fuel.  The  STTO  fuel  is  usually  a  constant  fuel  quantity 
for  a  particular  aircraft.  For  example,  1900  pounds  of 
fuel  is  used  for  the  C-141B  (Ref  l:p.2-2).  The  ground 
delay  fuel  is  based  on  time  and  is  usually  computed  from  a 
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Flight  Planning  Segments 


constant  fuel  consumed  per  unit  time.  For  the  C-141B, 
the  ground  idle  fuel  consumption  rate  is  60  pounds  per 
minute.  Known  ground  delays  are  always  accounted  for  in 
fuel  planning. 

The  second  part  of  fuel  planning  is  the  climb  fuel 
consumption.  The  ground  distance  traveled  during  the  climb 
can  be  obtained  from  the  performance  manual  for  the  aircraft, 
or  from  the  computer  flight  plan.  Using  distance  and  the 
time  required  to  climb  to  cruise  altitude,  the  fuel  con¬ 
sumed  in  the  climb  can  be  determined  from  the  aircraft 
performance  manual . 

The  third  part  of  fuel  planning  is  the  fuel  required 
to  cruise  to  the  destination  or  air-refueling  rendezvous 
point.  To  calculate  this  fuel,  time  at  cruise,  cruise 
altitude,  and  aircraft  gross  weight  are  required.  Again 
tables  in  the  performance  manual  for  the  aircraft  are  used. 

Once  the  rendezvous  point  is  reached,  the  fuel 
required  for  section  four,  the  rendezvous  and  refueling 
maneuver,  is  computed.  At  this  point,  it  is  important  to 
insure  that  each  aircraft,  the  tanker  and  the  airlifter, 
has  sufficient  fuel  reserves  to  continue  their  mission. 

For  the  airlifter,  sufficient  fuel  must  be  aboard  prior  to 
the  refueling  so  that  if  unable  to  refuel,  it  can  reach 
a  suitable  abort  location  and  land  with  the  proper  fuel 
reserves.  These  reserves  will  be  detailed  later  in  this 
section.  If  additional  fuel  beyond  that  required  to 
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accomplish  the  abort  is  available,  then  e.n  air-refueling 
will  be  attempted. 

The  air-refueling  track  is  usually  a  fixed  dis¬ 
tance,  over  which  the  refueling  must  be  accomplished  (Ref 
17) .  The  track  location  and  length  are  determined  prior 
to  flight  and  the  amount  of  fuel  transferred  from  the 
tanker  will  be  enough  to  allow  the  airlifter  to  reach  its 
destination.  The  tanker,  prior  to  the  refueling,  must  have 
enough  fuel  to  supply  the  airlifter  and  return  to  its 
recovery  base  with  the  proper  fuel  reserves.  If  either  the 
tanker  or  the  airlifter  cannot  meet  these  fuel  requirements 
the  air-refueling  is  aborted.  However,  the  abort  rate  for 
airlift  aircraft  is  only  between  10  and  15  percent  (Ref  3) . 
Fuel  consumption  during  this  stage  of  flight  differs  from 
cruise  fuel  rates,  and  must  also  be  taken  into  account. 
Differences  in  fuel  consumption  rates  are  due  to  formation 
flying  with  the  refueling  boom  extended.  Detailed  fuel 
computatic  are  discussed  in  Appendix  A. 

Once  the  refueling  is  completed,  the  airlifter  and 
tanker  cruise  to  their  destinations.  Cruise  fuel  is  deter¬ 
mined  in  a  manner  similar  to  the  cruise  to  rendezvous  point 

Once  the  aircraft  approaches  the  destination,  the 
enroute  descent,  approach,  and  landing  maneuvers  occur  and 
fuel  computations  are  completed  for  section  six.  Approach 
and  landing  fuel  is  usually  given  as  a  constant  figure  for 
all  aircraft.  For  example,  the  approach  and  landing  fuel 
for  a  C-141B  is  2500  pounds  of  fuel  (Ref  l:p.2-2). 
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The  fuel  reserves  for  an  international  MAC  flight 
are  used  to  insure  that  enough  fuel  is  aboard  the  aircraft 
to  complete  the  mission  as  planned  and  to  account  for 
unforeseen  changes.  These  changes  can  occur  at  various 
times  throughout  the  mission.  EX£unples  are  ATC  delays, 
changes  in  the  direction  and  velocity  of  the  wind,  weather, 
and  changes  in  the  temperature  deviation  of  altitude. 
Combining  these  unforeseen  delays,  a  large  quantity  of  fuel 
could  be  expended.  To  ensure  that  the  mission  is  not 
jeopardized,  AFR  60-16/MAC  Supplement  One  provides  guide¬ 
lines  for  the  aircrew  establishing  procedures  to  calculate 
fuel  reserves.  These  unforeseen  delays  are  accounted  for 
in  three  ways: 

1.  Enroute  fuel  reserves. 

2.  Alternate  fuel  reserves. 

3.  Holding  fuel  reserves. 

Enroute  fuel  reserves  are  added  to  the  normal 
flight  plan  fuel  load  and  consist  of  fuel  which  is  10  per¬ 
cent  of  the  fuel  used  to  fly  over  a  category  one  route/ 
route  segment,  not  to  exceed  one  hour  fuel  at  normal 
cruise  (Ref  l:p.2-2).  A  category  one  route  is  any  route 
where  a  navigation  aid  cannot  be  flown  directly  over  once 
every  hour.  This  generally  applies  only  to  those  portions 
of  the  flight  which  are  over  water. 

The  second  fuel  reserve  is  the  fuel  required  to 
divert  to  an  alternate  airfield  due  to  weather,  airfield 
closure,  or  other  circumstances.  For  flights  outside  the 


CONUS,  an  alternate  is  always  required  (f’.f'-'  l:p.8-l). 

Fuel  reserves  to  fly  to  the  alternate  airfield  are  computed 
as  follows; 

Fuel  (is  required)  for  flight  time  from  overhead 
destination  or  initial  penetration  fix  to  alternate, 
or  to  the  most  distant  alternate  when  two  are  required, 
at  the  speed  and  altitude  in  the  appropriate  fuel  plan¬ 
ning  publication.  Add  a  ten  percent  reserve  when  time 
to  alternate  exceeds  one  hour  for  turbojet.  .  .  . 
Alternate  reserve  plus  enroute  reserve  will  not  exceed 
one  hour  at  normal  cruise  [Ref  l:p.2-2]. 

The  final  fuel  reserve  to  be  calculated  during 
flight  planning  is  the  holding  fuel.  Once  an  aircraft  has 
arrived  at  a  destination,  there  are  many  reasons  why  it 
may  be  instructed  to  hold  by  ATC.  In  any  event,  if  holding 
is  necessary,  a  holding  reserve  fuel  can  be  calculated  as 
follows:  "Holding  fuel  will  be  45  minutes  fuel  for  turbojet 
.  .  .  computed  from  the  appropriate  fuel  planning  publica¬ 
tion  using  endurance  or  holding  charts  [Ref  l:p.2-2]." 

Fuel  planning  is  completed  when  all  of  the  fuel 
figures  are  added  to  provide  a  flight  plan  fuel  load.  This 
is  the  fuel  required  to  fly  the  intended  mission.  Each 
MAJCOM  provides  planners  and  aircrews  with  fuel  management 
procedures  to  ensure  a  safe  flight.  When  the  fuel  plan¬ 
ning  is  concluded  the  important  parts  of  flight  planning 
are  accomplished. 


The  scope  of  this  research  was  limited  to  scenarios 
which  require  refueling  of  the  airlifter  to  complete  the 
mission  and  where  in-flight  refueling  is  dictated  because 
enroute  refueling  bases  are  not  available.  Missions  which 
involve  more  than  a  single  tanker  and  airlifter  aircraft 
are  not  addressed.  Two  scenarios  have  been  selected  for 
investigation.  The  first  is  an  airlifter  flight  from 
McGuire  AFB,  NJ,  to  Tehran,  Iran.  The  second  scenario  is 
a  flight  from  Travis  AFB,  CA,  to  Yokota  AB,  Japan. 

The  following  assumptions  are  made  in  this  research 

1.  The  tanker  aircraft  will  take  off  and  recover 
from  the  same  base  which  will  be  limited  to  CONUS  or  U.S. 
possession  bases.  The  recovery  restriction  stems  from  the 
high  demand  for  tankers  by  SAC,  TAC,  and  MAC  during  a 
general  contingency  operation.  Use  of  a  single  tanker 
base  insures  the  availability  of  maintenance  and  staging 
crews  for  fast  turnaround  of  the  tanker  to  support  another 
refueling  mission.  This  base  may  be  other  than  the 
tanker's  home  station  if  the  tanker  is  prepositioned  close 
to  the  refueling  point. 


in  accordance 


2.  The  aircraft  operations  must  br- 
with  AFR  60-16,  General  Flight  Rules,  and  other  major  com¬ 
mand  regulations  pertaining  to  air-refueling. 

3.  The  maximum  allowable  cargo,  limited  by  the 
maximum  gross  weight  or  volume  of  the  cargo  compartment, 
will  be  loaded  on  the  airlifter.  This  operating  practice 
will  minimize  the  total  number  of  sorties  required  to 
deliver  a  particular  deployment  package. 

Conceptualization 

Two  models  were  developed  for  this  research.  The 
first,  "FLTPLN, "  is  deterministic  and  it  models  the  air¬ 
craft  flight  planning  as  described  in  Chapter  II.  It  com¬ 
putes  the  fuel  requirements  for  a  specific  set  of  mission 
input  parameters.  The  model  compares  the  fuel  requirements 
at  65  rendezvous  points  in  the  feasible  region  and  selects 
the  geographic  point  which  minimizes  the  total  fuel  con¬ 
sumed  by  the  airlifter  and  tanker  aircraft.  The  second 
model  is  a  SLAM  simulation  model  whose  inputs  are  provided 
by  FLTPLN.  The  model  simulates  the  actual  flight  (see 
Figures  1  and  2)  of  the  two  aircraft  substituting  a  stochas¬ 
tic  variable  for  the  wind  parameter  and  adding  the  proba¬ 
bilistic  occurrence  of  takeoff  delays.  The  mean  for  the 
wind  is  the  constant  value  used  by  the  flight  planning 
model.  The  variance  of  the  wind  and  the  probability  of 
delay  are  determined  from  real  world  events.  The  simula¬ 
tion  tests  the  operational  feasibility  of  the  input 
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SCeiARIO 


Fig.  2.  In-flight  Typical  Refueling  Scenario 


parameters  verifying  the  ri  jults  of  the  FT/'TLN  model.  The 
SLAM  model  was  also  used  to  determine  the  extreme  values 
of  variance  of  the  stochastic  variables  beyond  which  the 
flight  planning  results  were  always  operationally 
infeasible.  A  detailed  description  of  the  two  models  is 
in  Appendix  K. 

The  response  variable  under  consideration  in  this 
research  is  the  total  fuel  consumed  by  an  airlifter  and 
tanker  aircraft  during  a  mission  which  requires  an 
in-flight  refueling  of  the  airlifter.  As  described  in 
Chapter  II,  the  major  factors  which  affect  the  total  fuel 
are  the  fuel  consumption  rates  for  each  aircraft  and  the 
time  that  fuel  is  consumed  at  the  respective  rates.  The 
fuel  consumption  rates  are  different  for  each  aircraft 
and  vary  for  a  specific  aircraft  with,  changes  in  aircraft 
gross  weight  and  meteorological  conditions,  such  as  tempera 
ture  and  air  density.  These  meteorological  conditions 
change  proportionally  with  the  flight  altitude  of  the  air¬ 
craft  as  well  as  randomly  with  changing  weather  patterns. 
The  gross  weight  of  the  aircraft  is  derived  from  the  basic 
airframe  weight  plus  the  cargo  and  fuel  at  takeoff.  The 
gross  weight  changes  constantly  during  flight  as  the  fuel 
weight  is  reduced.  The  flight  time  is  a  function  of  the 
total  distance  and  groui^idspeed.  The  total  distance  is  the 
flight  path  of  each  aircraft  between  its  departure  and 
destination.  Groundspeed  is  the  sum  of  the  aircraft  true 
airspeed  and  the  component  of  the  wind  velocity  in  the 


direction  of  flight.  The  true  airspeed  varies  with  the 
mach  number  and  altitude  of  the  aircraft. 

The  meteorological  conditions  are  the  most  diffi¬ 
cult  to  predict.  The  computer  flight  plans  generated  by 
AFGWC  require  a  wind  and  temperature  update  every  300  NM 
along  the  planned  flight  path  to  accurately  compute  time 
and  fuel  results.  A  worldwide  weather  data  base  updated 
in  real  time  supports  this  interaction.  Because  this 
research  is  directed  at  a  general  policy  rather  than  a 
specific  solution,  the  scenario  weather  patterns  are  held 
constant.  The  temperature  and  air  density  are  fixed  at  the 
standard  day  values  for  the  cruise  altitude  and  the  wind  is 
a  constant,  set  to  263  degrees  at  55  KNOTS.  This  wind  is 
an  approximate  mean  value  for  east-west  flights  at  mid¬ 
latitudes.  The  cruise  altitude  is  selected  at  the  highest 
Air  Traffic  Control  Hemispheric  altitude  which  the  aircraft 
can  attain  for  its  gross  weight.  Thus,  all  meteorological 
conditions  are  controlled  in  the  models. 

Each  aircraft  has  a  fixed  basic  airframe  weight 
and  operationally  assigned  cruise  mach  number.  Therefore, 
only  cargo  and  fuel  weight  are  varied  on  the  airlifter 
and  only  fuel  weight  on  the  tanker.  The  fuel  consumption 
rates  for  each  aircraft  using  these  parameters  are  obtained 
from  the  common  fuel  data  base,  referred  to  as  FLYME  (see 
Appendix  H) . 

For  an  actual  flight,  the  flight  path  may  be  modi¬ 
fied  because  of  airspace  constraints,  national  boundaries. 
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or  weather  considerations.  Since  the  tota’  distance  is  the 
only  route  factor  pertinent  to  this  research,  the  flight 
path  is  defined  as  the  great  circle  route  between  the 
departure  and  destination  points.  As  a  result,  a  scenario 
of  departure  and  destination  bases  for  each  aircraft  com¬ 
pletely  defines  the  total  distance.  The  major  factors 
which  are  varied  in  this  research  are  therefore; 

1.  the  airlifter  scenario, 

2.  the  tanker  scenario, 

3.  the  aircraft  combination, 

4.  the  airlifter 's  cargo  load,  and 

5.  the  airlift  initial  fuel  load. 

In  order  to  investigate  the  interactive  effect  of  each 
factor,  a  full  factorial  design  is  employed  combining  all 
levels  of  each  factor  with  the  levels  of  all  other  factors. 
This  design  varies  the  levels  of  only  one  factor  at  a  time 
while  keeping  the  others  constant.  This  routine  is  repeated 
until  all  levels  of  all  factors  are  examined.  Table  2 
lists  the  design  factors  and  their  levels. 

Two  airlifter  scenarios  were  investigated,  one 
east  bound  with  a  tail  wind  and  one  west  bound  with  a  head 
wind.  The  attempt  was  to  offset  the  effects  of  holding  the 
wind  constant  in  the  model.  For  each  airlifter  scenario, 
four  tanker  departure  bases  were  considered  to  test  the 
impact  of  tanker  basing  on  the  optimal  rendezvous  point. 

The  tanker  bases  were  selected  with  one  near  the  airlifter 's 
departure  base,  one  to  the  north,  one  to  the  south,  and  the 


TABLE  2 


EXPERIMENTAL  DESIGN  FACTORS  AND  LEVELS 


i;  Factor 

- 


Level 


Airlifter  Scenario  a) 

b) 

Tanker  Scenario  a) 

b) 

c) 

d) 

e) 

f) 

Aircraft  Cci±)ination  a) 

b) 

c) 

d) 


McGuire  AFB,  NJ  to  Tehran,  Iran 
Travis  AFB,  CA  to  Yokota  AB,  Japan 

Castle  AFB,  CA 
Charleston  AFB,  SC 
Eielson  AFB,  AL 
Fairchild  AFB,  VJA 
Loring  AFB,  ME 
McGuire  AFB,  NJ 

C-5/KC-10 

C-5/KC-135 

C-141B/KC-10 

C-141B/KC-135 


Airlift  Cargo  Load 


a) 

1) 

100, 

,000 

lbs 

2) 

85, 

,000 

lbs 

3) 

70, 

,000 

lbs 

b) 

C- 

141B 

1) 

70, 

,000 

lbs 

2) 

55, 

,000 

lbs 

3) 

40, 

,000 

lbs 

Airlifter  rtexiraum  Fuel  Load 


Fuel  (1000  Its.) _ 

Cargo  Cargo 


C-5 

(1000  lbs.) 

C-141B 

(1000  lbs.) 

261.3 

100.0 

114.0 

70.0 

276.3 

85.0 

129.0 

55.0 

291.3 

70.0 

144.0 

40.0 

NOTE:  These  initial  maximum  fuel  loads  were 
decremented  by  20,000  lbs.  ten  times  or  until  infeasible 
fuel  loads  were  obtained. 
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fourth  was  Barksdale  AFB,  LA.  Barksdale  u  the  proposed 
base  of  assignment  for  the  KC-10  (Ref  23) .  This  allows 
comparison  of  the  KC-10  operating  from  Barksdale  against 
operating  from  a  base  closer  to  the  airlifter  departure 
base.  The  bases  considered  in  each  scenario  are  listed  in 
Table  2. 

For  each  scenario  and  aircraft  combination,  many 
rendezvous  points  are  feasible.  Rather  than  compare  all 
of  the  feasible  points,  a  small  rectangular  region  is  con¬ 
sidered.  The  rectangle  was  fixed  along  the  feasible  boun¬ 
dary  furthest  from  the  airlifter  destination.  The  size  of 
the  rectangle  is  20  degrees  of  latitude  by  24  degrees  of 
longitude.  The  region  was  determined  to  always  contain 
the  optimal  rendezvous  point.  All  other  points  have  higher 
total  fuel  values  and  the  values  increase  monotonically  in 
all  directions  away  from  the  optimal  point.  As  a  result, 
only  the  matrix  of  points  in  the  rectangle  were  required 
to  define  the  optimal  point.  For  rendezvous  points  near 
the  airlifter 's  departure  base,  the  leveloff  point  after 
initial  climbout  was  considered  as  the  first  feasible 
rendezvous  point.  These  cases  are  close  to  not  requiring 
an  air-refueling  and  are  treated  as  extreme  levels.  A 
single  abort  base  is  designated  for  each  scenario.  In  the 
case  of  an  aborted  refueling,  the  airlifter  must  have  suf¬ 
ficient  fuel  to  recover  to  the  departure  base  or  the 
abort  alternate,  whichever  is  closer.  Loring  AFB,  ME,  is 
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the  abort  alternate  for  Lh-^  east  scenario  ;  nd  Elmendorf  AFB, 
AK,  is  the  abort  alternate  for  the  west  scenario. 

The  structure  of  the  SLAM  model  exactly  follows 
the  sequence  of  flight  planning  steps  outlined  in  Chapter  II. 
The  wind  velocity  is  allowed  to  vary  and  delays  which  result 
in  additional  fuel  consumption  are  included.  The  delay 
data  is  based  on  the  actual  home  station  delay  rates  from 
Travis  AFB,  CA,  and  Dover  AFB,  DE,  for  the  C-5A,  Norton 
AFB,  CA,  and  Travis  AFB,  CA  for  the  C-141,  and  March  AFB, 

CA,  and  Plattsburg  AFB,  NY,  for  the  KC-135.  The  KC-10  is 
not  operational,  so  the  results  of  the  C-5A  data  were  used 
for  the  KC-10  also.  The  raw  sample  data  for  these  delay 
rates  are  in  Appendix  B. 

Multiple  runs  of  each  scenario  are  made  in  the  SLAM 
model  to  test  the  feasibility  of  the  FLTPLN  results.  The 
sample  size,  or  number  of  replications  per  flight  was 
determined  while  the  sample  runs  were  conducted.  The 
samp]»e  size  was  directly  determined  from  results  generated 
by  the  SLAM  model  output.  Initially,  runs  of  ten  replica¬ 
tions  with  a  variance  reduction  technique  were  used.  The 
variance  reduction  technique  was  antithetic  variates 
(Ref  18:484-485).  Using  antithetic  variates  and  twenty 
replications  reduced  the  variance  obtained  from  twenty 
replications  using  regular  Monte  Carlo  simulation  tech¬ 
niques  in  every  case  tested. 
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IV.  Verification  and  Validation 

The  steps  taken  to  validate  the  analytical  flight 
planning  model,  FLTPLN,  are  described  in  this  chapter. 

The  validation  included  verification  of  three  aspects  of 
the  model;  input  parameters,  internal  design,  and  tech¬ 
nical  accuracy  of  the  computations.  The  validity  of  the 
output  of  the  model  was  also  measured  against  the  objec¬ 
tives  and  assumptions  used  in  the  modeling  efforts.  The 
validity  criteria  was  set  at  5  percent  deviation  from  the 
performance  charts  for  the  FLYME  fuel  data  and  10  percent 
deviation  of  the  total  fuel  consumption  results  of  the 
model  from  manual  computations.  Final  validation  of  FLTPLN 
for  operational  suitability  was  accomplished  by  the  SLAM 
simulation  model. 

The  input  parameters  used  by  FLTPLN  are  constants 
or  derived  results  computed  in  other  model  subroutines. 

The  constants  are  either  levels  of  the  factors  varied  in 
the  experiment  or  parameters  extracted  directly  from 
the  respective  aircraft  performance  manual.  The  use  of 
operational  constants  such  as  cruise  mach  number  signifi¬ 
cantly  adds  to  the  external  validity  of  the  model  because 
these  are  the  values  in  actual  use  by  the  MAJCOMs.  The 
derived  results  include  fuel  data  from  the  FLYME  set  of 
subroutines,  great  circle  distance  and  course  computations 
from  subroutine  RHOTHTA  ( see  Appendix  H) ,  and  ground  speed 
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computations  from  within  FLTPLN  (see  Appcniix  C) .  The 
FLYME  fuel  data  was  technically  verified  by  comparing 
selected  values  to  the  aircraft  performance  charts.  The 
results  are  presented  in  Table  3.  The  average  deviation 
of  the  FLYME  data  from  the  performance  charts  as  a  per¬ 
centage  of  the  chart  values  is  4  percent.  This  deviation 
was  consistent  for  all  parameters  in  both  the  FLTPLN  and 
SLAM  models. 

Appendix  A  contains  a  sample  fuel  consumption  cal¬ 
culation  of  a  typical  C-141  mission  using  the  aircraft  per¬ 
formance  manuals.  Table  4  compares  these  manual  fuel  con¬ 
sumption  values  to  the  manual  computations.  These  are 
sufficiently  close  to  verify  that  the  aircraft  fuel  con¬ 
sumption  calculations  are  representative  of  the  respective 
aircraft . 

The  internal  design  of  FLTPLN  was  verified  by  pro¬ 
ceeding  manually  through  each  step  of  the  program  to  ensure 
that  values  were  used  as  the  designers  intended  and  not 
changed  or  internally  lost.  The  output  parameters  were 
also  cross-checked  to  verify  one  result  from  another. 

For  example,  the  onload  fuels  must  equal  the  sum  of  the 
route  fuels  plus  reserves.  From  these  tests,  the  program's 
design  was  technically  certified  as  accurately  computing 
and  printing  the  experimental  results. 

The  output  of  the  models  was  displayed  in  a 
5  X 13  matrix  consisting  of  total  fuel  consumption  values 
for  a  specific  rendezvous  point.  For  every  scenario  and 
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TABLE  3 


VALIDATION  OF  FLYME  DATA  BASE 


Input 

Parameters 

Clint) 

Time 

(Min) 

Climb 

Dist 

(NM) 

Clint) 

Fuel 

(1000  Ib) 

Cruise 

Fuel 

Rate 

(1000/Min) 

Holding 

Fuel 

Rate 

(1000/Min) 

C-5 

FLYME 

22 

133 

12.3 

.32 

.20 

500,000# 

CHT 

21 

129 

11.7 

.32 

.23 

FL350 

%  Dev 

+  .05 

+.03 

+  .05 

0 

-.12 

C-5 

FLYME 

25 

150 

16.1 

.42 

.23 

650,000# 

CHT 

25 

150 

15.3 

.42 

.26 

FL290 

%  Dev 

0 

0 

+.05 

0 

-.12 

C-141B 

FLYME 

15.4 

93.6 

6.6 

.20 

.16 

240,000# 

CHT 

15.1 

91.0 

6.3 

.21 

.16 

FL350 

%  Dev 

+  .02 

+.03 

+  .05 

-.05 

0 

C-141B 

FLYME 

15.8 

92 

7.6 

.25 

.19 

300,000# 

CHT 

14.8 

90.5 

6.9 

.26 

.19 

FL290 

%  Dev 

+.07 

+  .02 

+.10 

-.04 

0 

KC-10 

FLYME 

18.3 

127.5 

10.3 

.35 

.24 

420,000# 

CHT 

17.6 

121 

9.8 

.35 

.24 

FL370 

%  Dev 

+  .04 

+  .05 

+.05 

0 

0 

KC-10 

FLYME 

22.2 

154.5 

13.7 

.41 

.31 

575,000# 

CHT 

21.2 

151 

13.2 

.39 

.31 

FL310 

%  P^V 

+.05 

+.02 

+.04 

+.05 

0 

KC-135 

FLYME 

19.3 

128.6 

6.6 

.18 

.16 

180,000# 

CHT 

20.0 

127.0 

6.6 

.18 

.16 

FL350 

%  Dev 

-.04 

+.01 

0 

0 

0 

KC-135 

FLYME 

27.2 

171.5 

9.9 

.25 

.20 

250,000# 

CHT 

28.0 

171.0 

9.9 

.25 

.20 

FL290 

%  Dev 

-.03 

.00 

0 

0 

0 

NOTE: 

CHT  values  are 

obtained 

from  aircraft 

perform- 

ance  manuals. 
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combination  of  factors,  FLTPLN  determined  an  optimal 
rendezvous  point  for  which  the  total  fuel  requirement  was 
lower  than  any  other  point  considered.  This  can  be  seen 
in  the  FLTPLN  output  in  Appendix  G.  While  this  matrix  repre¬ 
sented  only  a  portion  of  the  feasible  region  of  points, 
the  optimal  point  was  always  surrounded  by  higher  values. 
Every  matrix  of  total  fuel  values  had  a  gradient  increasing 
monotonically  from  the  optimal  point  out  to  the  feasible 
boundary  or  to  the  edge  of  the  matrix.  Thus  the  minimum 
value  of  the  matrices  is  considered  the  absolute  minimum 
for  the  entire  feasible  region. 

The  output  was  further  validated  by  checking  that 
changes  in  the  input  parameters  resulted  in  consistent 
changes  in  the  output.  For  example,  higher  cargo  weights 
or  route  distances  resulted  in  higher  total  fuel  consump¬ 
tion  . 

Further  validation  of  FLTPLN  results  was  accom¬ 
plished  using  the  SLAM  model.  FLTPLN  output  (total  fuel 
consumption  per  aircraft)  was  compared  to  SLAM  output  and 
mission  success  rates  were  computed.  In  this  way,  FLTPLN 
output  could  be  shown  to  be  feasible  under  actual  real-time 
flight  conditions.  The  SLAM  validation  of  FLTPLN  was  done 
for  each  aircraft.  The  results  of  this  validation  are 
described  in  the  following  paragraphs. 

Each  of  the  twenty-three  C-5  flights  simulated  by 


the  SLAM  model  was  successful  under  the  normal  wind  factor 
using  probabilistic  ground  delays  of  zero  or  fifteen 


minutes.  The  SLAM  fuel  consumption  values  averaged  2.1 
percent  higher  than  those  calculated  by  the  analytical 
model.  Deviation  of  the  SLAM  fuel  consumptions  for  the 
C-5  are  shown  in  Figure  3  for  particular  cargo  loads, 
tanker  base,  and  refueling  tankers.  A  negative  value  means 
the  SLAM  model  computed  less  fuel  than  the  analytical  model 
and  a  positive  value  means  the  SLAM  model  computed  more 
fuel.  Even  with  this  increased  fuel  consumption,  the  C-5 
was  still  able  to  complete  all  flights  successfully.  To 
further  validate  FLTPLN  a  95  percent  confidence  interval 
based  on  the  SLAM  fuel  consumption  mean  and  variance  was 
constructed  for  each  aircraft.  Using  the  upper  limit  of 
the  confidence  interval  as  the  actual  C-5  fuel  consumption, 
and  as  a  Mission  completion  success  criteria,  87  percent 
of  the  flights  were  successful.  For  example,  if  the  avail¬ 
able  fuel  for  the  C-5  exceeded  the  upper  limit  of  the  fuel 
consumption  confidence  interval  the  mission  was  successful. 

Every  one  of  the  twenty-three  C-141B  flights  simu¬ 
lated  by  the  SLAM  model  was  successful.  The  SLAM  fuel  con¬ 
sumption  values  averaged  2.7  percent  higher  (see  Figure  4) 
than  that  calculated  by  the  analytic  model.  Using  the 
upper  limit  of  the  95  percent  confidence  interval  as  the 
actual  fuel  consumed  in  flight  and  as  a  mission  success 
criteria,  83  percent  of  the  C-141  flights  were  successful. 
Combining  the  success  rates  for  both  aircraft  the  analytic 
model  correctly  predicted  the  proper  ramp  fuel  load  for 
the  airlift  aircraft. 
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3.  C-5 — Deviation  of  SLAM  Results  Compared  to 

the  Analytical  Fuel  Consumption  Results 
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4.  C-141 — Deviation  of  SLAM  Results  Compared 

to  the  Analytical  Fuel  Consumption  Results 


One  hundred  percent  of  the  SLAM  tcrker  flights 

were  successful  using  the  SLAM  wind  variant  and  ground 

delays.  Also  each  of  the  flights  was  successful  using  the 

upper  limit  of  the  95  percent  confidence  interval  as  the 

success  criteria.  The  SLAM  KC-10  fuel  consumption  values 

averaged  less  than  a  1  pfercent  decrease  (see  Figure  5) 

from  the  analytical  model  result  and  the  KC-135  averaged 
\ 

less  than  a  2  percent  decrease  (see  Figure  6)  from  the 
analytical  model  result.  From  these  deviations  we  con¬ 
cluded  that  the  analytical  model  provides  results  that  are 
operationally  feasible  under  actual  flight  conditions  for 
all  four  aircraft. 

Overall,  the  tankers  burned  about  1  percent  less 
fuel  in  the  SLAM  model  than  in  the  analytical  model.  Con¬ 
versely,  the  airliftcrs  burned  about  2.5  percent  more  fuel 
in  the  SLAI4  model  than  in  the  analytical  model.  This 
difference  in  airlift  fuel  can  be  explained  from  the  nature 
of  the  two  aircraft  missions.  The  tanker  flies  to  a 
rendezvous  point  and  returns  to  his  departure  base  which 
is  relatively  close  compared  to  the  airlifter  destination; 
the  airlifter,  however,  flies  to  the  rendezvous  point  and 
then  continues  on  to  a  destination  which  is  thousands  of 
miles  longer  than  the  tanker  route  of  flight.  Conse¬ 
quently,  the  airlifter  has  more  time  to  be  affected  by 
increases  in  the  wind  factor,  and  these  increases  plus  any 
delay  in  fuel  consumption  would  cause  the  increased  fuel 
requirement.  There  are  two  reasons  the  tanker  SLAM  results 
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Fig.  5.  KC-10--Deviation  of  SLTVM  Results  Compared 
to  the  Analytical  Fuel  Consumption  Results 


are  less  than  the  analytical  consumption.  First  is  that 
an  increase  in  wind  factor  has  less  effect  on  the  tanker 
than  on  the  airlifter.  A  head  wind  increase  to  the 
rendezvous  point  becomes  a  tail  wind  increase  when  return¬ 
ing  to  destination.  Second,  is  the  difference  in  the  loiter 
fuel  computations.  In  FLTPLN  the  tanker  loiter  fuel  is 
always  based  on  a  fifteen-minute  loiter  duration;  however, 
in  SLAM  the  loiter  time  is  based  on  actual  aircraft 
arrivals  which  usually  requires  less  than  fifteen  minutes 
loiter  time.  For  the  tankers,  fifteen  minutes  loiter  fuel 
is  more  than  the  difference  between  SLAM  and  FLTPLN  fuel 
values.  However,  without  any  adjustments  for  loiter  fuel 
the  fuel  consumption  values  remain  within  our  10  percent 
validation  criteria.  Finally,  from  these  validation 
results  we  ooserve  that  the  SLAM  model  has  validated  the 
results  of  the  analytical  model.  The  SLAM  model  validation 
is  discussed  in  Appendix  I. 
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V.  Results  and  Analysis 

The  major  purpose  in  the  analysis  of  the  research 
results  was  to  identify  the  principal  factors  affecting 
the  determination  of  the  optimal  rendezvous  point  and 
takeoff  fuel  loads.  First,  these  principal  factors  were 
analyzed  in  terms  of  the  two  hypotheses  stated  in  Chapter  I. 

The  effect  of  cargo  load  on  the  optimal  rendezvous 
point  is  the  first  factor  to  be  considered.  The  first 
hypothesis  stated  that  the  optimal  rendezvous  point  would 
occur  at  the  maximum  flight  range  of  the  airlifter  from 
its  destination  base.  This  maximum  range  from  the  destina¬ 
tion  constitutes  a  boundary  of  the  region  of  feasible 
rendezvous  points  closest  to  the  airlifter  takeoff  base. 

The  dashed  lines  on  Figures  7  to  10  represent  this  boundary 
at  the  levels  of  the  cargo  weight  indicated.  The  optimal 
rendezvous  points  determined  by  FLTPLN  did  not  always  occur 
at  this  boundairy  range.  The  separation  between  the  ren¬ 
dezvous  points  and  the  feasible  boundary  is  indicated  on 
the  figures  by  the  arrows.  For  example,  in  Figure  7,  the 
100,000  pounds  cargo  load  on  the  C-5A  resulted  in  an 
optimal  rendezvous  point  location  244  nautical  miles  from 
the  respective  feasible  boundary.  The  results  for  the 
C-141  on  Figure  8  are  similar.  As  the  cargo  weight  is 
decreased,  the  range  of  the  airlifter  increases.  This 
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shifts  the  feasible  boundary  away  from  the  airlifter 
destination  base.  For  each  decrease  in  cargo  weight, 
the  optimal  rendezvous  points  shift  toward  the  airlifter 
takeoff  base  along  with  the  corresponding  feasible  boun¬ 
dary.  However,  the  optimal  rendezvous  point  is  located 
farther  from  the  hypothesized  location  with  each  decrement 
For  the  maximum  feasible  cargo  weight  of  an  airlifter,  the 
optimal  point  occurs  on  the  feasible  boundary,  as  is 
indicated  on  Figure  9  for  70,000  pounds.  Beyond  these 
cargo  weights,  the  shorter  range  of  the  airlifter  does 
not  permit  a  feasible  solution.  These  results  contradict 
the  first  hypothesis. 

The  second  factor  to  be  considered  is  tanker  base 
selection.  Five  of  the  six  air lifter/cargo  weight  combina¬ 
tions  for  the  east  scenario  resulted  in  the  same  optimal 
rendezvous  points  regardless  of  which  tanker  base  was  used 
However,  for  the  west  scenario,  most  tanker  bases  resulted 
in  a  unique  rendezvous  point  for  a  given  airlif ter/cargo 
weight  combination.  The  rendezvous  point  using  Eielson 
was  different  from  the  rendezvous  point  using  Fairchild 
or  Castle.  The  location  of  Eielson  far  along  the  western 
route  of  the  airlifter,  precipitated  this  result.  However 
each  of  the  bases  had  some  influence  on  the  rendezvous 
point  measured  by  the  distance  of  the  point  from  the  great 
circle  route  of  the  airlifter.  Table  5  lists  these  dis¬ 
tances  by  bases.  The  influence  of  the  tanker  base  con¬ 
sistently  increases  as  the  bases  are  located  farther  along 
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TABLE  5 


TANKER  BASE  INFLUENCE  ON  RENDEZVOUS  POINT 


Base 

Location 

Eielson 

758.7  NM 

Fairchild 

402.2  NM 

Castle 

93.7  NM 

the  airlifter's  route.  The  lack  of  variance  due  to  tanker 
basing  in  the  east  scenario  rendezvous  points  stems  from 
the  increased  fuel  requirements  for  the  airlifter.  Because 
of  the  large  fuel  requirement  to  reach  the  ultimate  des¬ 
tination  in  the  eastern  scenario,  little  flexibility  was 
present  in  determining  the  rendezvous  point.  Thus,  the 
location  of  the  tanker  base  affects  the  selection  of  the 
optimal  rendezvous  point  more  when  the  fuel  constraints 
on  the  airlifter  are  relaxed  by  enroute  tanker  basing  or 
shorter  total  distances  for  the  airlifter. 

The  second  hypothesis  stated  that  the  maximum 
allowable  takeoff  fuel  load  for  the  airlifter  would  result 
in  less  total  fuel  consumption  than  any  other  fuel  load. 
However,  for  most  scenarios,  the  optimal  airlifter  fuel 
load  was  the  minimum  required  to  reach  the  optimal 
rendezvous  point  and  still  maintain  sufficient  fuel 
reserves  to  abort  to  the  closest  recovery  base.  This  is 
indicated  by  the  lowest  point  on  the  graph  in  Figure  11. 
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Fig.  11.  Total  Fuel  Consumption  versus 
Initial  Airlifter  Fuel  Load 


Investigating  higher  fuel  onloads  for  the  se.me  scenario, 
the  optimal  rendezvous  point  was  found  to  remain  the  same, 
but  the  total  fuel  requirement  increased.  This  trend  can 
be  seen  in  the  sample  FLTPLN  output  data  in  Appendix  B. 

For  smaller  fuel  loads,  the  rendezvous  point  was  located 
closer  to  the  airlifter  takeoff  base.  This,  however, 
resulted  in  a  larger  total  fuel  requirement  for  the  mission. 

This  trend  is  also  shown  by  the  graphs  in  Figures 
12  to  15.  The  fuel  loads  of  the  airlifter  and  tanker  com¬ 
binations  are  expressed  as  percentages  of  their  respective 
fuel  capacities.  Each  curve  corresponds  to  a  different 
combination  of  cargo  weight  and  tanker  base.  As  the 
airlifter  fuel  load  is  reduced,  the  tanker  fuel  load  is 
increased  by  an  amount  proportional  to  the  tanker  fuel 
decrease.  The  total  fuel  represented  between  any  single 
point  on  a  curve  (Figures  12  to  15)  varies  only  a  few 
thousand  pounds.  Therefore,  the  transfer  fuel  is  the  only 
factor  which  changes  significantly.  The  slope  of  each  line 
represents  the  relative  efficiency  of  the  aircraft  for 
transporting  fuel  to  the  rendezvous  point.  Similarities 
in  the  C-5/KC-10  graphs  (Figure  12)  and  the  C-141B/KC-135 
(Figure  13)  graphs  can  be  explained  by  comparing  fuel 
consumption  rates.  For  example,  the  C-5  and  KC-10  have 
similar  fuel  consumption  rates,  as  do  the  C-141B  and  KC- 
135.  The  slopes  of  the  lines  for  these  combinations  are 
approximately  -1.2.  The  approximate  slopes  for  the  C-5/ 
KC-135  (Figure  14)  and  C-141/KC-10  (Figure  15)  which  have 
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vastly  different  fuel  consumptions,  are  -0.6  and  -2.6 
respectively.  Therefore,  for  the  latter  two  combinations, 
the  C-5  and  KC-10  utilize  a  smaller  percentage  of  their 
fuel  capacity  than  the  KC-135  or  C-141B  for  the  same  total 
amount  of  fuel. 

The  optimal  fuel  loads  marked  by  a  box  on  each 
curve,  corresponds  to  the  two  percentages  whose  sum  is  the 
smallest  of  any  two  fuel  load  percentages  on  a  curve. 

For  example,  the  McGuire  70,000  pound  line  in  Figure  12 
has  the  optimal  load  at  the  point  corresponding  to  a  Siam 
of  112  percent  for  both  aircraft  while  other  points  on  the 
graph  increase  up  to  123  percent  at  the  extreme  end  of  the 
curve. 

The  optimal  airlifter  fuel  load  did  not  always 
occur  at  the  minimum  level.  This  is  indicated  by  the 
boxes  on  Figures  13  and  15.  For  these  curves,  the  maximum 
airlifter  fuel  load  was  optimal.  However,  because  of  the 
differences  in  slopes  for  these  combinations,  the  sum  of 
the  percentages  is  still  the  minimum  for  all  values  along 
the  curve.  This  minimum  sum  criterion  contradicts  hypo¬ 
thesis  two.  The  relative  efficiencies  of  the  two  aircraft 
to  carry  the  fuel  to  the  rendezvous  point,  represented  by 
the  slope  of  the  curves  in  Figures  12  to  15,  are  the  deter' 
mining  considerations  for  optimal  takeoff  fuel  loads. 

The  main  result  of  this  research  is  the  identifica' 
tion  of  significant  fuel  savings  which  can  be  derived  from 
the  use  of  the  optimal  rendezvous  point  and  takeoff  fuel 
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loads.  The  exact  savings  depend  on  the  scenario  and  air¬ 
craft  combination.  Table  6  lists  the  average  fuel  savings 
by  aircraft  combination  for  the  optimal  rendezvous  point 
and  fuel  loads  compared  to  the  point  and  fuel  loads  defined 
by  the  two  hypotheses. 


TABLE  6 

AVERAGE  FUEL  SAVINGS —OPTIMAL  RENDEZVOUS  POINT 
VERSUS  HYPOTHESIZED  RENDEZVOUS  POINT 


Scenario 

Aircraft 

Combination 

Average  Fuel 
Savings 
(lbs . ) 

East 

C-5/KC-10 

25,200 

C-141B/KC-10 

8,300 

C-5/KC-135 

19,700 

C-141B/KC-135 

6,300 

West 

C-5/KC-10 

39,400 

C-141B/KC-10 

21,700 

C-5/KC-135 

34,700 

C-141B/KC-135 

10,300 

Since  the  original  two  hypotheses  were  contra¬ 
dicted  by  the  research  results,  the  hoped  for  general 
optimal  refueling  policy  was  not  derived.  The  results 
were  examined  to  see  if  any  general  hypotheses  were  sug¬ 
gested. 

The  results  are  summarized  below  by  describing  the 
effect  of  the  principal  factors  on  the  determination  of 
the  optimal  rendezvous  point  and  onload  fuels.  At  the 
maximum  feasible  cargo  weights,  the  optimal  rendezvous 
point  occurs  on  the  feasible  boundary  closest  to  the 
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airlift  takeoff  base.  For  lesser  cargo  v.-eights,  the 
separation  between  the  point  and  the  coresponding  feasible 
boundary  increases.  Enroute  tanker  bases  draw  the  loca¬ 
tion  of  the  optimal  rendezvous  point  closer  to  the  tanker 
base  than  do  inland  bases.  The  inland  bases  will  most 
often  result  in  the  same  optimal  rendezvous  point  despite 
vast  differences  in  the  tanker's  flight  distances.  The 
takeoff  fuel  loads  for  both  aircraft  are  optimized  by  mini¬ 
mizing  the  combined  percentage  of  fuel  capacity  used  by  the 
two  aircraft.  This  sum  is  greatly  dependent  on  the  spe¬ 
cific  aircraft  combination  used. 

Sensitivity  Analysis 

The  sensitiyity  analysis  systematically  yaried  the 
levels  of  selected  input  factors  and  was  used  to  verify 
the  consistency  of  the  results  of  the  models  for  the  entire 
feasible  range  of  the  input  variables.  The  analytic 
results  were  collected  for  extreme  values  of  airlifter  and 
tanker  total  flight  distances,  and  cargo  weights  for  each 
aircraft  combination.  Sensitivity  analysis  of  the  SLAM 
model  was  completed  for  the  extreme  values  of  ground 
delays,  wind  factor,  and  increased  in-flight  cruise  times. 

The  sensitivity  of  the  results  to  the  total  dis¬ 
tance  flown  by  the  airlifter  was  tested  by  generating 
scenarios  with  different  destinations.  For  the  east 
scenario  CONUS  bases,  an  airlifter  destination  at  Cairo, 
Egypt,  represented  the  minimum  distance  examined.  This 
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reduced  the  Tehran  flight  distance  by  570  ”M.  The  maximum 
extreme  distance  examined  used  Kabal,  Afghanistan,  as  a 
destination,  which  increased  the  airlifter  flight  distances 
by  510  NM.  The  Cairo  data  results  were  consistent  with 
the  Tehran  data  except  previously  unfeasible  cargo  loads 
became  feasible.  For  example,  120,000  pounds  cargo  levels 
for  the  C-5A  and  the  70,000  pounds  cargo  level  for  the 
C-141B  both  became  feasible.  For  the  Kabal  data,  the 
55,000  pounds  cargo  load  for  the  C-141B  became  infeasible. 
As  a  result,  few  airlifter  fuel  loads  were  feasible.  The 
rendezvous  points  were  located  closer  to  the  airlifter  des¬ 
tination  and  the  transfer  fuel  at  this  increased  tanker 
range  was  reduced  by  an  average  of  23,000  pounds.  These 
results  are  continuations  of  the  trends  established  by  the 
main  data  and  are  consistent  with  the  previous  results. 

A  southwest  scenario  from  Travis  AFB,  CA,  to  Clark 
AB,  Philippines,  provided  an  example  of  enroute  basing 
using  Hickam  AFB,  HI,  as  the  tanker  base.  The  rendezvous 
point  moved  an  average  of  1203  NM  from  the  great  circle 
route  for  all  runs  of  this  scenario.  The  use  of  Barksdale 
AFB,  LA,  by  the  KC-10  resulted  in  the  same  rendezvous  point 
as  Castle,  though  the  extra  2500  NM  route  segment  for  the 
tanker  increased  the  total  fuel  used.  Thus,  these  extreme 
examples  of  enroute  and  deep  inland  tanker  basing  were  con¬ 
sistent  with  previous  results  in  the  determination  of  the 
optimal  rendezvous  point. 


The  effects  of  cargo  loads  on  the  results  of  the 


analytic  model  were  consistent  with  the  Tehran  results  for 
all  feasible  values-  For  large  cargo  loads  which  pre¬ 
cluded  transferring  sufficient  fuel  to  complete  the  mission, 
no  feasible  solution  was  possible.  This  occurred  at  140,000 
pounds  for  the  C-5A  and  70,000  pounds  for  the  C-141B.  For 
small  cargo  loads,  the  allowable  fuel  load  was  sufficient 
to  conduct  the  mission  without  refueling.  Between  these 
extremes,  changes  in  the  cargo  load  consistently  produced 
the  previous  results  illustrated  in  Figures  7  to  10. 
Therefore,  the  entire  feasible  ranges  of  the  input  param¬ 
eters  produced  consistent  results  for  all  scenarios. 

Extreme  values  of  the  factors  indicate  feasible  boundaries 
but  do  not  change  the  trends  or  results. 

Although  the  primary  purpose  of  the  SLAM  model  was 
to  validate  the  analytical  model  for  operational  suitabil¬ 
ity,  it  was  also  used  to  conduct  part  of  the  sensitivity 
analysis.  The  SLAM  sensitivity  analysis  was  conducted  in 
three  major  areas;  (1)  increases  in  ground  delay  times, 

(2)  increases  in  wind  factors,  and  (3)  increases  in 
cruise  times. 

The  first  variable  investigated  was  the  ground 
delays.  In  the  SLAM  model,  the  delay  was  either  zero,  or 
fifteen  minutes.  To  determine  the  sensitivity  of  ground 
delay  time  on  the  results,  the  delay  was  increased  to 
twenty-five  minutes.  Increased  ground  fuel  consumption 
resulted;  for  example,  800  pounds  for  the  C-141B,  1200 
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pounds  for  the  C-5,  but  this  had  no  effect  on  the  mission 
completion  rates  for  any  of  the  aircraft  or  scenarios.  A 
larger  increase  in  the  ground  delay  time  was  not  explored 
because  from  the  experience  of  the  authors,  an  aircraft 
commander  would  not  run  the  engines  for  longer  than  twenty- 
five  minutes  on  the  ground  while  waiting  for  maintenance 
repairs  or  other  types  of  services. 

The  second  variable  investigated  was  the  wind 
factor.  The  SLAM  model  was  created  with  a  normally  dis¬ 
tributed  wind  variant  with  a  mean  of  55  KNOTS  and  a  standard 
deviation  of  10  KNOTS.  The  wind  was  first  changed  to  a 
mean  of  55  KNOTS  and  standard  deviation  of  20  KNOTS.  This 
caused  43  percent  of  the  missions  to  fail  due  to  inadequate 
fuel.  When  the  wind  standard  deviation  was  increased  to 
30  KNOTS  none  of  the  missions  flown  against  a  headwind 
were  successful.  It  is  not  surprising  that  the  missions 
failed  when  the  wind  standard  deviation  was  increased  to  20 
KNOTS  since  that  is  a  significant  increase.  Those  missions 
which  failed,  all  flew  the  west  scenario  to  Yokota,  had  a 
headwind  for  most  of  the  route.  The  effect  of  an  increase 
of  20  KNOTS  on  a  C-141  which  travels  3000  NM  after  the 
rendezvous  (this  is  representative  of  the  rendezvous  point 
to  destination  distance  for  the  west  scenario)  increased 
the  fuel  consumption  by  4400  pounds.  A  wind  increase  of 
30  KNOTS  increased  the  fuel  consumption  by  6400  pounds 
which  caused  all  C-141B  missions  to  fail.  In  actual 
flight  operations  large  wind  factor  deviations  are  not 


expected.  Combining  weather  satellite  data  and  pilot  mis¬ 
sion  weather  reports  which  are  forwarded  to  Air  Force 
Global  Weather  Central,  the  current  wind  forecasts  are  very 
accurate . 

The  third  variable  investigated  was  the  increase  in 
unplanned  cruise  time  for  each  aircraft.  To  examine  these 
delays,  the  cruise  time  was  increased  by  fifteen  minutes 
for  each  aircraft.  With  the  extra  cruise  time  added,  83 
percent  of  the  C-5  missions  were  completed  but  for  the 
C-141  only  74  percent  of  the  missions  were  completed.  For 
the  KC-10,  100  percent  of  their  missions  were  successful 
but  for  the  KC-135  only  83  percent  of  their  missions  suc¬ 
ceeded.  The  low  mission  completion  rates  seem  to  indicate 
that  a  fuel  reserve  is  needed  to  make  the  mission  comple¬ 
tion  rate  acceptable.  HQ  MAC  does  have  a  fuel  reserve 
for  aircraft  flying  on  overwater  routes.  If  these  fuel 
reserves  are  added  to  each  aircraft's  fuel  load,  the  mis¬ 
sion  completion  rates  increase  to  100  percent. 

Since  the  original  fuel  load  hypothesis  has  been 
rejected  and  the  optimal  airlifter  ramp  fuel  loads  are  most 
often  less  than  the  maximum  fuel  allowable,  adequate  fuel 
capacity  remains  for  a  fuel  reserve.  The  reserve  would  com¬ 
pensate  for  cruise  delays,  and  should  not  change  the  trends. 

The  sensitivity  analysis  performed  indicates  that 
variations  in  the  parameters  of  the  FLTPLN  model  used  in 
the  specific  scenarios  of  this  research  do  not  change  the 
basic  results  obtained. 
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I.  Conclusions  and  Recommendations 


The  conclusions  of  this  research  are  presented  as 
they  relate  to  the  questions  and  hypotheses  stated  in 
Chapter  I.  For  a  specific  mission  defined  by  the  input 
factors,  one  rendezvous  point  always  resulted  in  less  total 
fuel  consumption  than  any  other  point.  This  point  results 
in  a  significant  fuel  savings  compared  to  the  suboptimal 
points.  However,  the  optimal  point  was  not  always  located 
at  the  maximum  feasible  range  of  the  airlifter  from  its 
destination  as  proposed  by  hypothesis  one.  The  location 
of  the  optimal  rendezvous  point  can  only  be  determined  by 
considering  the  interaction  of  the  airlifter 's  total  dis¬ 
tance,  cargo  load  and  the  location  of  the  tanker  base. 
However,  definite  trends  in  the  results  emerged.  As  the 
airlifter  total  distance  to  destination  or  cargo  weight 
was  increased,  or  the  tanker  base  was  located  farther 
inland  of  the  airlifter  takeoff  base,  the  optimal  ren¬ 
dezvous  point  was  located  incrementally  closer  to  the  maxi¬ 
mum  feasible  range  boundary.  If  airlifter  dis-  ^  ce  to  des¬ 
tination  or  cargo  weight  were  reduced  or  '  i-ne  .  nker  base 
was -located  along  the  airlifter 's  route  of  flight,  the 
optimal  rendezvous  point  was  located  farther  from  the 
feasible  boundary  and  closer  to  the  tanker  base. 


The  second  hypothesis  is  also  rejected  in  favor  of 
a  proposal  to  determine  the  takeoff  fuel  loads  as  they 
relate  to  the  percentage  of  the  aircraft  total  fuel  capaci¬ 
ties.  Adjusting  the  transfer  fuel  until  the  sum  of  the  two 
percentages  of  fuel  capacities  is  the  minimal  sum  for  all 
feasible  fuel  loads  will  result  in  the  minimum  total  fuel 
consumed  by  both  aircraft.  This  minimum  sum  is  dependent 
on  the  aircraft  combination.  It  typically  occurred  at  the 
minimum  airlifter  fuel  loads  when  the  KC-lOA  was  used  and 
maximum  airlifter  fuel  loads  when  the  KC-135A  was  used. 

The  following  recommendations  are  submitted  as  a 
result  of  this  research: 

1.  The  planning  of  Strategic  Airlift  air-refueling 
missions  should  incorporate  the  results  of  this  research 
when  determining  the  rendezvous  point  and  takeoff  fuel 
loads.  Attempts  at  optimization  of  these  factors  will 
result  in  significant  fuel  savings  to  the  Air  Force. 

2.  To  aid  in  the  optimization  of  specific  airlift 
missions,  the  concepts  of  the  analytic  flight  planning  model 
used  in  this  research  should  be  expanded  to  include  three 
additional  dimensions: 

A.  A  real  time  weather  data  base, 

B.  Actual  route  segment  flight  planning  to 
replace  the  great  circle  route  planning  of  the  FLTPLN 
model;  and  finally. 


65 


C.  The  inclusion  of  the  possibility  of  multiple 
air-refuelings  of  a  single  airlifter  or  of  multiple  air- 
lifters  by  a  single  tanker. 

The  computer  flight  planning  programs  maintained 
by  AFGWC/DOY  provides  the  best  opportunity  to  consolidate 
all  of  these  dimensions  into  an  operational  air-refueling 
flight  planning  model. 
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Manual  Fuel  Consumption  Computation  of  a  Typ 
Airlift  Mission  Profile  Using  the  C-141B 


erformance  Manual 


This  appendix  contains  a  typical  fuel  planning  cal¬ 
culation  for  a  C-141B  mission.  Output  from  these  calcula¬ 
tions  will  be  compared  to  the  model  outputs.  (All  per¬ 
formance  figures  are  from  the  C-141B  performance  manual. 
Change  16,  24  Aug  1979.) 

The  mission  will  be  from  McGuire  AFB,  NJ,  to  Tehran, 
Iran.  The  tanker  will  be  a  KC-135  departing  and  returning 
from  and  to  McGuire.  The  following  steps  are  used  to  calcu¬ 
late  fuel  requirements: 

1.  INITIAL  CONDITIONS  (C-141B):  (obtained  from 

FLTPLN) 

RAMP  GROSS  WEIGHT  -  325,000  lbs. 

McGuire  to  rendezvous  entry  point  DISTANCE  - 
1350  NM. 

Rendezvous  track  =  250  NM. 

Rendezvous  exit  point  to  destination  DISTANCE  - 
1810  NM. 

TRANSFER  FUEL  -  53,000  lbs.  (JP-4) . 

WIND  FACTOR  -  55  KNOTS. 

2.  FUEL  CONSUMPTION  CALCULATIONS: 

The  following  section  calculates  mission  fuel  ‘ 
consumption  for  the  initial  conditions  in  (1) .  Fuel  con¬ 
sumption  calculations  consist  of  seven  areas: 

A.  START,  TAXI,  AND  TAKEOFF  FUEL  -  1900  lbs. 

B.  INITIAL  CLIMB  FUEL: 

(1)  INITIAL  PARAMETERS:  TAKEOFF  GROSS 
WEIGHT  (GW)  = 
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RAMP  GW  -  TAKEOFF  FUEL  = 

325,000  -  1900  =  323,100  lbs. 

CLIMB  ALTITUDE  =  33,000  ft, 

(2)  FROM  T.O.  C-141B  1-1  Chapter  4: 

CLIMB  FUEL  =  9050  lbs, 

CLIMB  DISTANCE  =  144  lbs. 

C.  CRUISE  FUEL: 

(1)  DISTANCE  from  level  off  to  rendezvous  point = 
1350  -  144  =  1206  NM. 

(2)  Fuel  flow  at  cruise  altitude  for  a  constant 
MACH  (0.74)  is  a  function  of  Gross  Weight  and  Temperature 
Deviation  from  a  Standard  Day  and  is  given  in  air  nautical 
miles  (ANM)  per  1000  lbs.  fuel. 

(3)  CRUISE  GW  =  TAKEOFF  GW  -  CLIMB  FUEL  = 

323,100  -  9,050  =  314,050  lbs. 

(4)  TRUE  AIRSPEED  (TAS)  AT  33,000  ft.  (.74) 

MACH  =  431  KNOTS 

GROUND  SPEED  (GS)  =  TAS  +  WIND  FACTOR  = 

431  +  55  =  486  KNOTS 

ANM  =  DISTANCE  *  TAS/GS  =  1026  NM  *  486  KNOTS/ 
431  KNOTS  =  1081  ANM 

(5)  FUEL  CONSUMPTION  RATE  =27  ANM/ 10 00  lbs. 

fuel 

FUEL  CONSUMPTION  =  1081  ANM/27  ANM/1000  lbs. 
fuel  =  40,000  lbs. 

D.  REFUELING  TRACK  FUEL  CONSUMPTION 

(1)  TRACK  DISTANCE  =  250  NM 

(2)  ENTRY  GW  =  CRUISE  GW  -  CRUISE  FUEL 

=  314,050  -  40,000  =  274,050  lbs. 

(3)  AIRSPEED  =  .75  MACH  9  25,000  ft.  altitude. 


(4)  FUEL  CONSUMPTION  RATE  =  2:0  lbs . /minute 

(5)  REFUELING  TIME  =  DISTANCE/AIRSPEED 

TAS  @  25,000  ft.  =  453  KNOTS 

GROUNDSPEED  =  TAS  +  WIND  FACTOR  = 

453  +  55  =  508  KNOTS 

TIME  =  DISTANCE/GROUNDSPEED 

TIME  =  250  NM  *  60/508  KNOTS  =29.5  min. 

FUEL  CONSUMPTION  =29.5  min.  *  260  lbs. /min 
7680  lbs. 

POST-RENDEZVOUS  CLIMB  TO  ALTITUDE 

(1)  INITIAL  CLIMB  GW  =  ENTRY  GW  -  TRACK  FUEL  * 

TRANSFER  FUEL 

274,050  +  53,000  -  7680  =  319,370  lbs. 

(2)  CLIMB  ALTITUDE  =  29,000  ft. 

(3)  From  the  C-141B  performance  manual; 

CLIMB  DISTANCE  =  27  NM 

CLIMB  FUEL  =  1400  lbs . 

CRUISE  TO  DESTINATION  FUEL 

(1)  CRUISE  GW  =  CLIMB  GW  -  CLIMB  FUEL 

=  319,370  -  1400  =  317,970  lbs. 

(2)  CRUISE  DISTANCE  =  RENDEZVOUS  EXIT  TO 

DESTINATION 

DISTANCE  -  CLIMB  DISTANCE  =  3810  -  27  = 

3783  NM. 

(3)  AT  29,000  ft.  TAS  =  438  KNOTS 
GS  =  438  +  55  =  493  KNOTS 

(4)  ANM  =  3783  *  438  /  493  =  3361  ANM 

(5)  FUEL  CONSUMPTION  RATE  =  29  NM/1000  lbs. 

fuel 

(6)  FUEL  CONSUMPTION  =  ANM/ RATE  =  3361  ANM/ 

29  ANM/1000  lbs.  =  115,900  lbs. 


G.  APPROACH  AND  LANDING  FUEL  IS  2500  lbs. 

The  total  planned  fuel  consumption  from  McGuire  AFB 
to  Tehran,  Iran  is: 

G 

^  (FUEL  CONSUMPTION)^  =  175,900  lbs. 
i=A 
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The  aircraft  ground  delays  in  the  SLAM  Simulation 


Model  were  based  on  the  data  shown  below.  The  data  was 
obtained  from  various  CONUS  Air  Force  bases  for  the  follow¬ 
ing  aircraft: 


AIRCRAFT:  KC-135 
PERIOD:  30  DAYS 
INSTALLATION:  MARCH  AFB,  CA 
TOTAL  TAKEOFFS:  93 
DELAYED  TAKEOFFS:  12 


DELAY  # 

LENGTH  OF  DELAY 
(MIN) 

DELAY  # 

LENGTH  OF  DELAY 
(MIN) 

1 

21 

7 

19 

2 

19 

8 

77 

3 

227 

9 

210 

4 

6 

10 

48 

5 

1 

34 

11 

74 

6 

236 

12 

1 

6 
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AIRCRAFT:  C-141 
PERIOD:  30  DAYS 
INSTALLATION:  NORTON  AFB,  CA 
TOTAL  TAKEOFFS:  156 
DELAYED  TAKEOFFS:  25 
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AIRCRAFT:  C-5 

PERIOD:  60  DAYS 

INSTALLATION:  DOVER,  DE 

TOTAL  TAKEOFFS:  143 

DELAYED  TAKEOFFS:  25 

DELAY  # 

LENGTH  OF  DELAY 
(MIN) 

DELAY  # 

LENGTH  OF  DELAY 
(MIN) 

1 

31 

14 

90 

2 

28 

15 

121 

3 

72 

16 

1152 

4 

30 

17 

243 

5 

35 

18 

44 

6 

37 

19 

36 

7 

29 

20 

162 

8 

187 

21 

1578 

9 

48 

22 

119 

10 

67 

23 

24 

11 

34 

24 

146 

12 

98 

25 

218 

13 

104 
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FLTPLN  Computer  Code 


C  THIS  IS  THE  HASTER  FLIGHT  PLftNHlNC  FROGRAH 
C 

PROCRAH  FLTPLN  (INPUT.OUTPUTI 
C 

OIHEHSION 

DIHENS  ION  LAT ( 7 1 . LONG ( 71 . PRL ( 7 ) f  HRN ( 7 1 . CARGO ( 4 ) 

D IHENS ION  LEG ( 5 1  .Fl'EL  (101. FEASEL  ( 5 . 1 3 1 . FUEL?! I N 1 !  0 ) 
DIMENSION  MACH(4).«XCMT(T).H](Fy£L(4).0FHTI4).RFFi4) 
DIMENSION  RHL(7).  RriACHIC).  STT0F(4)  .KK(!D) .LL113I 
REAL  LAT .LONG . HRN .MIGHT .HXFUEL . LEG . MACH 
INTEGER  AC.  TK 

DATA  PI. DEC. RAO/3. U1592t54. 57. 27577051.. 017453M52/ 
DATA  MACH  /.77.  .74.  .S2.  Ml 
DATA  RMACH/  .S2.  .74  / 

DATA  RFF  /.45.  .26.  .27.  .25/ 

DATA  HllGHT/712.5.323. 1.590. 0.270. 8/ 

DATA  Mr.FUEL/315.1.  151.45.  345.4.  163.0/ 

DATA  mil  354.8.  141.8.  240.47.  105.0/ 

DATA  STTOF/2.8.  1.0.  3.0.  2.0/ 

DATA  CARCO/100..85..76..70..55.,40./ 

C 

C  FUNCTION  STATEMENT  TO  CONVERT  RADIANS  TO  DEGREES 
ORADICOORI  =  {INT(COORH(COOR-lNT(COORtl/3.i}/D£C 
C  FUNCTION  STATEMENT  TO  CONVERT  DECREES  TO  RADIANS 
RDECIANGLEl  =  INr(ANCL£)+8.4*(AVGL£-INT (ANGLE)) 

C  GROUND  SPEED  CALCULATION  FROM  TRUE  AIRSPEED  AND  HIND 
CSISPD.ALT.TCt  -  TASP(SPD.ALT)iCOS(ASlN(HV*SIN(HD-TC) 
t  /TASP(SPD.ALT)))-UV»COS(HD-TC) 

C 

C  SCENARIO  INPUT  DATA  IN  DECREES.  MINUTES  TENTHS 
PRINT*."  EAST  SiCENARIO  t  I" 

PRINT*. “  TANKER  BASE  t  1" 

ELEVl  --0.0 
LAT(l)  =  38.216 
L0NC(11=  121.57 
LAT (21  =  35.27 
L0NC(2»=-139.126 
LAT(31  =  61.131 
L0NC(3t=  149.535 
LAT(4t  =  64.376 
L0NC(4I=  147.83 
ELEV2  ^0.8 

C  CONVERT  LAT/LiOG  DEGREES  TO  RADIANS 
DO  10  I  M.  4 
PRL(I)-'  DRAOlLATd)) 

MRf((I)=  DRADILONCID) 

If  CONTINUE 
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C  HIKD  DIRECT 10«  IK  RAOlAKi  AT  CD 
UD  : 

C  UIKD  VELOCITT  IN  N«  /  MINUTE  AT  KV 
UV  : 

C 

C  SELECT  KC-10A  TANKER  FIRST  THEN  KC-135A  TANKER 
DO  S99  NN  =  1.  2 
TK  =  3 

IF  (NN  .Ea.  2»  TK  =  4 
C 

C  SELECT  C-5A  AIRLIFTER  FIRST  THEN  C-141E  AIRLIFTER 
DO  938  KH  =  1.  6 
PRINT  m 

m  FORMAT  <ll(i/.120("»")) 

AC  =  1 

IF  (HH  .GE.  4)  AC  =  2 

PRINU,"  AC  =“iAC.“  TK  ="iTK."  CARGO  =".CARCO(HHI 
C  CALCULATE  RESERVE  FUELS 
RESV  =  3«./  AC 
REST  :  15,»(5.-TK) 

AZFN  =  OPNTIACI  +  CARCO(««) 

C  OETERHINE  FEASIBLE  EOUNDARt  RANGE 

CKFUEL  s  AHINU«)(FUEL(AC).(f1IGMT(ACl-A2FWI! 

CALL  RHOTHTA(PRL(l).KRN{U.PRL(2).HRH(2I.TDIST.CC.TCI 
ACU  =  HXCWKAO  -  (CKFUEL-RESVI/2. 

ALT2  =  CALT(AC«.ACiTCI 

CALL  FLT»EIAGH.MACH(AC) iALT2.4.CKFF,AC) 

LEG(2)  =  (CKFUEL-RESVHCS(«ACH(AC).ALT2.TC)/CKFF 
C  RPT  AT  level  off  IF  WITHIN  498  NN  OF  AIRLIFTER  EASE 
IF  (LEG(2I.LT.  (TDIST-4e3.)I  CO  TO  Id 
PRINT* t“  RPT  AT  LEVEL  OFF" 

AGU  -  «1CWT(AC) 

CALL  FLT  NE ( ACH . ELEV 1 . 25 . 1 2 1 CLND 1 ST  r  AC ) 

PRINT*."  CLKBIST  -".CLNBIST 
LEG (2)  =  TOiST-CLNDlST-350. 

2#  LEGdt  =  TDIST-LEC(2) 

C  COMPUTE  COORDINATES  OF  EOUNDART  RPT  ON  GREAT  CIRCLE 
CALL  LATLONG(PRL(l).HF.N(l).LEC(ll.GCifRL(7).MRNI7)) 

C  SET  SEPARATION  BETWEEN  RPT 
SN  :  5.  *  RAD 

IF  IPRLI7).CE.FRLI4))  SN  =  -SN 
EH  =  2.  *RAD 

IF  (SINIGO.CE.D)  EH  =  -EH 
C  SET  INITIAL  AIRLIFTER  FUEL  LOAD  AT  NAIINUN 
AFUEL  =  CKFUEL 
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C  DEVELOP  MATRICES  FOR  10  A1RLIF1ER  FLEL  LC.‘.D. 

DO  i0f  in  :  ]•  10 
PfiLIS)  =  PRL(7) 
nRNIS)  :  nRNI7) 

FUELniMdMl  =  777.777 
ACUT  :  AniNUAZFIl4AFUELiniCUT(ACM 
AFL  mT(ACUT*ACiTCI 

C  COMPUTE  CLIMB  VALUES  FOR  CURRENT  AIRLIFTED  FUEL  LOAD 
CALL  FLTMElACHTiELEVl.AFL. 1 iCLMT.AC) 

CALL  FLTME ( ACHT. ELEVl t AFL . 2 iCLMD . AC 1 
CALL  FLTNE (ACUT  >ELEV1 . AFLi 3 iCLMF . AC ) 

C 

C  GENERATE  5  I  13  MATRII  OF  RPT 
DO  80  Ldil3 
DO  70  KdiS 

C  COMPUTE  ROUTE  SEGMENT  DISTANCES 

1S0  C ALL  RHOTHT A ( PRL ( 1 ) . MRN ( 1 1 >  PRL i 5 ) I MRN 1 5 ) ) LEG ( 1 1 1 CC  t  RML 1 1 ) ) 
CALL  LATL0NG1PRL(5) >MRN(S) >250.1 (GC-PI ) iPRLIS) rMRNIS) ) 

CALL  RH0THTA(PRL(5))MRN(51iPRL(2l)MRNI2)tLEC(2)iCC.RML(21) 
CALL  RHOTHTA(PRL(G) fMRNtS) iPRL(3) iKRNI3) iLECI3) >CCiRML(3l ) 
CALL  RWTHTA(PRL(4)>MRNI4hPRL(S)>MRNIG)>LEC(4))CC)RKL(4)l 
CALL  R»)THTA(PRL{S)>MRN(5hPRL(4)iMRNi4)iLEC(5)>CCiRnLIS)) 
FEASBLlKiLl  •  777.777 
C 

C  AIRLIFTER  POST  RENDEIVOUS  TO  DESTINATION 
FAV  =120./  AC 

ACU  :  AniNl(AZFU>RESVi2*FAV>ni(GUT{AC)) 

DO  30  JU  Ml  3 
CLD  :  0.0 
CLF  :  0.0 

AFL  :  CALT(AGUfACiRHL(2)} 

IF  (AFL  .EQ.  25.)  CO  TO  25 
CALL  FLTRE(ACUi25.iAFLi2iCLDiAC) 

CALL  FLT«E(ACHi25.iAFLi3iCLFiAC) 

CLD  =  CLD»CS(MACH(ACIiAFLiRML(2))/TASP(MACH(AC)iAFLI 
25  ACM  :  AI1INt(AZFU^RESVtFAViMlCUT(AC)) 

CALL  FLTnE(ACUiHACH(AC)iAFL>4>AFFiAC) 
FUEL(21dLEC(21-CLD)»AFF/GSlMACH(AC)>AFLiRML(2)) 

ACU  =  At1INl(AZFU+FUEL(2)4CLF4RESViniCUTlAC)) 

FAV  =  FUEL(2)  /  2. 

30  CONTINIE 
AFL2  =  AFL 

FUEL(2)  =  FUEL(2t  +  CLF 

IF  (FUEL(2)  .CT.  (CKFUEL-RESV) )  CO  TO  40 


83 


C  AIRLIFTER  DEPARTURE  TO  RERDLZVOUS 
LEG(l)  :  LEC(t)  *  25t. 

AFL  :  CALT(ACyT>AC.RKLn)) 

ATIM  :  ILECID-CEtID)  /CSinACH(AC) )AFLiRHLIl)) 

ACU  :=  ACUT  -  CinF 

CALL  FLTnE{ACUiHACH(AC) iAFLi4iAFF)AC) 

ACU  =  ACM  -  (ATIH  »  AFF/  2.1 
CALL  FLT!1E(ACUiHAi:H(AC)iAFLi4iAFFiAC) 

FUELIU  =  ATIK  ♦  AFF  4  CLAF 
ATIH  =  ATIH  4  CLHT 
C 

C  AIRLIFTER  ABORT  TO  ALTERNATE 
ACU  :  ACUT  -  FUEL  in 
AFL  :  CALT(ACUiAC.RHL(3)) 

CALL  FLIKElACUiHACHIAC) lAFLiAiAFFiAC) 
ACU=ACU-(LEC(3»4AFF/CS(HACHlACItAFL.RflL(3n)/2. 

CALL  FLTKE ( ACU iHACH I AC) i AFL 1 4 • AFF . AC > 
FUEL(3):LEC(3)4AFF/CS(HACH(AC)iAFL>RKL{3)) 

C 

C  AIRLIFTER  ABORT  TO  DEPARTURE  BASE 
ACU  :  ACUT  -  FUEL  ID 
RHL(6)  =  RHLd)  -  PI 
AFL  :  CALTIACUtACiRHLan 
CALL  FLTHEIACUfHACHIAC) iAFL)4)AFFtACI 
ACUsACU-(LEC(n*AFF/CS(HACHIACI  .AFL.RHLItI )  )/2. 

CALL  FLTHE(ACU>HAi:H(AC)iAFLi4(AFF>AC) 
FUEL(6)=LEC(n*AFF/CS(HACH(AC).AFL.RHL(ni 
C 

C  AIRLIFTER  RENDEZVOUS  TRACK  FUEL 

FUEL ( 7) =258 .4RFF { AC ) /GS (RRACH ( AC I . 25 . . RHL ( 1) ) 

IF  ( IFUEL(n4AHlNl (FUELCJI .FUELI6I I4RESV4FUELI7) ) 
t  .CT.  AFUEL)  CO  TO  C8 
C 

C  TRFL  IS  THE  TRANSFER  FUEL  AT  RENDEZVOUS 

TRFL  =  FUEL(n4FUELI7)4FUEL12)4RESV  -  AFUEL 

IF  IAFUEL4TRFL-FUEL(n-FUEL17l  .CT.  CKFUELl  CO  TO  it 
C 

C  TANKER  POST  RENDEZVOUS  TO  DESTINATION 
FAV  =  (5.Z-TKI48.5 
TCUT  =  0PUT(TK)4REST 
TCU  =  AHINUTCUT  4  2.»FAV.HIGHT(TKI ) 

DO  4il  JJ  =  1,  3 
CLD  =  0.0 
CLF  =  0.0 

TFL  =  CALTiTCU.TK.RHLISn 
IF  ITFL  .E8.  25.)  CO  TO  35 
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CALL  FHI1E(TCHi25. i1FLi2.CLD.TK) 

CALL  FLFHElTCW.25.iTFLi3iCLF.TK) 

CLD  *  CLD*CS(HACKITK».TFLiRf!L(5ll/TASP(f!ACH(TKI.TFLI 
35  TCU  :  AniMKTCUT  -»  FAVi  HICUKTKM 
CALL  FL1!1E(TCHt«ACH(TKhTFLi4.TFF,Ti<l 
FUEL(51=  (LEC(5t-CLDI»TFF/CSlRACH(TK)iTFL.R(<LI5)' 
TC«=AHINl(TCMT+FUELt5HCLFt  RKCHKTKtl 
FAV  =  F1JEL(51  IZ. 

4«  CONTINUE 

FUELlSt  =  FUEL(5)  4CLF 
C 

C  TANKER  RENDEZVOUS  TRACK  FUEL 

FUEL (8) =250 . ♦RFF (TK) /CS (RKACH ( AC I . 25 . i RHL ( 1 > ) 

TON  =  At1INl(TCW+FUEL(8).  HJCUTITKII 
CALL  FLT«E(TC«.KACH(TK1 i25. tStWAlTF.TK) 

FUEL{81  =  FUEL(81  +  15.*  MAITF 
C 

C  TANKER  DEPARTURE  TO  RENDEZVOUS 
FAV  =  (4.3-TKH5.4 

TCNT  =  OPWT(TKHREST4FU£L(51+TRFL+FUELi8) 

TCN  =  AflINHTCMT+Z.iFAVi  RTGMTITK)) 

DO  50  JJ  =  1*  3 
TFL  =  CALT(TCHfTK,RRL{4n 
CALL  FLTNE(TCU.ELEV2,TFL.liCLT,TK) 

CALL  FLTHE(TCNtEL£V2iTFLi2.CLD.TK> 

CALL  FLTME(TCH.ELEV2.TFLi3tCLF.TK) 

TCN  --  AtlINKTCMT+FAV.  NICNTITKl) 

CALL  FLTt1ElTCNi«ACH{TK>.TFL.4.TFF.TI  ) 

TTIR  =  ILEC(4)-CLB)/CS(HACH(TK>.TFL.RRL14)I 
FUEL (41  =  TTIN  *  TFF 

TCN=  ANINUTCNT  +  FUEL(41  ♦  CLFi  MTCNTITKll 
FAV  =  FljEL(41  /  2. 

50  CONTINl£ 

FUEL{41  =  FUEL(4)  +  CLF 
TCNT  =  TCNT  +  FUEL (41 

IF  (TCNT  .CT.  «XCNT(TKI1  CO  TO  Sf 
C 

C  COMPUTE  TANKER  TAKE  OFF  FUEL  LOAD  AT  TFUEL 
TFUEL  =  FUEL(41+TRFL+FUEL(8)*FUEL(5)*REST 
IF  (TFUEL  .CT.  NXFUEL(TKI1  CO  TO  ti 
C 

C  STORE  TOTAL  FUEL  CONSURPTICN  IN  FEASBL  NATRU 
FEASBL ( K  f  L ) =  FUEL ( 1 ) +PLEL ( 2 1 *FUEL ( 7 1 
C  ♦FUEL(4HFUEL(51+FUELI8) 
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C  SELECT  KIKWJM  VALUE  AKD  SAVE  PA'^'ArEIERS 

FUELniHdKI  :  AnlNKFUELniNUnii  FEASCukiD) 

IF  (FUELHINdKI  .N£.  FEASSLIKiLd  CO  TO  ti 

KKIItl)  :  X 

LLdll)  -  L 

00  55  I  Ml  5 

Ud)  :  LECd) 

nd-»5>  =  RMLdI 

Ild-dl)  =  FOELII) 

Ild-dSl  *  FUELd+5» 

55  CONTIHUE 

lldZd  =  AFUEL  +  STTOFIACl 
n(22)  :  TFUEL  +  STTOF(TK) 

I)(<23)  =  TRFL 
IH2d  =  AFL2 
II(25J  =  e. 

imi>)  ■-  TTIW  ♦  15.  -  ATIM 
IF  (11(20)  50)  58i  58 

56  »(25)  =  'n(26) 

nm)  - 1. 

58  U(27)  -  PRL(6) 

11(28)  :  nRN(6) 

C 

C  INCREMENT  LATITUDE  COLUHN 
6f  PRL(5)  =  PRL(5)  ♦  SM 

70  COKTIHUE 

C  IHCREKENT  LOHGITUDE  AFTER  FIFTH  LATITUDE 
PRL(5)  5  PRL(7) 

HRN(S)  •  miS)  *  EU 
80  CONTINUE 

IF  (FUELHlNdH)  .EQ.  777.7771  CO  TO  100 
C 

C  PRINT  mm  VALUES  AND  OUTPUT  FOR  HINIHUN  RPT 
DO  WO  X  Ml  5 

PRINT  8001  (FEASBL(XiL)iLdil3) 

800  FORMAT  (Uil3(2IiF7.3)) 

900  CONTINUE 

LAT(6)  :  RDECdKZTXDEC) 

L0NC(6)  --  RD£C(ni28)}DEC) 

PRINT»i“  MMiFUELHIN(I«)tKK(IMl.LLdM)f 
(  “  RPT  EMTRT  ="tLAT(6liL0NC(6) 

PRIHT»i"  LEGS  ="i(ll(I)iIMi5) 

PRIMT»i"  TRUE  COURSES  =Mnd)iI=6il0) 

PRINT*!"  FUELS  -“ilIKDiIdlilS) 

ACUT  :  0PUT(AC)*CARGQ(MR)*n(2n 
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PRINT*."  AFUEL  =".n(Zl).“  ACMT  =".A&WT. 
t  “  AFLZ  :".n(24).“  ATOT  =‘‘.n(Zi) 

TCUT  :  OPUTITK)  *  HIZZI 
PRINT*."  TFUEL  -“.niZZh"  TCNT  =".TCMTi 
(  "  TRFL  =".n(231."  TTOT  =“.11(251 

C  DECREMENT  AIRLIFTER  TAKEOFF  FUEL  LOAD 
10i  AFUEL  =  AFUEL  -  2«. 

C  SELECT  OPTIMAL  RPT  FROM  NININUN  VALUES 
FUELIN  =  FUELHINtn 
DO  300  N  =  I.  10 

FUELIN  =  ANINliFUELIH.FUELHlN(NI) 

IF  (FUELIN  .NE.  FUELNIN(N))  CO  TO  300 
IN  =  H 
300  CONTINUE 

IF  (FUELIN  .NE.  777.7771  CO  TO  400 
PRINT*."  NO  FEASABLE  SOLUTION" 

CO  TO  998 
400  PRINT  700 

PRINT*."  FUELNIN(".IN.“.“.KK(IN).".".LL(INI.")  :".FUELNIN(IN1 
LAT(7t  =  R0EC(PRL171*BECI 
L0NC(7h  RDEC(NRN(7)*0EC} 

PRINT*."  POINT  ONE  =".LAT(7I.L0NC(71 
PRINT  700 

998  CONTINUE 

999  CONTINUE 
STOP 
END 
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c 

SUBR0UV1«E  RHOTHWIPl.PIW.PZ.PZN.RHOtTHETA.RXL) 

C  0)KPUT£S  ''.RtAT  CIRCLE  DISTRHCE  AND  COURSE 
60  =  1.5»79S3?7 
IF  (PI  .GT.  8D)  PI  s  6D 
IF  (P2  .CT.  601  PC  =  60 

0=ftC0S(SIN{PU«SI«(P2l4C0S(Pl)*C0StP2HiCD£iPCn-PlHl) 
RH0=D43437. 74677 

THETA=AC0S1  (SI«(P2)-S1N1P1  >4C0S(D)  1  /SIMiiJ)  /COSIPl)  1 
IF(SIN(P2M-PIH).CE.«}  THETA=  6.2£3185308-THETA 
RHL  -  ATAN({PIH-P2I1)/(AL0C(TAN;8.7853?81635+P2/211 
t-AL0C(TAN<«.7e539816354Pl/2) > 1 > 

IFIRRL.LT.I)  RKL=3.14l592t544PfL 
IF(SIM(P2H-P1H1 .CE.O)  R«L=3.14l5926544ftHL 
RETURN 
EMO 
G 

SUBROUTINE  LATL0NClPt.PlH.RH0.THETA.P2,P2«l 
C  COHPUTES  LAT/LONC  GIVEN  OIST  AND  COURSE  FROR  POINT 
R  -  RHQ  /  3437.74677 
P2  =  ASlH(SIN(Pn  »  COSIR)  ♦  COSlPl)  * 
t  SINIRI  4  COS(TH£TAI) 

D  *  AC0S((C0S(R)-SIN(PI)4SI»1P2)) 

I  /COS(PU/COSIP2)) 

IF  ISINITHETA)  .CE.  0.)  0  -  -0 
P2H  s  PIN  4  D 

IF  (P2N.CT.I.  .AND.  SINIP2H).LT.»-)  P2H=P2n-t. 283185308 
IF  {P2H.LT,».  .AND.  SIN(P2H).CT.B  1  P2N=6.283I853084P2H 
RETURN 
END 


88 


SLAM  (Simulation  Language  for  Alternative  Modeling) 
is  a  new  FORTRAN-based  simulation  language  which  allows 
simulation  models  to  be  created  in  three  world  views: 

1.  Network 

2.  Discrete 

3.  Continuous 

A  SLAM  model  consists  of  a  set  of  interconnected 
symbols  that  describe  the  operation  under  study.  SLAM  pro¬ 
vides  network  symbols  (see  Figures  D-1  to  D-10)  which  can  be 
used  to  build  models  and  which  can  be  translated  into  input 
statements  for  computer  processing.  SLAM  symbols  and  input 
statements  used  for  the  SLAM  Flight  Simulation  Model 
described  in  Chapter  IV  are  explained  here.  The  following 
symbols,  statement  formats,  and  definitions  are  taken  from 


Introduction 

to  Simulation  and 

SLAM  by  A.  Alan 

and  Claude  D, 

.  Pegden  (1979) ,  pp 

.  435-551. 

Network  Element 

Figure 

1. 

ASSIGN  NODE 

D-1 

2. 

CREATE  NODE 

D-2 

3. 

ENTER  NODE 

D-3 

4. 

EVENT  NODE 

D-4 

5. 

GOON  NODE 

D-5 

6. 

MATCH  NODE 

D-6 

7. 

QUEUE  NODE 

D-7 

8. 

TERMINATE  NODE 

00 

1 

Q 

9. 

REGULAR  ACTIVITY 

o 

1 

VO 

10. 

SERVICE  ACTIVITY 

D-IO 
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NODE  TYPE: 


ASSIGN 


FUNCTION : 

INPUT  FORMAT: 

SPECIFICATIONS : 


SYMBOL : 


Fig.  D- 


The  ASSIGN  node  is  used  to  assign  values 
to  SLAM  variables  (VAR)  at  each  arrival 
of  an  entity  to  the  node.  A  maximum  of  M 
emanating  activities  are  initiated. 


ASSIGN,  VAR=value,  VAR=value, . . . ,M; 


ENTRY  OPTIONS 

VAR  ATRIB  ( INDEX), XX(INDEX)  ,11,  where 

INDEX  is  a  positive  integer  or 
the  SLAM  variable  II. 

value  an  expression  containing  con¬ 

stants,  SLAM  variables,  or  SLAM 
random  variables. 

M  positive  integer. 


ASSIGN  Node  Description  Summary 
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NODE  TYPE: 


CREATE 


FUNCTION: 


INPUT  FORMAT: 

SPECIFICATIONS: 


SYMBOL : 

TF 


Fig.  D- 


The  CREATE  node  is  used  to  generate 
entities  within  the  network.  The  node  is 
released  initially  at  time  TF  and  there¬ 
after  according  to  the  specified  time 
between  creations,  TBC,  up  to  a  maximum 
of  MC  releases.  At  each  release  a  maximum 
of  M  emanating  activities  are  initiated. 

CREATE, TBC, TF,MA,MC,M; 


ENTRY 

TBC 


TF 

MA 

MC 


OPTIONS 

constant,  SLAM  variable,  or 
SLAM  random  variable. 

constant. 

positive  integer. 

positive  integer. 

positive  integer. 


CREATE  Node  Description  Summary 
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NODE  TYPE: 


ENTER 


FUNCTION: 


INPUT  FORMAT: 

SPECIFICATIONS: 

SYMBOL: 


Fig .  D 


The  ENTER  node  is  provided  to  permit  the 
user  to  enter  an  entity  into  the  network 
from  a  user -written  event  routine.  The 
node  is  released  at  each  entity  arrival 
and  at  each  user  call  to  subroutine 
ENTER (NUM) .  A  maximum  of  M  emanating 
activities  are  initiated  at  each  release. 

ENTER,  NUM,M; 

ENTRY  OPTIONS 

NUM  positive  integer. 

M  positive  integer. 


3.  ENTER  Node  Description  Summary 
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NODE  TYPE: 


EVENT 


FUNCTION : 


INPUT  FORMAT: 

SPECIFICATIONS 

SYMBOL : 


Fig . 


The  EVENT  node  causes  subroutine  EVENT 
to  be  called  with  event  code  JEVNT  at  each 
entity  arrival.  This  allows  the  user  to 
model  functions  for  which  a  standard  node 
is  not  provided.  A  maximum  of  M  emanating 
activities  are  initiated. 

EVENT, JEVNT,  M; 

:  ENTRY  OPTIONS 

JEVNT  positive  integer. 

M  positive  integer. 


-4.  EVENT  Node  Description  Summary 


94 


NODE  TYPE :  GOON 


FUNCTION:  The  GOON  node  provides  a  continuation 

node  where  every  entering  entity  passes 
directly  through  the  node. 

INPUT  FORMAT:  GOON,M; 

SYMBOL : 


Fig.  D-5.  GOON  Node  Description  Suiranary 
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NODE  TYPE: 

FUNCTION: 


INPUT  FORMAT: 

SPECIFICATIONS: 

SYMBOL : 


MATCH 

The  MATCH  node  is  used  to  delay  the  move¬ 
ment  of  entities  by  keeping  them  in  QUEUE 
nodes  (QLBLs)  until  entities  with  the 
same  value  of  attribute  NATR  are  resident 
in  every  QUEUE  node  preceding  the  MATCH 
node.  When  a  match  occurs,  each  entity 
is  routed  to  a  route  node  NLBL  that  corres 
ponds  to  QLBL. 

MATCH, NATR, QLBL/NLBL,QLBL/NLBL, .  .  .  ; 


ENTRY 

OPTIONS 

NATR 

positive. 

QLBL 

a  queue  node 

label . 

NLBL 

a  node  label 
node. 

for  any  type  of 

Fig.  D-6 .  MATCH  Node  Description  Summary 


NODE  TYPE: 


QUEUE 


FUNCTION : 


INPUT  FORMAT: 

SPECIFICATIONS: 


SYMBOL : 


Fig.  D 


The  QUEUE  node  is  used  to  delay  entities 
in  the  IFL  until  a  server  becomes  avail¬ 
able-  The  QUEUE  node  initially  contains 
IQ  entities  and  has  a  capacity  of  QC 
entities . 

QUEUE (IFL) , IQ, QC, BLOCK; 


ENTRY 

OPTIONS 

IFL 

integer  between  1  and 

MFIL. 

IQ 

non-negative  integer. 

QC 

integer  greater  than 
to  IQ. 

or  equal 

SLBLS 

the  labels  of  MATCH  nodes 
separated  by  commas. 

7.  QUEUE  Node  Description  Summary 
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NODE  TYPE :  TERMINATE 

FUNCTION :  The  TERMINATE  node  is  used  to  destroy 

entities  and/or  terminate  the  simulation. 

All  incoming  entities  to  a  TERMINATE 
node  are  destroyed.  The  arrival  of  the 
TCth  entity  causes  a  simulation  run  to 
be  terminated . 

INPUT  FORMAT;  TERMINATE ,TC ; 


SPECIFICATIONS :  ENTRY  OPTIONS 

TC  positive  integer. 


SYMBOL : 


Fig.  D-8.  TERMINATE  Node  Description  Summary 
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ACTIVITY  TYPE: 

FUNCTION : 


INPUT  FORMAT: 

SPECIFICATIONS : 


SYMBOL: 

Fig.  D-9. 


REGULAR 


A  REGULAR  activity  is  any  activity 
enamating  from  a  node  other  than  a  QUEUE 
node.  The  REGULAR  activity  is  used  to 
delay  entities  by  a  specified  duration, 
perform  conditional/probabilistic  testing, 
and  to  route  entities  to  non-sequential 
nodes . 


ACTIVITY/A, duration,  PROB  or  COND,NLBL; 


ENTRY  OPTIONS 

A  positive  integer. 


duration  constant,  SLAM  variable,  SLAM 
random  variable. 

PROB  or  probability:  constant  between  0 
and  1 . 


COND  condition:  value  .OPERATOR. 

value  where  value  is  a  constant, 
SLAM  variable,  or  SLAM  random 
variable  and  OPERATOR  is  LT,  LE, 
EQ,  GE,  GT,  OR  NE. 

NLBL  the  label  of  a  labeled  node 

which  is  at  the  end  of  the 
activity. 


■> 


REGULAR  Activity  Description  Summary 
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ACTIVITY  TYPE: 


SERVICE 


FUNCTION: 

INPUT  FORMAT: 

SPECIFICATIONS : 


SYMBOL : 


Fig.  D-10. 


The  SERVICE  activity  is  any  activity 
emanating  from  a  QUEUE  node.  The  service 
activity  is  used  in  conjunction  with  the 
QUEUE  node. 

ACTIVITY (N) /A, duration, PROB,NLBL; 


ENTRY  OPTIONS 


N 


positive  integer. 


A 


positive  integer  between  1 
and  50. 


duration  constant,  SLAM  variable, 

SLAM  random  variable. 


probability  constant  between  0  and  1. 

NLBL  label  of  a  labeled  node. 


SERVICE  Activity  Description  Summary 
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C  SLflH  COMPUTER  CODE 

GEH*HflRCOTTE.THESISill/ll/80i9iTEStNOtTESiHO,!;Oi 

LI«lTSi2i5f5i 

net; 

CREATE  1 1 1 1 1 1 1 ; 

GOONili 

EVENTii; 

term; 

ENTER* i; 

ASl  ASSIGNflUi 
ACT*>>AS3; 

ACTillt(lt)itAS4; 

ENTER(2: 

ASSIGN* I 1=2; 

ACT***AS4* 

ACT*I(I((12;**AS3* 

AS3  ASSICN*ATRIB(l):X]((12)iATRlB(S)--l.; 

ACT*15*ATR!6(2).E&.1.0R.ATRlB(2).E3.3*C0Ui 
ACT*15*ATRIB(2).E0.2.0R.ATRIB(2).EQ.4*CG33i 
GOll  GOON*l 

ACT*1S*.18*ASI1* 

act**.82«asn; 

C033  COON*l 

ACT*15t.ltASn* 

ACT**.8*ASN 

ASH  ASSIGN*n(15)-t5.* 

ACT  M*  ASS 

ASH  ASSICN*XX(lS)-0.: 

ACT***AS5 

ASS  ASSIGN  rATRl  B  (4)  rUSERF  (I )  *  Cl  1 17 )  1 .  *  1 1 

ACT*3.*II.E0.2..AHD.XI(U8I.EQ.1)AS7; 
ACT**Il.E8.2..AKD.XT(18I.NE.l.iDi; 
ACT*3.t*AS7 
01  GOON* I 

ACT*nn2(+3**AS7; 

AS7  ASSlGN*ATRlB(4)=ATRIB(4t+USERF(3) I 
ACT*USERF(5)**ASS; 

AS4  ASSIGN.ATRIB(l)=n(ll)*ATRIBI5)=2.; 

ACT*l5*ATRI8(2).E8.3.0R.ATRlB(Zt.EQ.4*G022; 
ACT*15*ATRIB(2) .EQ.l.OR.ATRIBIZ) .EQ.2iG044) 
C022  GOON* I 

ACT*15*.18*AS8; 

ACT* *.82*  C044* 

C044  C00N*1 

ACT*15*.3*ASQ: 

ACT  *  * .7 lASR 

ASa  ASSIGN*Xm6)M5.; 

ACTif*ASZ 
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ASR  ASSIGNiX)((16)=e.t 
ACTmiASZ 

ASZ  ASSICNiATRIB(4)^USERF(Z).XX(18)M.ii; 
ACT)3fII.EQ.l..AND.XX(17l.£0.1.tAS/,; 
ACTiiIl.EQ.t..AND.tX(17!.NE.l.iDZ; 

ACTi3iiAS6; 

02  GOONfl 

ACTtXX(ll)+3.tASfe! 

AS&  ASSICNfATRie(4)=ATRlB(4)^USERFI4); 

ACTiUSERF(6)5 

AS8  ASSIGNtATRIBi4hATRie(4)tUSERFI8)iII(2e):USERF(10); 
ACTtUSERFdZi; 

ASA  ASSIGNiATRIB(4)rATRlB(4KUSERF(14l)It(22)=TN0U; 

ACT) ( (82 

AS9  ASSIGNiATRIB{4>=ATRlB(4)tUSERF(7),XXll9)=USERF{9)) 
ACTfUSERFUli; 

ASB  ASSIGN)ATRlB(4)=ATRlB(4»+USERFd3lrXX(21)=TN0Wi 
Q1  QUEUEdItO; 

82  8UEUE(2lf8)‘ 

NTl  RATCHi2f81/G05)82/C0i8; 

GOS  GOONfli 

ASC  ASSIGKiATRIB(4|:ATRlB(4)^l.fti; 
ACTi)USERFd7I.LT.0.)Tl! 

ACTfffASX 

ASX  ASSIGH)ATRIB(3):USERF(19): 

ACTfATRIB(3)«)ASEi 
T1  TERHid 

ASE  ASSIGN)ATRIB(4):ATRIB(4)«USERF(2U«XX(27)‘ATRIBI3)*5i 
ATRIB(3)=USERF123)? 

ACTfATRI6(3)) 

ASC  ASSIGN) ATRI6 (4) =ATRIB (4) ^USERF (25)1 
T3  TERR! 
coif  G00N)1) 

ACT! 

ASD  ASSICN)ATRIB(4)=ATRIB(4)+USERFd6))i; 
ACT))USERFd8).LT.f.)T2; 

ACT)))AST) 

AST  ASSICN)ATniB(3l=USEfiF(l9i; 

ArT)ATRIB(3l))ASF) 

T2  TERR) 

ASF  ASS:CN.ArH;B(41=H:Rii;!dHlSERF(2B))XX(27)--ATRIB(3)t5) 
ATRIE(i)-USERF(22): 

ACT)ATRI[.-:)t 

ASH  ASSICN)AT'iIB:4>=ATR1B(41+USERF(241! 

T4  TERR! 

ENDNET! 

IN1T)9!9?3) 
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i 

iMntf>990;  I 

NOt<TRiTRACE>d»9»0; 

SEDSi9375Z95(l) 

SinULATE 

SEE0S*-9375295U) 

NO«TRfCLEAR.0 
nOttTRrTRACEi9>900 
SinULATE 
SEE0St9375295(2) 
nONTRiCLEAR.A 
nOt<TRtTRACEt0r900 
SIMULATE 

SEEDSf-9375295(2) 

ROHTR>CLEARfS 
HOMTR>TRACEi0i9l» 

SinULATE 
SEEOSf 9375295 (3) 

HONTRiCLEARifl 
HONTRiTRACw  yf9M 
SINULATE 
SEE0Si9375295I4) 

HONTRtGLEARtfl 
nONTRfTRACE>8>90f 
SIMULATE 

SEEDS 1-9375295(3) 

HOMTRtaEhRig 
MOMTRiTRACEi0r90l 
SIMULATE 

SEEOSi -9375295 (4) 

HOMTR>CLEARi0 
MONTR.TRACEf0r90f 
SIMULATE 
SEE0Sr9375295(5) 

MOMTRfCLEAR.0 
HOMTRtTRA€Er0>900 
SINULATE 
FIN 


J 
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Appendix  F 

User  Function  Summaries 
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used  to 


¥■ 

«> 

1.  User  function  one  J[US^F1)  (Figure  F-l)\  is 
calculate  the  fuel  required  ^uring  all  airlifter 

ground  operations.  Normal  ground  operations  consist  of 
start,  taxi,  takeoff  (STTO)  and  delays.  STTO  fuel  is  a 
constant  value  for  a  particular  aircraft  shown  in  Figure 
F-26.  USERFl  returns  fuel  consumed  in  thousands  of  pounds 
of  fuel  using  the  equation: 

USERFl  =  STF  +  XX(15)*  FF  where 

XX(15)  =  delay  time  (min) 

STF  =  STTO  fuel  (1000  lbs) 

FF  =  ground  idle  fuel  consumption  (1000  Ibs/min) 
Ground  idle  fuel  consumption  for  all  aircraft  is  shown  in 
Figure  F-26  and  obtained  from  the  aircraft's  performance 
manuals . 

2.  User  function  two  (Figure  F-2)  is  used  to  cal¬ 
culate  the  fuel  required  during  all  ground  operations  for 
the  tankers.  Calculations  are  similar  to  USERFl  and  will 
not  be  repeated.  XX (16)  is  the  tanker  ground  delay. 

3.  User  function  three  (Figure  F-3)  calculates 
the  takeoff  fuel  from  brake  release  until  landing  gear 
retraction  for  the  airlifter.  Values  used  are: 

USERF  =  1000  lbs  fuel  for  the  C-5A 

500  lbs  fuel  for  the  C-141B 

^All  of  the  figures  referred  to  in  this  appendix 
appear  at  the  end  of  the  appendix. 
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4.  User  function  four  calculater-  te  takeoff  fuel 
for  the  tanker  (Figure  F-4) .  Values  used  are: 

USERF  =  1000  lbs  for  the  KC-10 
500  lbs  for  the  KC-135 

5.  User  function  five  (Figure  F-5)  is  used  to 
return  the  initial  climb  time  for  the  airlifter  computed 
by  FLYME  using: 

XX (10)  =  Airlifter  current  gross  weight 

XX (32)  =  Cruise  altitude  calculated  by  sub¬ 
routine  CALT  based  on  XX (1)  and 
IBIRD. 

XX(31)  =  Airlifter  departure  base  elevation. 

6.  User  function  six  (Figure  F-6)  calculates 
tanker  time  (sec)  in  the  initial  climbout.  Calculations 
are  similar  to  user  function  five  except  for  the  following 
input  variables  to  FLYME: 

XX  (9)  replaces  XX (10)  as  the  current  tanker 
gross  weight. 

XX (42)  replaces  XX (32)  as  the  tanker  cruise 
altitude . 

XX(41)  replaces  XX(31)  as  the  tanker  departure 
base  elevation. 

7.  User  function  seven  (Figure  F-7)  computes  the 
airlifter  initial  climb  fuel  using  FLYME  and  the  same  input 
variables  as  USERF5.  The  value  returned  is  in  thousands 

of  popnds  of  fuel.  Also  XX(35)  is  set  equal  to  the  climb 
fuel . 
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8.  User  function  eight  (Figure  F-3)  computes  the 
initial  tanker  climb  fuel  using  FLYME  with  the  USERF6 
input  parameters.  Also  XX (40)  is  set  equal  to  the  climb 
fuel . 

9.  User  function  nine  (Figure  F-9)  computes  the 
airlifter  climb  distance  (NM)  using  FLYME  and  USERF5  input 
parameters.  USERF9  also  sets  XX(19)  equal  to  the  climb 
distance . 

10.  User  function  ten  (Figure  F-10)  computes  the 
tanker  initial  climb  distance  (NM)  using  FLYME  with  input 
parameters  of  USERF6.  USERFIO  also  sets  XX (20)  equal  to 
climb  distance. 

11.  User  function  11  (Figure  F-11)  calculates  the 
airlifter  cruise  time  (min)  from  leveloff  to  rendezvous 
entry  point.  Variables  used  are; 

a.  XX (1)  =  distance  from  airlift  departure  to 

rendezvous  entry  point. 

b.  AS  =  airlifter  true  airspeed  at  altitude 

XX(32)  (NM/MIN) 

USERFll  =  XX (33)  =  CRUISE  TIME  = 

XX(1)  -  XX(19)  ^ 

AS  +  RNORM(l,.16) 

RNORM  =  Cruise  altitude  winds 

(y  =  1  NM/MIN;  a  =  . 16  NM/MIN) 

The  parameter  for  these  cruise  winds  was  obtained 


from  Global  Weather  Central  AFGWC/DOY,  Offutt  AFB,  NE . 
The  winds  are  modeled  as  constant  between  25,000  ft  to 
41,000  ft.  The  wind  variance  of  .16  NM/MIN  was  taken 


as  an  average  of  wind  velocity  changes  o'  ■  the  routes 
selected . 

12.  User  function  12  (Figure  F-12)  determines  the 
tanker  cruise  time  from  leveloff  to  the  rendezvous  entry 
point  in  the  same  manner  as  USERFll.  XX (43)  is  set  equal 
to  the  cruise  time  and  the  user  function  returns  time  in 
minutes.  An  undefined  variable  used  in  USERF12  is  XX(4) 
which  is  the  distance  between  the  tanker  departure  point 
and  the  rendezvous  entry  point. 

13.  User  function  13  (Figure  F-13)  computes  the 
airlifter  cruise  fuel  from  leveloff  to  rendezvous  entry 
point.  The  steps  used  to  compute  the  fuel  are; 

a.  Compute  the  initial  cruise  gross  weight: 

XX  (10)  =  Ramp  gross  weight  (XX (75))  -  ATRIB(4) 
ATRIB(4)  =  total  fuel  consumed  up  to  leveloff. 

b.  Compute  the  fuel  consumption  rate  (V)  at  cruise 
altitude  XX(32).  The  rate  is  returned  from 
FLYME  in  thousands  of  pounds/minute. 

c.  Cruise  fuel  =  V  *  XX(33)  fuel/min  *  min  = 

1000  lbs  fuel. 

Finally,  USERF13  decrements  the  aircraft  gross  weight  by 
the  cruise  fuel. 

14.  User  function  14  computes  (Figure  F-14)  the 
tanker  cruise  fuel  from  level  off  to  the  rendezvous  entry 
point  in  an  analogous  manner  to  user  function  13.  The 
user  function  returns  cruise  fuel  in  thousands  of  pounds. 

15.  User  function  15  (Figure  F-15)  computes  the 
airlifter  loiter  fuel.  The  computation  is  computed  as 
follows : 
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a.  Determine  loiter  time  (ATIME) .  This  is  deter¬ 
mined,  after  both  aircraft  have  reached  the 
match  node,  MATl,  by  computing  TNOW  -  XX (21), 
where  XX (21)  is  the  arrival  time  of  the  air- 
lifter  and  TNOW  is  the  arrival  at  MATl  of  the 
last  aircraft.  If  the  expression  is  negative 
or  zero,  the  user  function  is  returned  as  zero 
and  there  is  no  loiter  time.  If  the  expression 
is  positive,  then  airlifter  loiter  occurred  and 
the  loiter  fuel  is  computed. 

b.  Compute  the  loiter  fuel.  The  loiter  fuel  con¬ 
sumption  rate  is  determined  from  FLYME  and 
returned  as  V.  The  loiter  fuel  is  then  com¬ 
puted  as  "ATIME  *  V"  and  returned  in  thousands 
of  pounds . 

16.  User  function  16  (Figure  F-16)  computes  the 
tanl<er  loiter  fuel  in  a  manner  analogous  to  user  function 
15.  XX (22)  is  the  tanlcer  arrival  time  at  Q2. 

17.  User  function  17  (Figure  F-17)  calculates  the 
rendezvous  abort  fuel  for  the  airlifter  (1000  lbs  fuel). 

If  this  user  function  is  less  than  zero  the  airlifter  does 

« 

not  have  the  required  fuel  to  return  to  the  abort  base  and 
a  refueling  is  not  possible.  Therefore,  the  mission 
becomes  infeasible  and  the  model  is  terminated.  The  user 
function  is  computed  by: 

XX(8)  =  ATRIB(4)  -  XX(23)  where 

XX(8)  =  Airlifter  ramp  fuel 

XX (23)  =  Flight  plan  fuel  to  fly  from  rendezvous 

entry  point  to  the  abort  base,  calculated 
by  FLTPLN. 

18.  User  function  18  (Figure  F-18)  serves  the  same 
purpose  for  the  tanker  as  USERF17  serves  for  the  airlifter 
and  returns  the  remaining  tanker  fuel  after  the  abort  fuel 
has  been  subtracted.  XX (7)  is  the  tanker  ramp  fuel 
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and  XX (24)  is  the  flight  plan  fuel  to  fly  from  rendezvous 
entry  point  to  the  abort  base,  calculated  by  FLTPLN. 

19.  User  function  19  (Figure  F-19)  computes  the 
refueling  track  time.  The  track  is  250  NM  at  25,000  feet 
altitude.  The  time  is  computed  as  follows: 

a.  Track  true  airspeed  (TASTRK)  is  computed  from 
SUBROUTINE  TASP. 

b.  Compute  track  time; 

2  50  NM 

^TASTRK  +  WIND^  ^  MINUTES 

TASTRK  =  True  airspeed  during  refueling. 

20.  User  function  20  (Figure  F-20)  determines  the 
tanker  track  .fuel  consumption.  The  air-refueling  track 
fuel  consumption  rates  for  each  aircraft  were  obtained 
from  each  aircraft's  performance  manual.  The  rates  are 
shown  in  Figure  F-26.  Track  fuel  consumption  is  computed 
as: 

USERF20  =  (TRACK  TIME  *  FUEL  CONSUMPTION  RATE) 

=  Thousands  of  pounds  of  fuel. 

21.  User  function  21  (Figure  F-21)  computes  the 
airlift  refueling  track  fuel  consumption  similar  to  USERF20. 

22.  User  function  22  (Figure  F-22)  determines  the 
tanker  cruise,  approach,  and  landing  time  from  rendezvous 
exit  to  destination  as  follows: 

a.  Compute  an  average  gross  weight  (GW)  from 
rendezvous  exit  to  landing. 

AVG  GROSS  WT  =  (RENDEZVOUS  EXIT  GW  -  TRANSFER 
FUEL  +  TANKER  ZERO  FUEL  +  RESERVES)  r  2 


r 


The  gross  weight  at  the  end  c  the  flight  is 
computed  as  XX(95)  +  XX(46)  where  XX(95)  = 
zero  fuel  wt  and  XX (4 6)  =  tanker  fuel  reserves. 
Using  an  average  gross  weight  v;ill  permit  the 
tanker  to  initially  climb  higher  than  allow¬ 
able  at  his  track  exit  gross  weight.  In  this 
way  the  model  approximates  a  step  climb  profile. 

b.  Cruise  altitude  to  destination  is  determined 
from  SUBROUTINE  CALT  using  step  A  gross  weight. 

c.  Climb  time  (CT)  from  25,000  ft  (refueling 
altitude)  to  cruise  altitude  is  determined  from 
SUBROUTINE  FLYME . 

d.  The  climb  distance'  (CD)  is  determined  from 
FLYME. 


e.  The  climb  fuel  is  determined  from  FLYME  and 
set  equal  to  XX(60). 

f.  The  gross  weight  is  decremented  by  the  climb 
fuel . 


g.  The  cruise  distance  (CRDIST)  from  the  level- 
off  point  to  tanker  destination  is  computed  as 
the  distance  from  the  rendezvous  exit  point  to 
destination  minus  the  climb  distance  (CD)  . 


h.  The  time,  XX(62),  is  computed  as 

CRDIST 


XX (62)  = 


TASK  -  WIND  FACTOR 


+  15  minutes 


The  15  minutes  is  a  constant  time  all  aircraft 
use  for  approach  and  landing. 


i.  USERF22  =  XX(62)  +  CT  =  CRUISE  TIME  +  CLIMB 
TIME 


23.  User  function  23  (Figure  F-23)  computes  the 


airlift  time  (min)  from  rendezvous  exit  point  to  landing 


at  the  airlift  destination.  Computations  are  analogous 


to  user  function  22.  Variables  used  are  defined  in 
Figure  K-3,  Appendix  K. 
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24.  User  function  24  (Figure  F-24)  computes  the 
tanker  fuel  consumed  from  the  rendezvous  exit  point  to 
destination.  Computations  are  as  follows: 

a.  Determine  the  approach  and  landing  fuel  for 
the  tanker.  This  value  changes  by  aircraft; 
values  used  were  taken  from  AFR  60-16  MAC 
Supplement  One,  except  for  the  KC-10  which  was 
estimated  by  HQ  SAC/D08. 

b.  Determine  cruise  fuel  consumption  rate  (R) 
from  FLYME  using  the  average  gross  weight 
computed  by  USERF22. 

c.  Compute  cruise,  approach,  and  landing  fuel: 

USERF24  =  CRUISE  TIME  *  FUEL  CONSUMPTION  RATE 
+  APPROACH  AND  LANDING  FUEL  +  CLIMB  FUEL 

=  XX(62)  *  R  +  AAL  +  XX(60) 

25.  User  function  25  (Figure  F-25)  computes  the 
airlifter  fuel  from  rendezvous  exit  to  destination.  Compu¬ 
tations  are  analogous  to  user  function  24.  Variables  used 
are  defined  in  Figure  K-3  in  Appendix  K. 


FUNCTION  USERF(IFN) 

COMMON/ SCOMl/ATRIB (100) ,DD(100) ,DDL(100) ,DTNOW, 
III, MFA, MSTOP , NCLNR, NCRDR, NPRINT , NNRUN , 
INNSET,NTAPE,SS (100) ,SSL(100) , TNEXT , TNOW , XX ( 10 0 ) 
COMMON  QSET 

GO  TO  (1,2,3,4,5,6,7,8,9,10,11,12,13,14 
1,15,16,17,18,19,20,21,22,23,24,25) ,IFN 
C  GND  FUEL  FLOW  &  TAXI  FUEL-AIRLIFT 
1  IF(ATRIB(2)  .EQ.l.  .OR.ATRIB(2)  .EQ  .  3  , )  FF=: ,  12 

IF(ATRIB(2) .EQ.2. .OR.ATRIB(2) . EQ . 4 . ) FF= . 08 
IF{ATRIB(2) .EQ.l . .OR.ATRIB (2) .EQ. 3. ) STF=1. 8 
IF(ATRIB(2) .EQ.2. .OR.ATRIB(2) . EQ . 4 . ) STF=1 . 4 
USERF=STF+XX (15) *FF 
RETURN 

Fig.  F-1.  USERF  (1) — Ground  Fuel  Flow  & 

Taxi  Fuel-Airlifter 


C  GND  FUEL  FLOW  &  TAXI  FUEL-TANKER 
2  IF(ATRIB(2) .EQ.l. .OR.ATRIB (2) .EQ.2.)FF=.l 

IF(ATRIB(2) .EQ.3. .OR.ATRIB (2) .EQ.4)  FF=.12 
IF(ATRIB(2) .EQ.l. .OR.ATRIB(2) . EQ . 2 . ) STF=1 . 5 
IF(ATRIB(2) .EQ.3. .OR.ATRIB (2) . EQ . 4 . ) STF=2 . 0 
USERF=STF+XX(16) *FF 
RETURN 

Fig.  F-2.  USERF  (2)— Ground  Fuel  Flow  & 
Taxi  Fuel-Tanker 


C  TAKEOFF  FUEL  AIRLIFT 

3  IF(ATRIB(2) .EQ.l. .OR.ATRIB(2) . EQ . 3 . ) USERF=1 . 0 

IF(ATRIB(2) .EQ.2. .OR.ATRIB (2) . EQ . 4 . ) USERF= . 5 
RETURN 

Fig.  F-3.  USERF  (3) — Takeoff  Fuel  Airlift 


C  TAKEOFF  FUEL  TANKER 

4  IF(ATRIB(2) .EQ.l, .OR.ATRIB (2) . FQ . 2 . ) USERF= . 5 

IF(ATRIB(2) .EQ.3. .OR.ATRIB (2) . EQ . 4 . ) USERF=1 . 0 
RETURN 


Fig.  F-4.  USERF  (4)— Takeoff  Fuel  Tanker 


O  Ln  on 


C  CLIMB  TIME  AIRLIFT 
5  I0PT=1 

IBIRD=1 

IF(ATRIB(2) .EQ.2. .0R.ATRIB(2) . EQ . 4 . ) IBIRD= 2 
51  XX(10)=XX(75)-ATRIB(4) 

XX(32)=CALT(XX(10) , IBIRD, XX ( 8 1) ) 

CALL  FLYME(XX(10) ,XX(31) ,XX{32) , lOPT , V, IBIRD) 

USERF=V 

RETURN 

Fig.  F-5.  USERF  (5) — Climb  Time  Airlift 


CLIMB  TIME  TANKER 
lOPT  =  1 

IF(ATRIB(2) .EQ,1.0R.ATRIB(2) .EQ.2.)G0  TO  55 
IBIRD=3 
GO  TO  56 
55  IBIRD  =  4 
6  XX(9)=XX(76)-ATRIB(4) 

XX(42)=CALT(XX(9) , IBIRD, XX (84) ) 

CALL  FLYME  (XX ( 9 ) , XX ( 4 1 ) ,XX ( 4 2) , lOPT , V, IBIRD) 
USERF=V 
RETURN 

Fig.  F-6.  USERF  (6) — Climb  Time  Tanker 


CLIMB  FUEL  AIRLIFTER 
IOPT=3 
IBIRD=1 

IF(ATRIB(2) .EQ.2. .OR.ATRIB(2) .EQ.4.) IBIRD=2 
58  CA^^  FLYME(XX(10) ,XX(31) ,XX(32) ,IOPT,V, IBIRD) 
USERF=V 
XX(35)=V 
RETURN 

Fig.  F-7.  USERF  (7)— Climb  Fuel  Airlifter 


C  CLIMB  FUEL  TANKER 
8  IOPT=3 

IBIRD=3 

IF(ATRIB(2) .EQ.1.0R.ATRIB(2) .EQ.2) IBIRD=4 
60  CALL  FLYME(XX(9) ,XX(41) ,XX(42) ,IOPT,V, IBIRD) 
USERF=V 
XX(40)=V 
RETURN 

Fig.  F-8.  USERF  (8) — Climb  Fuel  Tanker 


C  CLIMB  DISTANCE  AIRLIFTER 
9  I0PT=2 

IBIRD=1 

IF(ATRIB(2) .EQ.2. .0R.ATRIB(2) ,EQ.4) IBIRD=2 
62  CALL  FLYME(XX(10) ,XX(31) ,XX(32) ,IOPT,V,IBIRD) 
XX(19)=V 
USERF=V 
RETURN 

Fig.  F-9.  USERF  (9) — Climb  Distance  Airlifter 


C  CLIMB  DISTANCE  TANKER 
10  IOPT=2 

IBIRD=3 

IF(ATRIB(2) .EQ.l. .OR.ATRIB(2) .EQ.2. ) IBIRD=4 
64  CALL  FLYME(XX(9) ,XX(41) ,XX(42) ,IOPT,V,IBIRD) 
XX(20)=V 
USERF=V 
RETURN 

Fig.  F-10 .  USERF  (10) — Climb  Distance  Tanker 


C  CRUISE  TIME  AIRLIFTER 
11  AMACH=.77 

IF{ATRIB(2) .EQ.2. .OR.ATRIB(2) . EQ . 4 ) AMACH= . 74 
AS=TASP(AMACH,XX(32)  ) 

USERF=( (XX (1) -XX (19) ) ) / (AS+RNORM ( 1 . 0 , .16,1) ) 

XX (33) =USERF 
RETURN 

Fig.  F-11.  USERF  (11) — Cruise  Time  Airlifter 


C  CRUISE  TIME  TANKER 
12  AMACH=.82 

AS=TASP ( AMACH , XX ( 4  2 ) ) 

USERF=( (XX (4) -XX (20) ) ) / (AS+RNORM ( 1 . 0 , .16,1) ) 

XX(43)=USERF 

RETURN 

Fig.  F-12.  USERF  (12) — Cruise  Time  Tanker 
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C  AIRLIFT  CRUISE  FUEL 
13  I0PT=4 

IF(ATR1B(2) .EQ.2.0R.ATRIB(2) ,EQ.4)  GO  TO  69 

AMACH=.77 

IBIRD=1 

GO  TO  70 

69  IBIRD=2 
AMACH=. 74 

70  XX(10)=XX(75)-ATRIB(4) 

CALL  FLYME(XX(10) ,AMACH,XX (32) , lOPT , V, IBIRD) 
USERF=V*XX(33) 

XX  (10)  =XX  (10)  -V 
RETURN 

Fig.  F-13.  USERF  (13)— Airlift  Cruise  Fuel 


C  TANKER  CRUISE  FUEL 
14  IOPT=4 

IF(ATRIB(2) .EQ.1.0R.ATRIB(2) .EQ.2)GO  TO  71 
IBIRD=3 
AMACH= . 8  2 
GO  TO  72 
71  IBIRD=4 

AMACH=.82 

XX(9)=XX(76) -ATRIB(4) 

CALL  FLYME(XX(9) ,AMACH,XX (42) , lOPT ,V, IBIRD) 
XX(9) =XX(9) -V 
USERF=V*XX (43) 

RETURN 

Fig.  F-14.  USERF  (14) — Tanker  Cruise  Fuel 


C  LOITER  FUEL  AIRLIFTER 
15  IOPT=5 

ATIME=TNOW-XX(21) 

IF (TNOW-XX (21) ) 50,50,54 
50  USERF=0 
RETURN 
54  IBIRD=1 

IF(ATRIB(2) .EQ.2. .OR.ATRIB(2) .EQ.4.)  TBIRD=2 

CALL  FLYME(XX(10) , ATIME , XX ( 3 2 ) , lOPT , V, IBIRD) 

USERF=ATIME*V 

XX (10) =XX ( 10) -USERF 

RETURN 

Fig.  F-15.  USERF  (15) — Loiter  Fuel  Airlifter 
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C  LOITER  FUEL  TANKER 
16  IOPT=5 

ATIME=TNOW-XX{22) 

IF(TNOW-XX (22) >52,52,53 

52  USERF=0. 

RETURN 

53  IBIRD=3 

IF(ATRIB(2) .EQ.l. .OR.ATRIB(2) . EQ . 2 . ) IBIRD=4 
76  CALL  FLYME (XX(9) ,ATIME,XX(42) ,IOPT,V,IBIRD) 
USERF=ATIME*V 
XX(9)=XX(9) -USERF 
RETURN 

Fig.  F-16.  USERF  (16) — Loiter  Fuel  Tanker 


C  ABORT  CALCULATION  AIRLIFTER 
17  USERF=XX(8) -ATRIB(4)-XX(23) 

RETURN 

Fig.  F-17.  USERF  (17) — Abort  Calculation  Airlifter 


C  ABORT  CALC  TANKER 
18  USERF=XX(7) -ATRIB(4) -XX(24) 

RETURN 

Fig.  F-18.  USERF  (18) — Abort  Calculation  Tanker 


C  REFUELING  TRACK  TIME 
19  AMACH=.62 

IF(ATRIB(2) .EQ.2. .OR.ATRIB(2) . EQ . 4 . ) AMACH= . 74 
RALT=25.0 

TASTRK=TASP (AMACH ,RALT) 
USERF=250./(TASTRK+RNORM(.9, .16,1) ) 

RETURN 

Fig.  F-19.  USERF  (19) — Refueling  Track  Time 


C  TANKER  TRACK  FUEL 
20  R=.27 

IF(ATRIB(2) .NE.3.0R.ATRIB(2) .NE.4.)R=.25 
USERF=ATRIB(3) *R 
XX (9) =XX (9) -USERF 
RETURN 


Fig.  F-20.  USERF  (20) — Tanker  Track  Fuel 


AIRLIFTER  TRACK  FUEL 
R=.26 

IF(ATRIB(2) .NE.2. .0R.ATRIB(2) .NE.4)R=,45 
USERF=R*ATRIB(3) 

XX ( 10 ) =XX ( 10 ) -USERF 
RETURN 


Fig.  F-21.  USERF  (21) — Airlifter  Track  Fuel 


TANKER  TIME  TO  DESTINATION 
IOPT=l 
IBIRD=3 

IF(ATRIB(2) .EQ.1.0R.ATRIB(2) .EQ.2.)  IBIRD=4 
XX(9)  =  (XX(9) -XX(6)+XX(95)+XX(46)  ) /2 
XX(70)=CALT(XX(9) , IBIRD , XX ( 8 2 ) ) 

CALL  FLYME(XX(9) ,XX(93) ,XX{70) , lOPT , CT , IBIRD) 

IOPT=2 

CALL  FLYME(XX{9) , XX ( 93 ) , XX ( 70) , lOPT, CD, IBIRD) 

IOPT=3 

CALL  FLYME(XX(9) ,XX(93) ,XX(70) ,IOPT,XX(60) , IBIRD) 
AMACH=.82 

IF(ATRIB(2) .EQ.1.0R.ATRIB(2) .EQ.2.)GO  TO  86 

IBIRD=3 

GO  TO  87 

IBIRD=4 

XX(9)=XX(9) -XX(60) 

CRDIST^XX (5) -CD 

XX(62)=  CRDIST/(  (TASP  (AMACH,  XX  (  70)  )-RNORM(.83,  .16,1) )  )+15. 

USERF=XX(62) +CT 

RETURN 


I 


Fig.  F-22.  USERF  (22) — Tanker  Time  to  Destination 


C  AIRLIFT  TIME  TO  DESTINATION 
23  IOPT=l 

IBIRD=1 

IF(ATRIB(2) .EQ.2. .OR.ATRIB(2) ,EQ,4.) IBIRD=2 
95  XX(10)  =  (XX(10)+XX(6)+XX(94)+XX(45)  )  /  2. 

CALL  FLYME(XX(10) ,XX(93) ,XX(73) , lOPT , CT , IBIRD) 

IOPT=2 

CALL  FLYME(XX(10) ,XX(93) ,XX(73) , lOPT, CD , IBIRD) 

IOPT=3 

CALL  FLYME(XX(10) ,XX (93) ,XX ( 73) ,IOPT,XX(61) , IBIRD) 
AMACH= . 7  7 

IF ( IBIRD . EQ . 2 ) AMACH= . 74 
97  XX (10) =XX (10) -XX(61) 

CRDIST=XX(2) -CD 

XX(63)=  CRDIST/  (  (TASP  ( AMACH ,  XX  (  73 )  ) +RNORM(  .83,  ,16,1) )  )+15. 
S=TASP (AMACH, XX (73) ) 

USERF=XX(63) +  CT 
RETURN 

Fig.  F-23.  USERF  (23) — Airlift  Time  to  Destination 


C  FUEL  CRUISE  &  LANDING  TANKER 
24  IOPT=4 

AMACH=.82 

IF(ATRIB(2) .EQ.1.0R.ATRIB(2) .EQ.2.)GO  TO  150 

IBIRD=3 

AAL=5. 

GO  TO  151 

150  IBIRD=4 
AAL=2. 4 

151  CALL  FLYME(XX(9) ,AMACH,XX(70) ,IOPT,V,IBIRD) 

R=V 

USERF=(XX(62) -15) *R+AAL+XX ( 60) 

RETURN 

Fig.  F-24.  USERF  (24) — Fuel  Cruise  &  Landing  Tanker 


120 


C  FUEL  &  LND  AIRLIFTER 
25  I0PT=4 

IF(ATRIB(2) ,EQ.2.0R.ATRIB(2) .EQ.4,)G0  TO  152 

AMACH=  .77 

AAL=5.2 

IBIRD=1 

GO  TO  153 

152  AMACH=.74 
AAL=2 . 5 
IBIRD=2 

153  CALL  FLYME(XX(10) ,AMACH,XX(73) ,IOPT,V,IBIRD) 

R=V 

USERF==  (XX  (63)  -15)  *R+AAL+XX  (61) 

RETURN 

Fig.  F-25.  USERF  (25)— Fuel  Cruise  &  Landing  Airlifter 


START,  TAXI,  TAKEOFF  FUEL  CONSUMPTION  VALUES 


Aircraft 

C-5 

C-141 

KC-10 

KC-135 


STTO  (lbs) 
2800 
1900 
3000 
2000 


GROUND  IDLE  FUEL  CONSUMPTION  RATES 

Aircraft  RATE 


( Ibs/min) 


C-5  120 
C-141  80 
KC-10  120 
KC-135  100 


AIR-REFUELING  FUEL  CONSUMPTION  RATES 
Aircraft  RATE 


{ Ibs/min) 
( Ibs/min) 


C-5  450 
C-141  260 
KC-10  270 
KC-135  250 


\ 


Fig.  F-26 .  Aircraft  Fuel  Consumption  Values 
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Ai'7HKn'>  A1T/>I10V  >*v^0tz 


ION 


•  2h? 


FLYME  is  the  control  program  for  a  set  of  sub¬ 
routines  which  provides  the  fuel  data  for  both  the  flight 
planning  and  simulation  models.  The  fuel  data  bases  for 
the  C-5A,  C-141B,  and  KC-135  aircraft  were  provided  by  the 
Air  Force  Global  Weather  Central  and  is  the  same  data  used 
to  generate  operational  computer  flight  plan  fuel  require¬ 
ments.  Fuel  data  for  the  KC-10  was  available  from  the 
KC-10  System  Project  Office,  AFSC,  WPAFB,  OH  in  the  form 
of  experimental  performance  charts.  Data  points  from  these 
charts  were  compiled  into  FLYME  to  permit  the  same  sub¬ 
routine  calls  to  address  all  four  aircraft. 

FLYME  transfers  control  to  one  of  the  four  air¬ 
craft  siabroutines  labeled  C-5A,  C-141B,  KC-10,  and  KC-135. 
The  aircraft  subroutine's  coefficients,  in  dimensioned 
variables,  compute  the  feel  data  from  six-degree  polynomi¬ 
als  evaluated  in  subroutines  POLLY  and  POLY.  Five  options 
are  available  for  each  aircraft.  Option  one  computes  time 
in  the  climb  to  cruise  altitude  in  minutes.  Option  two 
computes  the  ground  distance,  in  nautical  miles,  traveled 
in  the  climb.  Option  three  computes  the  fuel  consumed  in 
the  climb  in  1000s  of  lbs.  Option  four  computes  the 
cruise  fuel  rate  in  1000s  lbs  per  minute.  Option  five  com¬ 
putes  the  rendezvous  loiter  fuel  in  1000s  lbs  per  minute. 
Each  option  contains  input  and  output  parameter  checks  and 
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error  statements  printed  by  subroutine  EHKors.  The  data 
is  restricted  to  cruise  altitudes  between  25,000  and 
41,000  feet  and  a  standard  day  temperature  with  a  tempera¬ 
ture  deviation  of  zero. 

The  C-141B  data  is  the  same  data  used  for  the  E-3A 
aircraft  with  a  3  percent  degrade  in  all  fuel  computations. 
The  current  C-141B  fuel  curves  have  not  been  prepared  by 
AFGWS  at  the  time  this  data  was  obtained.  However,  the 
data  provided  with  the  3  percent  degraded  data  matches  the 
performance  charts  in  the  C-141B  technical  orders. 

FLYME  calls  two  other  subroutines,  CALT  and  TASP . 
GALT  returns  the  highest  altitude  at  which  an  aircraft 
can  cruise  for  a  given  gross  weight,  and  Air  Traffic  Con¬ 
trol  (ATC)  hemispherical  altitude  structure.  The  calcula¬ 
tions  are  derived  for  a  given  gross  weight  linear  regres¬ 
sion  of  the  performance  ceiling  charts  for  each  aircraft. 

TASP  returns  the  true  airspeed  of  the  aircraft  at 
the  given  cruise  altitude  and  mach  .u.iber  derived  from 
the  equation; 

TAS  =  MACH  NUMBER  *  SPEED  OF  SOUND  AT  ALTITUDE 

The  listing  of  the  computer  code  for  these  sub¬ 
routines  follows: 
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SUBROUTINE  FUHE  (Vl.V2tV3.I0PT.V.IBIRDl 
IF  IIBIRB  .CE.  I  .AND.  IBIRD  .LE.  41  00  TO  5» 
CALL  ERROR  ("IBIRD  ERfiO".‘R  IN  FLIKE"! 

RETURN 

Sf  CO  TO  (liB)2Mi3F4i4e0)>IBIRD 
III  CALL  C5A  (VhV2.V3.10PT.V) 

RETURN 

211  CALL  CUIB  (VliV2iV3iI0PTiVl 
RETURN 

311  CALL  KCIIA  (Vl.V2iV3.10PIiVl 
RETURN 

411  CALL  KCt35A  (VlfV2»V3.10PTiVl 
RETURN 
END 
C 
C 
C 

C  RETURNS  TRUE  AIRSPEED  IN  NN/  BIN 
FUNCTION  TASPIANACH.ALTl 
DIMENSION  SSND(7) 

DATA  SSNO/9.563331  9.61167i  9.69667» 
t  9.783331  9.84833.  9.95333.  11.13447/ 

DO  II  I  :  1.  7 
J  *  I 

IF  {ALT.CE.(39.-2*I)>  CO  TO  21 
II  CONTINUE 
21  TASP=ANACH*SSNO(J» 

RETURN 

END 

C 

FUNCTIM  CALTICNT.IBIRD.TCI 
DIMENSION  ALPKA(4).BETA(4) 

DATA  ALPHA/41 . I 2378.41 . 79153.54 . 58588 . 41. 13495/ 

DATA  BETA/-.e45S554.-.f95f7B2.-, 6435342. -.1073148/ 
HIALT  =  ALPHAdBlRDI  ♦  EETAIIBIRDI  ♦  CWT 
DO  tl  I  M.  4 
ALT  --  41.  -  4*1 

IF  (TC  .CE.  3.141591  ALT  =  43-  4»I 
IF  (HIALT. CE.ALT!  CO  TO  26 
II  CONTINUE 
21  CALT  -  ALT 
RETURN 
END 
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SUBROUTINE  ERROR  (HSCtNSUB) 

PRINT*, NSC. NSUB 
STOP 
END 
C 

c 

c 

FUNCTION  POLriC.H.I) 

DINENSION  C(N) 

SUN  --  «.« 

DC  1«  IH^C.H 
J=H-IH+2 

If  SUN  :  r  *  (C(JI  4  SUHI 

poLt  =  c(i)  4  sun 

RETURN 

END 

C 

C 

C 

FUNCTION  POLLTIC.K.L.I.n 
DIMENSION  CIK.L) 

sun  :  f.f 

DO  Zf  in^ZiL 
JO  s  L-IH42 

zi  sun  =  T  ♦  (POLTicii.jji.K.m  4  sum 
POLLT  5  POLTICd.n.K.n  4  sun 
RETURN 
END 

c 

c 

c 

SUBROUTINE  C5ft  (V1,VZ,V3,I0PT,V1 
DIMENSION  R15CTI (5,51 ,A15CTZ(5,5) ,Rt5CFI (5,5) 
DIMENSION  AISCRt (5.5) ,A15CRZ(5,5) ,A15CFZ(5,5) 
DtnENSION  A15E1I0C),A15E1J13S) 

DIMENSION  A)5EIE(3C) .AISEIFISS) ,AI5E1C(3S) ,A15ElH(3t) 
DIMENSION  AI5ElA(3S),Al5ElB13i),A15EiC(3t),A15ElD(3Sl 
DIMENSION  A15tl(34.1f),PdEl(2e),  A15H(9,3) 
EOUIVftllNCE  (S15E1A,A15'1(1,1)),  (A15EIB,A15E1(1,Z)) 
EQUIVALENCE  '.IISEIC-Al?' ;  diS!),  (A15E10>A(5E1(1,4)) 
EQUlVALEii’t  ;v.5ElE,A;::'(N5i',,  (A15lIFiA'i5E'.('.4)I 
EQUIVALENCE  ir.'.5ElC,Al5.;i(I,7)),  (A15£1H.h15E1  (1,8)) 
EQUIVALENCE  ISIEEII,A15EI(1,9)),  (AI5£IJ,Al5EUl,lfll 
C 
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DIST  FOUR  ENCIME  COEFFICIENTS 


C 

C 


t 

i 

i 

i 

t 

t 
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t 
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i 

t 
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i 

i 

i 

t 
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I 

I 

t 

t 
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DATA  A15CT1/  .57427a7E402. 

.It7S5422E-65. 
.87JS4921E-DI. 
-.i77s58t7£-10, 
.65:if5442E-i’5. 
.243ZC33CE-?.. 
-.94273221£-r-. 
.454489tlE-§5. 
-.fc33!H522E-14/ 
DATA  A15CT2/  .547542l9E+02f 
-.l3()4774fcE-63. 
.17g54432£*0t» 
-.155il7inE-Di» 
.34091933E-D2. 
-.34&842?6£+8Di 
•.8447i539E-04. 
-.23314340E-03. 
.2142309fE-l«/ 
DATA  AISCFW  .9621S572E4g2, 
.24597864E-0ii 
-.2e843g73E4fS, 
.3334l333E-«9» 
-.84854549E-f4i 
.152«1547E+0g. 
•.48l«55«g£-§8i 
.133J3874E-i4f 
•.14g}4llZE-l3/ 
DATA  Al5CF2/*.40738972E+02. 

*.47785562E'«5t 
.64395325E40gi 
-.583I«999E-D7i 
.11091388E-02i 
-.154440I9E400. 
.2l5t8g54E-06. 
-.4984t420E-04, 
.38g56788E-ll/ 
DATA  A15CRI/  .44897424E4f2. 

.15t43fc5tE-02. 

-.2831890lE40li 

,4l25207lE-05. 

-.10247921E-01, 

.20g«4468E+00, 

-.45438Z95E-05. 

.l74(lve8OE-03. 

-.1322l863£-i9/ 


.340?22'OE+2e 
-.523-a-E-C9 
-.1874^;;-;-03 
.ifets:'"  E-tei 
-.2t624i:;-'-.-08 
149594  Mt- 03 
.42U512;M2 
-.l4307733E-g7 

-.2307g833E401 

.2282C511E-08 

-.185787i9E-lt 

.1051518iE402: 

-.14248281E-04 

.U21gS24E-0l 

-.l484927IE-08i 

.25464057E-05 

.4Ce21259£-01 

-.t277t840E-09- 

.t0823261E-g3 

•.370t9U9E401 

.ie28g328E-0C 

-.1210fc921E-g2 

.l8i75t77E-n 

-.35575162E-f7 

.51338421E400 

.17370157E-07 

-.706S8844E-02 

.35193831E401 

-.5I409909E-05 

Mmmz-n, 

-.35559627E-09I 

.5120344lE-ft' 

.44392409E4I1 

-.ii759124E-l5i 

.80245119E-01- 

-.47324t32E401 

.n873t35E-f3 

-.15339955E-01 

.17480e55E-»7 

-.43fc874i3E-05i 


-.94784741E-03 

'.15142192E402 

.13194550E-04 

458849431-02 

.1S817445E-I2 

.44197544E-04 

-.42545431E-I3 

.17444845E-10 

.28947473E-01 

-.57847f84E402 

.885851gIE-l4 

-.31197773E400 

.28528g73E-07 

-.17712504E-i3 

.787U422E-02 

-.r2152450E-07 

-.14874350E-03 

.25493794E402 

-.74947174E-f4 

.29754803E-fl 

-.44437003E-10 

.33924581E-D5 

-.17377349E-02 

.38945935E-10 

-.74424177E-03 
-.215I4840E402 
.33459370E-04 
-.10334714E+09 
.84739747E-08 
4723981 4E-04 
.17443395E-02 
-.22595754E-08 

-.12549410E400 

.352628i2E4«2 

-.94009847E-03 

.372294UE40fl 

-.48424332E-04 

.4l02e842E-03 

-.23448958E-82 

.43437131E-07 
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c 

c 

c 

c 

c 

c 

c 


DATA  ftl5CR2/-.lM85249E*C2i  -.5gfc88M3E+91,  .?1317535E-ei. 
i  -.47085761E-e3i  .ee535574E-86i  .48::e8850E+02. 

&  -.18i83532E+81.  .24188‘,45E-»1 .  -.13!J520S'E-03t 

i  .274'}4952E-e4i  .4!73n72E+Sl.  -.43913t3tE-01i 

t  .118«4»94E-07i  .37fc93tnE-05.  -.13293508E-07, 

t  -.42739199E480.  .9t429417E-e2.  -.48830744E-04. 

i  .12779489E-et.  .13743494E-09,  .t7et4995E-e2t 

i  -.14844954E-03i  .99727724£-0t,  -.14487723E-08» 

&  -.35185374E-11/ 

0)NSTANT  ALTITUDE  CRUISE  COEFFICIENTS 

A=23«  B=25H  C=27H  D=29N  E=31« 

F=33H  C=35«  H^37H  I=39H  J=41M 


DATA  A15E1A/ 

t  .34542733E402I-.45837S38E+03 
i  .12Ui25tE+05.*.5l979964E+04 
t  .8295926tE+01,-.10881802E+02, 
A-.912201»72E'04.-.88490494E-03 
(  .00000000E400.  ,00000000E400i 
i  .4ifl2l}03E'04.  ,00000000E400 
i-.U850194E-09,  .n5l6fc74E-09i 
i  .0e000t40iE400.  .tf0000000E400 
i  .00000000E400*  .000000f0E400i 
DATA  At5ElB/ 

i  . 3692339iE+03 i - .34399470E+04 
t  .Z5433474E+05.-.9l34423tE434i 
t  .43409:}47E+8l.-.38744322E+ini 
t-.40998729E-03.  .15584510E-02- 
i  .00000000E400>  .|)0000000E400 
t-.47227338E-04i  .f?3000000E+fl0 
t-.17424758E-09t  .20044t02E-«9 
t  .00000000E400.  .fl0000000E400< 
t  .00000I400E400)  .00000000E4001 
DATA  A15E1C/ 

(  .45587749E+03.-.5970l378E+04i 
t  .37U2890E405t-.t2793483E405i 
t  .37259380E401,-,t832ei88E+aii 
t- .  47049495E-  03 .  .  33  ?  ■'  V749E-02 
t  .f0000?79E+;;:4-  :;-;00£+00 
t  .3520259CE-i'i.  .,’j:::3'0E+a0 
(  .53t50935E-e?.-.iE:931S4E-«9i 
t  .00000a00E40».  .r00003r0E4a0< 
t  .100000008400.  ..i0000000E400i 


.35939765E404.-. 101846848405, 
.24923058E400,-.24911629E401, 
.49621873E401,  . 000000008400, 
.21114a37E-02,-.12173669E-02, 
. 519042788-06, -.84674241E-04, 
.000000008400,  .00000000E400, 
.000000008400,  .000000008409, 
.134635938-13,  .000000008400, 
.000000008400,  .000000008400/ 

.139341068405, -.272295358405, 
.212296358400, -.174767498401, 
.933329138400,  .000000008400, 
-.284407248-02,  .195803868-02, 
.179678238-06,  ,13420929E-06, 
.000000008400,  .000000008400, 
.000000008400,  .000000008+00, 
.136714328-13,  .:J00?0000E+00, 
.000000008400,  ,000000008400/ 

.227271698+05, -.418335708+05, 
.270397108+00, -.209087948401, 
-.76937575E-01,  .003000008+00, 
-.542141548-02,  .255319338-02, 
-.480638048-06, -.526788478-07, 
.000000008400,  .003000008+00, 
.000000008400,  .000000008+00, 
.176158178-12,  .r3300000E+00, 
.000000008+00,  .001000008400/ 
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DATA  A15E1D/ 

t  .92t»4244E+«3.-.84«ll76E+94 
t  .48e0347tE-t«5.-.U13838Sr+tf5i 
i  .3i«59(I84E+01 .  15i03997E+0l 
t-.246l7458E-02,  .5473t937E-g2. 
i  .002000008408.  ,00000000E400 
i  .436945Z7£-0t.  .«e20t020£4fl0 
i-.14378'?57£-08.  .9t07974££-09 
t  .80000000£400)  .00002000£400, 
t  .0000e:i00E400,  .00000000E408, 
DATA  A15EtE/ 

t  .U539089£464.-.149678t8E+S5i 
t  .747t0864E405.-.25t57140£+85i 
t  .10789359£402,-..42977e36E4«l, 
fc-.28883011£-i2.  .7028l398£-a2i 
(  .00000880E4e8)  .(}e000000E400 
l-.l5449815£-e5.  .ii0063000E4a0 
t-.298l2734£-08.  .59255904E'a9, 
i  .ee000a00£480i  .fi-3009f00£480, 
&  .80800a00E408.  .0230000fE400i 
DATA  A15£1F/ 

&  .20139264£404,-.187@224iE405i 
(  .99865248E405.-.327«3lie£+a5i 
i  .28932324£402.-.2173283lE4i;2i 
t-.28542462£-«2.  ,985A5l40E'«2i 
t  .08000a80E480i  .08000000E400 
t-.lSlS1388£-i5)  .l!0000000E480 
t-.2744e:373E-88.  .3478577S£'g9. 
(  .80000a00E4f8.  .00e00000E400i 
(  .f8000880E480.  .80000000E40I, 
DATA  Ai5ElC/ 

(  .28498(?53E404.-.26945449E405i 
(  .13863563E486.-.44537107E405i 
t  . 5 1 758489E482 , - . 34757t5:E402 1 
4.-.4859M99E-02I  .17735313E-01' 
t  .fe008088E48ei  ,00000000E408i 
t- .59440357E-05.  .03C5f000E400: 
1-.11842«27E-I7.  .38334l88E-«9i 
(  .8fe00008E4ff!  ,00000000£400< 
t  .8000C03?E4i0i  .00000§f0E4a0i 
DATA  Al^ElH/ 

&  .39282855E42J,-.344{3444E-;Ti 
i  .I7445';8i;£40!.-,59253t3:£4:j: 
t  .77245a9l£+i’2.-  ■•2t484nE4rf2i 
t-.  18107550E-02 .  - .  .3i95JSb£-'3l  i 
t  .88008a08E40e,  '  .t000E4a8< 
t  .59782m£-04.  ..ij200000E4«8, 


i  .31547536£405.-.55822944£4fl5. 

.£58191t4E4iii,-.2729459SE401, 
-.25354725£4e0,  .02308023E4a0, 
-.5705423i£-02.  .Z4719t5iE-02. 
■  .249£4545E-3TT-.2t343423E-;J5. 
^  .e:-:;;'7£402, 

,  .0cii;;2JgE4e0.  .S03e£:^t2£490, 
.24ie3304£-12.  .022gC00»E490, 
.00000200£400,  .00200000£4a0/ 

.5305P.558£4e5.-.9l498270E4g5. 
.147t4374£4gl,-.416070?4£4ai, 
.3109i440£408.  ,00080000£400. 
-.1187fc832£-0l.  .68525515E-32. 
.2t828714E-05.  .n51t27tE-ii5. 
.00000a00E400.  .00000000E430I 
.00000J00E400,  .a0000000E4a0, 
.94fc44581£-12.  .0C003000£430, 
.80000a00E400.  .00000200£400/ 

.4772935i£4g5.-.lt818523£4(}4T 
.28423329£481,-.150072t8E432. 
.505837y2£401,  .a200000g£400, 
-,t74l7374E-81.  .95731841E-92. 
•  .18245041E-05.  .1956S27t£-05. 
.003330«0£4g0.  .g3.\“r,0SIE+«0. 
.000f30T0E4ef,  .§5fJ2200E40g, 
.99517707E-12.  .80230g30£400, 
.00080028E48I) 

.948848:i£405.-. 1470387gE+«4. 
.57t494'8E48l,-.27435599E492. 
.7129R942E401,  .0 J320008E498, 
-.3546e957E-ei.  .23999188E-01. 
.73997283E-85,  .85733339E-85. 
.000000i0£40g,  .C§3?30?iE4«l, 
Miizmimy  .g::j?3!;0E490, 
.5397e737£-lh  .j;2J2;:;E49f, 
.fv;S00:0E4g0,  .i0gT00g0E+8f/ 

.t2e3l914E404.-.23f.31798E484» 
.5044712;E40l,-.::?.:7t04E+92. 
.47544480E402,  .■';;]000g£4«8, 
.5S299573E-91 , - .55734284E-81 . 
.J4454525E-04.-.32033487E-04. 
,80080a00£480i  .00k.)e080£408i 
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t-.385i98ii3E-«8i-.3303482e£-B7. 

t,  .M8«0M8E+e8»  .a«e8«08«E+M.  .1339£i4v£-l?.  MimiZMh 

&  .eeeemE^ei)  .(i0888888Et3«>  Aimmimi 

DATA  A15En/ 

i  .4739«746E+8«r-.43B5154E+J5.  .15327383E+3ii- .2t48Z74ZE+84» 
i  .Z2l522e8£+fiA.-.72523l4D£-»SS.  .97e223S4E*?l i-,4?384915E+T2. 
i  .13869259£403i-.12736934£+83.  .45499131£402i 
i  .65em59E-e3f-.l^7l933l8£-^}^t^.U^93ll7E-Sl,-.lm7^t7E-^^f 
(  .e3088l{88E488i  ,80g8808eE4af.  ,70382(r38E-e5i  .27444919E-84) 
t  .17727124E-04>  .802e0CeeE480.  .e808e088E4e8)  .80008088E488, 
t-.233579A4E-l7.-.32i3m0£-87f  .0f080;)08E400,  ,8fi2.3gg25£480, 
i  .esmumm,  .meemE-tih  .24850308E-181  .8020i'000E488) 
i  MmmEmt  .messm*if0i  Mteig^eE^es,  jgsmgeEH)^/ 
DATA  A15E1J/ 

<.-.13825tf5fc£403.  .t2171563E4S4i-.19473890E4e4,  .92843448E483. 
(  .M«i0Mg£4g»,  .W0eggf8£4a0,-.97332943E4e0.  .lt383971E4fllf 
i-.?543e:3e3£480>  .8888000i£488f  .D888aa80£4gg,  .88ge8008E4a0i 
t  .59334393E-83.  .14846428E-a3.  .8W0g300£400,  .800g0g08£4a0, 
(  jegfemE-ng,  .00mg§§E*08f.u797iziE-si,  jgggsmEm, 
t  .00E08;>08E4fet  .0g82008«E480,  .e0888a88E4e0(  .80800880E480I 
&  MmmEHS,  .900gemE488«  .8Mi}0a88£408,  .ae80i-000E4a0, 
i  .memsEmt  jmgsmiSit  .gmmeEm,  .emuEiSE^sti 
i  .M00»mE*g8,  jsmmE*f9f  .dsemigEm,  .6smm£*»t/ 
c 

DATA  P15Et/  23.li  25.1.  27.1.  29.0.  31.0. 
i  33.0.  3$.0t  37.0.  39.8)  41.0/ 

C 

DATA  At5H/ 

(  .180001  .18187.  .18887.  .20187.  .20333.  .28833. 

I  .22333.  .22333.  .23888,  .24500,  .24847.  ,25500. 

t  .28887,  .27038.  .28187,  .28833,  .29333,  .99999, 

t  .31008,  .32038,  .99999,  .33333,  .38030,  .99999, 

(  .38008,  .99999,  .99999/ 

C 

C 

IF  (lOPT  .GE.  1  .OR.  lOPT  .LE.  51  CO  TO  450 
CALL  ERROR  ("lOPT  ERROR","  IN  CSA"! 

RETURN 

C 

C 

450  CO  TO  (  500,800,700.880.9001.  lOPT 
C 
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C5A  CLIKB  TIME.  RANGE.  FUEL  COMPUTATIONS 


C 

C 
C 

C  Vl=CWT  V2=IN1T  ALT  y3=TERM  ALT 
C 

C  AT  5M  OUTPUT  V  =  TIKE  TO  CLIMD  (MINUTES) 

C 

SH  GUS=AMAn(Vii354.) 

PAS=V2 

TD=20. 

PAT--y3 

PAS^AKAJKPASiB.) 

CWT=CWS- ( {.0ei25tCMS-.2375)4PAT) 

IF  ((CHT.lT.335.).OR.(CMT.CT.770.)  .OR. 
i  (C«S.LT.35tf.l.OR.(C«S.CT.77e.)  .OR. 
i  (PAT.LE.PAS))  CO  TO  11111 
TIMES  =  t. 

IF  (PAS.EO.e.)  CO  TO  5S6 

TIKES=POLLI (A15CT2.5.5.P0LLt (A15CT1 .5.5.CNS.PAS) .TD) 
TIMES=AMAn(TIKES.2.5) 

550  IF  (CMT.LT.  i90.)  CO  TO  575 
C«T=CKT-130. 

575  T£STC«=CHT 

TEMP=P0LLT1A15CT1.5.5.CUT.PAT) 

TIKET=P0LLr(A15CT2.5.5.TEMP.TD) 

TIK£T=A«A1(1(TIMET.4.5) 

V=TIMET-TIMES 

IF  (IV.LT.0.).OR.(V.CT.55))  CO  TO  22222 
RETURN 
C 
C 

C  AT  600  OUTPUT  V  s  CLIMB  RANGE  (NMSI 
C 

600  CUS^AMAlKVl/ie.  35.4) 

PAS=V2 

TD=20. 

PAT=V3 

PAS--AMAI1(PAS.0.) 

CNT=CNS-((.0125«CIJS-.2375)K.UPAT)) 

IF  IICWT.LT.33.5).0R.(GHT.CT.77.).0R. 
t  (CNS.LT.35.)  .OR.iGHS.GT.77.).OR. 

t  (PAT. LE. Phi))  CO  TO  mil 

RANCES^0. 

IF  (PAS.EQ.  e.f)  CO  TO  o50 

RANGES=F0LLttH!5CR2.5.5.P0LLT(A15CR1.5.5.GNT.PAS).TD) 

RANCES=AMAr.(F:A.iCES,25.) 

650  TEMP=P0LLr(A15CR1.5.5,GNI.PAT) 
RANCET=P0LLT(A15CR2.5.5.TEMP.TD1 
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V:RANCET-RAKCES 

IF  (tv.LT.i.i.oR.iy.CT.ssen  co  to  iim 
RETURN 
C 
C 

C  AT  m  OUTPUT  V  =  CLIMB  FUEL  limt  OF  LBS> 

C 

709  CUS=AMAll(Vli3S4.) 

PAS=VZ 

TD=Z9. 

PAT=V3 

PAS--AHAn(PASi9.) 

G«T=CUS-  ( ( .091 25*C«S-.  23751  iPAT) 

IF  ( ICUT.LT.335.) .OR. (CMT.GT.770. 1 .OR. 
i  ICUS . LT . 359 . I . OR . (CMS . CT . 779 . 1 . OR . 

t  iPAT.LE.PASI)  CO  TO  Hill 
FUELS=B. 

IF  (PAS.EQ.  9.91  CO  TO  759 

FUELS=POLLr (A15CF2i5f5i (POLL! (A15CF J .5.5.GHS.PAS1 1 »TDI 
FUELS:ABAn(Fl£LS.2.51 
759  IF  (CMT.LT.  (.91.)  CO  TO  775 
CUT=CHT-51. 

775  TESTCM=CHT 

TEMP=P0LL!(A15CFl.5i5.GMT.PAT) 

FU£LTmLr(A15Cf2t5.5i(POLL!(A15CFl>5i5.CWSfPAT))»TI» 

FUELT--AnAIl(FUELT.3.5t 

V-FUELT-FUELS 

IF  (IV.LT.9.).0R.1V.CT.55.))  GO  TO  22222 
RETURN 
C 

C 

C  C5A  CONSTANT  ALTITUDE  CRUISE 
C 

C  VhCUT  V2=MACH  V3=ALT  OUTPUT  V  =  1909  LBS  /  WIN 
C 

890  CM  =  AMA71(Vli354.) 

AMACH  :  V2 
PA  =  V3 

IF  (PA.LT.0.9  .OR.  PA.CT.45.  .OR. 

(  AMACH.lt. 6. Z  .OR.  AHACH.CT.9.9  .OR. 

(  CH.LT.359.  .OR.  CH.GT.775.) 

«<  CO  TO  iiltl 
DO  859  J  =  2.  10 
I  =  J-1 

IF  (PA.LT.  PlSEUjn  CO  TO  875 
859  CONTINUE 

875  Tl=P0LL!(Al5El(ltn.6.6.AMACH.CHI 
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TZsP0LmA15Ei(lfI+n.t,iiAHACHiCHI 
V=T1  +  IPA-PISEKD)  »  (T2-T1)/!P15E1II+1)-P15E1(I)) 
IF  (V.LT.9.  .OR.  V.CT.50.)  CO  TO  lllll 
V  =  TASP(A«ACHiPAI  /  V 
RETURN 
C 
C 

C  C5A  RENDEZVOUS  HOLDiNC  FUEL 

C 

C  VUCWT  V3-ALT 

C  OUTPUT  V  =  1001  LBS/HIN 

C 

900  CUT  :  AnAllliVh354.) 

PA  =  V3 

IF  (CUT.LT.35i.  .OR.  CUT. GT. 775.  .OR. 

&  PA.LT.f.e  .OR.  PA.CT.45.)  CO  TO  lllll 
DO  920  I  :  If  9 
J  -  I 

IF  (CUT  .LE.  (3Z0.4IM0.))  CO  TO  940 
920  CONTINt£ 

940  00  960  I  :  li  3 

K  -•  I 

IF  (PA  .LE.  (25.+IM.11  CO  TO  980 
960  CONTINUE 
980  V  :  AiSH(UiKt 

IF  (V  .LT.  0,  .OR.  V  .CE.  .999991  CO  TO  22222 
RETURN 

lllll  PRINT!."  lOPT  ='MOPT."  CUT  =".V1 
CALL  ERROR  (“  INPUT"."  C5A") 

RETURN 

22222  CALL  ERROR  ("OUTPUT  ERR". "OR  IN  C5A"I 
RETURN 
END 
C 
C 

Ct«««**! 

c 

c 

c 

SUBROUTINE  Ct413  (V1.V2.V3.10PT.V) 

C 

C 

OINENSION  C1A(1.:.C2A(16I.C3A(16I 
DINENSION  CA(!6-3).El(16.8).Fl(8)iHFFi8) 

DINENSION  ElAllv) . £18(16) .EtC(16).ElDll6) 

DINENSION  E1E(!10>E1F(16).£!CI16).E1HI16) 

DINENSION  018(2). 028(2) .D38IZ)tDE(2i3) 


c 

EQUIVALENCE  (CiAtCA(iil))  >  (CZAiCAd.Z)) 

EQUIVALENCE  (C3AiCA(l>3) )  i  (DE(lil) fDlE! 

EQUIVALENCE  ifcBd i2) iD2£)  .  (DE(h3).D3E) 

EQUIVALENCE  iElAiEHiil))  .  (ElB.EKliZI) 

EQUIVALENCE  (ElC.El(li3)l  t  (ElDiEl(li4tl 
EQUIVALENCE  (ElEtEl(li5))  •  (ElF.Elil 16) ) 

EQUIVALENCE  (EtCiEl(l>7))  r  (ElHiEKliEI) 

C 

C  CLINE  COEFFICIENTS 
C 

DATA  CIA/ 

i  .4180952tE+02.-.53315043E+g8.  .2187t803E-e2i-.275t5421E-85, 
i-.8fc7S4842E+81i  .ni35442E+38i-.43t08378E-83.  .577l8t97E-8t. 
&  .44696824E+80.-.5i37g250E-82.  .22183138E-e4.-.29584408E-§7. 
Q-.t648S989£-g2.  .848902e4E-84.-.33228785E-0t.  .460l8t32E-89/ 
C 

DATA  C2A/ 

i  .10035453E+03t-.9452l497E+S8.  .19945408E-02.-.3894fc845E-8t» 
fc- .  19228782E+02.  .  17217714E+153.- .2ii8294BE-03. - .  10t45l22£-88, 
i  .79824788E+g0.-.55805l73E-82i-.10t89811E-65.  .24885471E-!)7, 
i-.90l3l1<.lE-82.  .4537816gE-04.  .18272715E-04t-.51003i39E-89/ 
C 

DATA  C3A/ 

.5825741  iE-tf4t 

&  .10373707E+8l.-.l4134949E-8h  .m08771E-04.-.97339575E-87, 
t-.48l<1585E-0W  .700<,8203E-83.-.33188495E-85t  .4t838939E-88. 
(  .559#5559£-03t-.81442520£-85.  .3744fc303£-07.-.44544775E-10/ 
C 
C 

C  OUTPUT  PARAMETER  LIHITS 
C 

DATA  D1B.D2E>D3B/  0.t  40.>  0.)  240. i  0.i  13.4/ 

C 

C  CONSTANT  ALTITUDE  CRUISE  COEFFICIENTS 

C 

C  A=25N  B=27«  C=29(1  D=31H 

C  E=33M  F=35H  C=37n  H=39I1 

C 

DATA  mi 

t  .39I44847E+03,-.15574907E484.  .240t4898E+04.-.12901585E+84. 
i-.37444194i.+01.  .14477557E482.-.25888t47E+02.  .t25t80l7E+S2, 
t  .ll«9416E-Bli-.53777762E-fll.  .85958771E-01  .-.450492BlE-01i 
l-.853i4424E-85i  .45845470E-84.-.79l398g9E-04.  .43972fc23E-84/ 
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c 

DATA  mi 

I  .5U924l2E+«3.-.l<r723254E+04t  .iee33207E+C4.-.48?7i804E+l»3. 
t-.65l5t544E+01.  .23482t24E+32i-.2tlt2veSE+e2.  .85415492E+91. 
(  .2t8i3472E-01.-.19330544E+l}8i  .12390882E+00.-.45274434E-i}h 
l-.35491305E-04(  .138380t7E-«3i-.lt990271E-03.  .4449709tE-«4/ 
C 

DATA  EIC/ 

t  .28949557E403.-.495820e5E+93i  .41718559E+03.-.19434718E+S3. 
t-.50477514E+0li  .t7444305E+tf2f-.19722525E402.  .49444219E+4W 
i  .25822704E-0h-.102759e4E+«0i  .13254223E+00.-.55888790E-!M» 
4-.40238875E-04.  .l45l0t37E-(}3.-.22l85442E-03.  .98219017E-04/ 
C 

DATA  mi 

i  .2l070841E+04t-.90042572E404i  .13304971E+05.-.44442359E+04, 
t- .  29394(?97£+02.  .  I2850585E4)}3  .-.187907 1 2E+03 .  .  9 1 87984 1 E+92 . 
i  .13199537E+00.-.58534003E+00.  .84239993E+00.-.42290730E+00. 
i-.19004980E-03.  .84411525E-03.-.12452789E-02.  .4t091292E-83/ 
C 

DATA  ElE/ 

i  .94574555E+03.-.3U55804E+04,  .35877400E+04.-.14375898E+04. 
t-.14021388E+02.  .5757ie55E+!J2.-.48477408E+02.  .24970184E+02. 
i  .83720121E-01 .3i892524E4(}0,  .40158349E+00.- . 14752504E+fl0. 
fc-.140l9347£-83.  .54454304E-03.-.70079584E-03.  .29909214E-03/ 
C 

DATA  EIF/ 

t  .225l9(l03E+03.  .21778855E+93.-.l3n5424£+04.  .925l98l5E4a3. 
1-.44744943E+01.  .l4899445E492.*.58483454E48l.-.31432492E+tft. 
(  .42731702E-01.-.13992524E400,  .14344784E400.-.45329795E-ai. 
4-.77l9@904E-i4.  .27475383E-03.-.31878584E-03.  .12e8l834E-g3/ 
C 

DATA  Etc/ 

t,  .87808213E404.-.28374209E405.  .28571235E405.-.837240i4E4»4. 
4-.14980194E403.  .51387421E483.-.54433S28E403.  .19514249E483. 
t  .80444275E400,-.28944827E4«1.  .33813903E401,-.12724897E401, 
t-.13887352E-32.  .5t2Z4322E-02.-.41882045E-02.  .2439t29lE-02/ 
C 

DATA  EIH/ 

4-.5345lfl/7E484.  .3t549174E405.-.545e4409E405.  .3213e929E405. 
t  .61tt3514E402.-.?8l23027E483.  .71183330E403.-.41571353E483. 
t-.2l3709Z4E400,  .14413042E481,-.284t3717E40t.  .t7208428E4lM. 
4  .22027924E-03r-.:78:9y3£-02.  .37514410E-02.-.23425345E-02/ 
C 

DATA  Ft/  25..  27..  29..  3t..  33..  35..  37..  39./ 

C 

DATA  HFF/  .t2t47,  .t3333.  .t4333.  .t5447, 

(  .t4833.  .t7833.  .t9333i  .20833  / 
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IF  IIOPT  .CE.  I  .OR.  lOPT  .LE.  5)  CO  TO 
CALL  ERROR  C'lOPT  £RROR“i"  IN  C141E“I 
RETURN 
C 
C 

i»6a  GOTO  (U0E>iie0in0eii2efii3E*ii  iopt 

c 

c 

C  CUIB  aiHB  TIME.  RANGE.  FUEL  COMPUTATIONS 
C 

C  Vl=CUT  V3--INIT  ALT  V3=T£RM  ALT 
C 

C  OUTPUT  V  =  CLIMB  TIME  (MINI 

C 

C  J=2  OUTPUT  V  =  CLIMB  RANGE  (NH) 

C 

C  J=3  OUTPUT  V  =  CLIMB  FIEL  U000  LBS) 

C 

1100  GU  :  AMAlKVl.Ul.) 

PAS  =  V2 
PAT  s  V3 

CUT  •-  CU  -  ((PAT-PAS)».000943»GM) 

TO  i  20. 

J  »  IOPT 

IF  ICH  .LT.  141.  .OR.  CM  .CT.  340.  .OR. 
t  PAS  .LT.  0.0  .OR.  PAS  .CT.  45.  .OR. 

i  PAT  .LT.  0.0  .OR.  PAT  .CT.  45.1  CO  TO  11111 

T1  :  POLLTICAd.J).  4.  4.  CM.  PAS) 

T1  -  AHAXK  Tl.  0.0) 

TZ  =  POLLTICAd.J).  4.  4.  CUT.PAT) 

V  =  (T2  -  Tl) 

IF  (J.EQ.3)  V  -  .98  »  V  »  .97 

IF  (V  .LT.  DBd.J)  .OR.  V  .CT.  DBIZ.JI)  CO  TO  22222 
RETURN 
C 
C 

C  CONSTANT  ALTITUDE  CRUISE  COMPUTATIONS 
C 

C  Vl-CMT  V2=NACH0  V';5=ALT  OUTPUT  V  =  1000  LBS/MIN 
C 

1200  CRH  :  AMAllIVl.UI.) 

AMACH  =  V2 
PA  =  V3 

IF  (CRM  .LT.  I4t,  .OR.  CRM  .CT.  340.  .OR. 

t  AMACH  .LT.  .2  .OR.  AMACH  .CT.  .8  .OR. 

t  PA  .LT.  0.0  .OR.  PA  .CT,  45.)  CO  TO  1)111 


DO  1250  0  =  1.  e 
IF  (FlU)  .LE.  PA»  I  =  J 
1250  CONTINUE 

V  :  P0LLt(EUhI)i4i4tAHACHiCRUIM.e4 

IF  (V  .LT.  16.  .OR.  V  .CT.  £0.)  CO  TO  22222 

V  =  TASPIABACH.PAI  /  V  »  .97 
RETURN 

C 

C 

C  HOLDING  FUEL  AT  FL256 
C 

C  V1=C«T  V2=m  V3=ALT 
C 

C  OUTPUT  V  =  1066  LBS  /  NIN 
C 

130(1  CUT  :  ANAUlVltUl.) 

PA  =  V3 

IF  (CUT. LT. 141.  .OR.  CUT. CT. 346.  .OR. 
i  PA  .LT.  6.6  .OR.  PA  .CT.  45.)  CO  TO  11111 
DO  1320  1  =  b  8 
J  =  I 

IF  (CUT  .LE.  (146.41*26.))  CO  TO  1340 
1320  CONTINUE 
1340  CONTINUE 

V  =  HFFIU) 

IF  (V  .LT.  6.  .OR.  V  .CT.  .221  CO  TO  22222 
RETURN 

mil  PRINT*i“  lOPT  ^MOPTi"  CUT  ='SV1.“  PA  ="tV3 
CALL  ERROR  ("  INPUT'S"  C141B") 

RETURN 

22222  CALL  ERROR  ("OUTPUT  ERR". "OR  IN  C141"l 
RETURN 
END 
C 
C 
C 

C  l*tl**f**f*f***»»«4tf}*ft««*4**»ft****4***4*ltt**l***t*tt* 

c 

c 

SUBROUTINE  KC16A  (VI.V2.V3.ICPT.V) 

C 

DIMENSION  TKCL«(7.3).Z25(7,3) 

DIMENSION  TKCRU(4.7),TKH(14) 
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c 

C  KC10A  CLItIB  COEFnCtENTS 
C 

DATA  TKCLfl/  t7.3i  17.4i  i8.3t  19. At  Z0.8f  22. 2t  19. 2i 

t  118..  124..  127.5.  132. t  1U.3.  154.5.  130.5. 

I  8.8.  9.5.  10.3.  11. 1.  12.5.  13.7.  12.0/ 

C 

DATA  225/  7.5.  8.fc,  10.2.11.7,13.4.15.6.15.8. 
i  48.5.  56..  64..  73.5.  87.,  99.5,  99.5, 

i  5.1,  5.9,  6.75,  7.7,  8.9,  10.4,  10.4/ 

C 

DATA  TKCRU/  19.4.  20.6,  23.4,  26.0,  28.4,  30.7, 

i  20.2,  22.5,  24.9,  27.1,  29.1,  31.1, 

i  18.1,  20.1,  22.1,  23.9,  25.7,  27.3, 

i  16.5,  18.3,  19.9,  21.4,  23.9,  24.3, 

4  20.4,  21.7,  23.1,  24.3,  25.5.  26.7, 

t  19.3,  20.5,  21.6,  22.6,  23.6,  24.5, 

4  19.3,  19.3,  20.2,  21.1,  21.9,  22.6/ 

C 

DATA  TKH/  19.6,  19.9,  26.3,  20.8,  21.6,  21.8,  22.5, 

4  23.3,  24.0,  24.5,  25.3,  25.8,  26.4,  27.1/ 

C 

C 

IF  (lOPT  .CE.  I  .OR.  lOPT  .LE.  5)  CO  TO  1495 
CALL  ERROR  <“IOPT  ERROR","  IN  KC10A") 

RETURN 

C 

C 

1495  GOTO  (1500,1500,1500,1700,18001,  lOPT 

C 

C 

C  KC10A  CLIHB  COMPUTATIONS 
C 

C  VUCNT  V2=1NIT  ALT  V3=TERM  ALT 
C 

C  J  =  1  OUTPUT  V  =  CLIMB  TIME  (MINI 
C 

C  J  -  2  OUTPUT  V  -  CLIMB  RANGE  1NM> 

C 

C  J  =  3  OUTPUT  V  =  CLIMB  FUEL  (1000  LES) 

C 

1500  J  =  lOPT 

GNT  =  AMU  (VI  ,248.1 

PAS  =  n 

PAT  =  V3 

IF  (CHT.LT.248.  .OR.  CUT. GT. 591.  .OR. 

4  PAT.LT.PRSt  CO  TO  mil 
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DO  1550  I  :  7 

K  =  I 

IFlPftT  .GE.  (43.-1*2.11  CO  TO  Uif 
155(1  COMTINUE 
l&0(i  START  :  0. 

IF  (PAS  .LT.  20.1  GO  TO  U5f 
START  =  225(K.J1 
1650  V  =  TKCLHlKiJt  -  START 

IF  ((V.LT.3.  .OR.  V.CT.23.  l.AND.  J.EQ.l)  CO  TO  22222 

IF  (IV.LT.30.  .OR.  V.CT.155.  l.AND.  J.EQ.21  GO  TO  22222 

IF  (IV.LT.1.5  .Oft.  V.CT.14.  l.AND.  J.ECl.31  CO  TO  22222 

RETURN 
C 

C  KC10A  CRUISE  COflPUTATICNS 
C 

C  Vl-CNT  V2=ni[  V3:ALT 
C 

C  OUTPUT  V  =  1000  /  HIN 
C 

1704  CUT  :  AHAIKVI. 240.1 
PA  »  V3 

IF  (GUT.LT.248.  .OR.  CUT.GT.591.  .OR. 

(  PA  .LT.  26.  .OR.  PA  .GT.  45.1  CO  TO  11111 
C 

DO  1720  I  =  h  7 
J  =  I 

•  IF  (PA  .CE.  (43. -1*211  CO  TO  1740 
1720  CONTINUE 
♦  1740  DO  1760  I  Ml  6 
K  =  I 

IF  (CUT.CE.(360*{4e*Jl-(20*Ull  GO  TO  1780 
1760  CONTINUE 

1780  V  -  TftSP(.82iPftl/TKCRU(KiJl 

IF  (V  .LT.  .2  .OR.  V  .CT.  .51  CO  TO  22222 
RETURN 
C 

C  KC10A  RENDEZVOUS  HOLDING  FUEL 
C 

C  VlsCUT  V2Mn  V3Mn 


OUTPUT  V  --  1000  LBS  /  WIN 


1800  CUT  :  AnAI[l(Vli240.1 

IF  (CUT  .LT.  248.  .OR.  CUT  .CT.  5911  CO  TO  11111 
BO  1640  I  r  1,  14 


IF  (C«T  .GE.  GO  TO 

18fl»  COMTINUE 
18fc«  V  =  k.lll  TKHtJ! 

IF  (V.LT.  .Zt  .OR.  V.GT.  .21)  CO  TO  IZZZZ 
RETURN 
C 

llin  PRINTi,“  lOPT  ^MOPTi"  CUT  =“,V1,'‘PA  --“.V3 
CALL  ERROR  (“  mPUT"i“  KCISA") 

RETURN 

22222  PRINT* f"  lOPT  CWT  PA  =‘'»V3»"  V  =‘‘»V 

CALL  ERROR  C  OUTPUT","  KCieA'M 
RETURN 
END 
C 
C 
C 

Ctft«titt 

c 

SUBROUTINE  KCi35A  (Vl,V2,V3,10PT,V) 

C 

C 

DIMENSION  A2C8A(34),A2C9A(34),A2C18A(3t) 

DIMENSION  A2F3(3.2»,  A2F2(4,4,2) 

DIMENSION  F2FI (18),  A2H1(4,7,2),  A2F1(7,3,18) 

DIMENSION  A2FlA(7,3I,A2FlBI7,3),A2Fli:{7,3!,A2FlD(7,3), 
W2FlE(7,3>,A2FlF(7,3I,A2FlC17,3),A2r'lH17,3),A2FlI(7,3l, 
M2FIJ{7,3l,A2Ftf((7,31,A2FlL(7,3l,A2FlM(7,3l,A2FIN<7,3)i 
4A2FI0('/,3),A2FIP(7,3),A2F10I7,3) 

EQUIVALENCE  (R2FtA,A2Fl(l,t,t)l,(A2FlB,A2Fl(l,l,2)), 
i  (A2FIC,A2F1(1,1,3)),(A2FID,A2F1(1,1,4)), 

t  {A2FlE,A2Fl(l,li5>),(A2FlF,A2Fl(l,l,61l, 

t  IA2FlG,A2Fl(I,l,7)),(A2FIH,A2FIIl,h8)), 

t  (A2FlI,A2Fl(til,9))>(A2FlJ,A2Flll,l>I0l) 

C 

C  CONSTANT  ALTITUDE  CRUISE 
C 

DATA  A2F1A/  .7039SE+02,  .15779E+«3,-,34584E+03,-.22t53E+e3 
i  .53343E+g3,-.15145£*03,  0.00, 

(  -.12288E+»l,  .2tl82E+81,-.339fi8E-e0,-.li085E+Sl 

(  .l28G!5E+81,-.87724£-88,  0.00, 

t  .39967E-03,  .3e30iE-02f-. 12897E-0; ,  .75422E-02 

&  .9l7ei?E-f2,-.73220E-02,  0.00/ 

DATA  A2F1B/  .25l27E+03,-.52458E4e3,  .29179E+02,  ,l0933E+e4 
t  -.n432£+04,  .27593E+03,  0.00, 

t  -.2WIIE+0I,  .712£tE+0t.-.194l2E+01,-.74513E+01 

t  .39220E401,  .140]7E*01,  0.00, 

t  .38753E-02,-.75104E-02,-.80590E-03,-.18701E-01 

t  .8l28:3E-0l,-.39947E-Bl,  0.00/ 
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(  .29337£+04.-.8S3I2E4e3t  tJl, 

t  -.U2WE400.-.«544eE+01.  .17073E+e2.-.2E518£+02i 

(  .23'.lSE+02i-.9t9t5E+01t  «.00i 

(  -.29285£-02i  .12494E-01.  .118UE-01.-.18285E409t 

t  .l3fc71E-8i.-.5545£E-0t.  0.00/ 

DATA  A2F1D/  .37544E+03.-.11455E+04.  .15523E+04.-.95008E+03. 
i  .70t24E403>'.58030E^03>  0.00> 

i  -.31465E+01t  .e0893Ef01.  .34599E-00.-.180e5E+02. 

(  .15781E*02i-.23145E40h  0.00. 

i  -.17231E-02.  .1811iE-0l.-.24390E-01.-.t0848E-01. 

(  .i4904E-00.-.83233E-01.  0.00/ 

DATA  A2F1E/  .2t77SE+83.-.10178E404.  .21887E404.-.17785E+04. 
i  -.39574E+03.  .l4393£*04.-.7i040E+03. 

t  -.14745E+01.  .38478E+01.-.17543E+01.-.45214E+01t 

t  -.417fc7E+01.  .2044tE+02.-.12100£+02. 

&  -.40467E-02.  .ie524£-ei.  .2127lE-@l.-.12339E+00. 

i  .1324aE-00.-.291UE-0l.-.15790E-01/ 

DATA  A2F1F/  .38353E+02.  .82<.>34E+03.-.25352E+04.  .38980E+04. 

t  -.15S43Emi  .88489E402)  0.00. 

i  -.4694IE401,  .18500E+02.-.30075E+02.  .14853E+02. 

(  .7835tE+8l.-.824UE+01.  0.00. 

i  -.t2904E-0l.  .794l5E*01.-.12fc53E-00.-.t0159E-82. 

i  .1492:3E-0I.-.85378E-01.  0.00/ 

DATA  A2F1C/-.U850E+03.  .12985E404,-.24248E+04.  .34314E404. 

(  -.31385E404f  .H425E404.  0.00. 

t  .95392E-00.-.95819E481.  .17752E402.-.17308E402. 

t  .17847E402,-.95525E401,  0.00. 

t  .43332E-03.  .29842E-0l.-.i4fc08E-8l.-.19098E-01. 

t  .U90IIE400.-.47919E-01.  0.00/ 

DATA  A2F1H/  ,2t270E403.-.84408E403.  .18182E404.-.12944E404. 
t  -.33707E403.  .49413E403.-.12880E403. 

(  -.15025E481,  .t798SE401,-.|4204E401.  .29452E401, 

t  -,13479E4g2.  .2e024E402,-.U49tE402, 

8  -.82385E-02.  .t09t5E-01.-.97733E-01,-.19848E-02. 

(  .84582E-0t.-.22308E-0l.-.17244E-0l/ 

DATA  A2FII/  .34995E403.-.3t2t0E403.-.55631E403.  .57083E403. 
t  .85440E403.  .494478402,-. 114948404, 

t  -.72100E401,  .11838E402,  .791448401, -.17944E402, 

(  .44275E-00, -.743288401,  .144918+02, 

(  .308348-01, -.32439E-01,-.44045E-01,  .10251E-01, 

i  .8032.18-01,  .129528-00,-.  144098-00/ 

DATA  A2F1J/  .278738. 03, -.418418+03,  .383298+03, -.113098+04, 
t  .301178+04, -.222128+04,  0.00, 

t  .573558-00,-. 709878+01,  .119378+02,-. 107878+02, 

(  .347248+01, -.289718+01,  0.00, 

t  -.874148-03,  .374788-01,-. 103448+00,  .949318-01, 

t  -.214238-01, -.103048-01,  0.00/ 


c 

c 

c 

c 


c 

c 

c 

c 


DATA  F2F1/23.0.  25. 0i  27. Bi  TC.8.  32.B. 

(  33,01  35.01  38.0>  40.01  43.0/ 

DATA  A2F2/  2.5322E+0li-7.1468E+8{.  1 .4859E+02.-8.8t23E401i 
(  -t.t354E+00.  4.7214E+30.-l.9222E+e0,-t.9480E+0tft 

(  5.9305E-62.-e.0547E-03.-2.256i£-01.  2.45e6E-01. 

t  -7.1gllE-04.  2.2244E-04,  3.436i£-03.-3.92l2E-03i 

i  -2.29t4E40t,  l.0035E+01t  4.2423E+0!,-5.7393E+00f 

i  1.4534E400»  2.7132£+B0»-3,19iEE4C0,-4.7t27E+e3. 

t  -3.94t8E-02.  5.8748E-03i-l.0305E-0l,  2.9979E-01, 

i  2.4324E-04f  5.U79E-04»  1.1279E-03»-3.4574E-03/ 

DATA  A2F3/  8.5849E-02.  8.7513E-01.  t.3284E-05, 

t  3.97g2E-02.  9.3499E-ei,  4.4710E-g5/ 


BEST  ENDURANCE  (CONSTANT  ALT  I TUBE  I 


DATA  A2H1/ 


l.fc779E-01f 

-3.2913E-02I 

3.4754E-03. 

-4.1t33E-05. 

•8.55l!»£-64. 

3.6958E-g7. 

-4.0l34E-09t 

1.2338E-02. 

5.U07E-02, 

-8.9942E-04. 

-l.5tME-05i 

-7.U94E-B8I 

9.t835E-l0. 

(>.35I1E-Uf 


4.W41E-g4, 

5.737tE-04. 

■5.4775E-05. 

1.4629E-etr 

4.9i91E-08. 

■2.3980E-09. 

1.5045E-11I 

•t.37eAE-03. 

-2.7S77E-04. 

2.5439E-04t 

9.3940E-B8. 

7.3507E-09. 

•2.17ttE-12. 

•9.3425E-13. 


3.2904E-0t 

-2.5602E-0£ 

2.2130E-07 

5.2e48E-09 

•1.2tt5E-10 

2.3e22E-12 

9.9904E-14 

4.5348E*05 

■9.2074E-07 

2.838ZE-09 

8.82i-7E-l0 

2.097BE-11 

5.4770E-13 

•3.51ME-15 


-6.9t89E-09, 

3.3941E-09. 

-2.4929E-1S, 

3.2877E-I2t 

1.1304£-t3f 

9.4055£-15i 

-4.2445E-14I 

2.4750E-09, 

-7.1474E-09, 

2.5989E-19. 

1.3121E-12. 

-1.4532E-13. 

-l.39fc3E-15. 

9.263fcE-l7/ 


280  KCAS  CLIMB  DATA  (MRTiTEMP  DEV  =  01 
8  =  TlMEi  9  =  RANGEi  10  =  FUEL  MILAGE 


DATA  A2i:8A/ 

t.l884279ZE+02.  .4I3638I7E+0I.-.55344857E400I  .192612lgE-B!» 
l-.29292:379E-B3..1582932lE-05.-.(>3323548E+00i-.89573058E-01i 
l.59e25186E-02i-.13521327E-03.  .105t4593E-B5.  .0300000eE+00» 
4.77333210E-a2i  .47702502E-03.-.1989t540E-04,  ,2S8tl935E-04. 
t . 0000001I0E+00 .  . giWBggBBE+eO. - . 4375535 1 E- 04 .  - .  1 37 13 1 72E-05 1 
t.22594734E-07.  .03000e00E+03.  .006g0t’90E+00.  .03000000E+03. 
(.1U43423E-04I  .&3801229E-09.  .00000030E400.  .03000030Ei03> 
i.mmmmf  .0(f000030E403>~.11782181E-09)  .03000000E403, 
(.00000030E400r  ,03000000E400|  .00000030E400I  .00000000E400/ 
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DATA  AZCi’A/ 

fc.3M85743E+«2.  .325t8397E+02i-.370t7i32E+ei.  .U247979£+e«. 
fc- .23552449E-82. . 13631375E-04 i - . 1 44534748+81 . - .55894354E+03 . 
8.34877233E-0l»-.93074950£-03i  .74295143E-t'5i  .83gi)6g98E+03r 
t.28797579E-6U  .3n99l42E-02t-.il842322E-e3t  .13343148E-05. 
8.«08e00«8E+e8.  .e!J080008E+0«i-. 132533948-03. -.772408878-05. 
6.13443422E-04.  .0.9000000E+00.  .000000IJ0E+00.  .07000000E+0l>. 
t.38847141E-04.  .71828750E-0O,  .00000000E+f0.  .03000800E+efl. 
i .  800000(10£+00.  .  i(lB00000E+03 .  - .  42401 239E-  09 .  .  04000000E+0tf . 
&.00000000£+00>  .000000008+00.  .000000008+00.  .000000008+00/ 
DATA  A2i:i0A/ 

t. 947790488+00.  .104107078+01. -.940145458-01.  .313007938-02. 
t- .444419428-04. .223728448-04.- .398774278-01 . - . 154799228-01 . 
A. 104442448-02. -.234049258-04.  .147807438-04.  .000000008+03. 
A.543137488-03.  .837474238-04. -.374187488-05.  .433303388-07. 
4.000000308+00.  .03000e00£+63.-.32533517E-C5.-. 209923128-04. 
4.438177318-08.  .030000008+03.  .000000338+00.  .030000008+03. 
4.893448458-08.  .198331838-09.  .000003308+00.  .030000008+03. 
4.006000308+00.  .03300000E+03.-.92200077E-11.  .030000C0E+63. 
4.000000308+00.  .033000008+03.  .000300338+00.  .030030008+03/ 
C 

IF  (lOPT.CE.l  .OR.  I0PT.LE.5I  CO  TO  1950 
CALL  ERROR  ("lOPT  ERROR"."  IN  KC135A"I 
RETURN 

1950  CO  TO  (2030.2100.2200.2300.2400).  lOPT 
C 

c 

c 

C  CONSTANT  ALTITUDE  CRUISE 

C 

C  V1=CMT  V2=HACH  V3=ALT 
C 

C  OUTPUT  V  i  1000LES  /  HIN 

C 

2300  CU  =  AKAX1(V1.100.) 

AHACH  --  V2 
PA  =  V3 

2310  IF  (PA  .LT.  0,0  .OR.  PA  .CT.  50.  .OR. 

4  ANACH  .LT.  0.2  .OR.  ARACH  .CT.  0.9  .OR. 

4  CU.LT.  100.  .OR.  CU.CT.  300.)  CO  TO  11111 
DO  2320  J  =  2.  13 
I  =  J-l 

IF  (PA  .LT.  r2Fl(.J)l  CO  TO  2330 
2320  CONTINUE 

2330  Tl  --  P0LLT(ALFl(l.l.l).7.3.ARACH.CNt 
T2  =  P0LLT(A2Fl(l.I.I+ll.7.3.A«ACH,CHI 
T  s  Tl  +  (PA-F2FllIt)  ♦  (T2-Tl)/(F2Fl(I+l)-F2FllI)) 
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I 


IF  IT  .LT.  P0LLT(fliF2(l.l.lti4.4iAMACHiPft!)  GJ  ; .  '.m 
T  =  P0LT(ft2F3(l.n.3iT) 

CO  TO  235« 

234(1  IF  (T.LT.P0LLT(A2F2(lil.2).4.4tA«ACH.Pftl)  CO  TO  2C5B 
T  =  P0LT(A2F3(lf2).3,Tl 

23511  IF  n  .LT.  14.  .OR.  T  .CT.  80.  1  CO  TO  22222 

V  =  TASPIAMACHtPA)  I  T 
RETURN 

C 

C 

C  BEST  ENDURANCE  (CONSTANT  ALTITUDE) 

C 

C  V1=CHT  V2--HACH  V3=ALT 

C 

C  OUTPUT  V=  IBSBLBS  /  HIN 

C 

2400  CU  =  AnAl(l(Vi>100.) 

PA  :  V3 
I  =  I 

IF  (PA  .CT.  30.0  )  I  =  2 

AHACH  =  P0LLr(A2Hl(l.l.I),4i7.GW.PA) 

CO  TO  2310 
C 

c 

C  ZeOKCAS  CLINE  TINE.RANCE.FUEL  (NRT) 

C 

C  Vl=CWT  V2--INIT  ALT  V3=TERH  ALT 

C 

C  OUTPUT  V  =  TINE  TO  CLINE  (NIN) 

C 

2000  CU:  ANAIl(Vltl00.) 

PAS  =  V2 
PAT  s  V3 

T1  =  P0LLT(A2C8A.4.ttPAS.CH) 

T2  =  P0LLT(A2C8Ai4i6iPATiCHI 

V  =  TZ  -  Tt 

IF  (V  .LT.  0.0  .OR.  V  .CT.  40.  I  CO  TO  22Z22 
CO  TO  2250 
C 

C  OUTPUT  V  =  RANGE  IN  CLINE  (NH) 

C 

ZlfO  CU  =  ANAinVl,!’?.) 

PAS  =  V2 
PAT  =  W 

T1  =  P0LLT(A2C'A>o.4.PAS.CH) 

T2  =  P0LLT(A2C?A.4.4iPATfCNI 

V  :  TZ  -  T1 


c 

C  OUTPUT  V  =  FUEL  TO  CLIMB  U«0«LES) 

C 

2ZM  CU  :  AMAH (VI 1 100.) 

PAS  =  V2 
PAT  =  V3 

Tl  :  POLLT(A2C10A.ii6.PAS.CU) 

T2  =  P0LLT(A2ClBA.S.SiPATtCW) 

V  =  T2  -  Tl 

IF  (V  .LT.  0.0  .OR.  V  .CT.  17.  )  GO  TO  22222 
2250  IF  (CM  .LT.  100.  .OR.  GM  .GT.  300.  .OR. 

i  PAS  .LT.  0.0  .OR.  PAT.GT.  50.  )  GO  TO  11111 

RETURN 
C 
C 
C 

mil  CALL  ERROR  (“INPUT  ERRO"i"R  TO  KC135“) 

RETURN 

C 

22222  CALL  ERROR  ("OUTPUT  ERR" t "OR.  KC135A") 

RETURN 


This  appendix  describes  the  procedure  used  to 
validate  the  SLAM  simulation  model.  Validation  of  the 
SLAM  model  is  broken  down  into  two  steps:  (1)  verifica¬ 
tion,  and  (2)  external  validation.  Verification  insures 
that  the  model  behaves  as  intended  and  designed;  this 
includes  internal  mechanical  verification.  External  vali¬ 
dation  tests  the  agreement  between  the  model  output  data 
and  the  true  output  obtained  from  the  real  world.  The 
validity  of  the  SLAM  model  was  only  evaluated  in  terms  of 
the  model's  purpose  which  is  to  validate  the  analytical 
model's  output  under  the  real  world  conditions  of  stochas¬ 
tic  wind  factors  and  ground  delays. 

Verification  of  the  SLAM  model  ensured  that  the 
model  produced  mechanically  correct  results  and  that  the 
results  were  calculated  as  intended.  This  verification 
step  was  done  independently  three  times  during  the  crea¬ 
tion  of  the  model  due  to  discovery  of  internal  errors. 
During  the  verification  process  errors  in  design  and  pro¬ 
gramming  were  found  and  corrected.  After  each  design  or 
programming  change,  the  verification  process  was  repeated 
until  the  SLAM  model  output  was  verified  to  be  correct 
and  -the  internal  mechanical  aspects  of  the  model  were  also 
determined  to  be  correct. 

The  verification  process  consisted  of  two  stages; 

1.  Verification  of  the  random  wind  variables  and 
probabilistic  ground  delays. 
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2.  Internal  mechanical  verif  icati.:>r. . 

The  first  stage  verified  that  the  climb  and  cruise 
winds  as  specified  in  the  SLAM  model  were  actually  pro¬ 
viding  sample  winds  for  those  distributions  and  that 
delays  occurred  appropriately.  For  example,  the  cruise 
wind  used  a  normal  distribution  (u  =  .917,  o  =  .16  NM/MIN) . 
To  determine  that  this  SLAM  wind  distribution  actually 
provided  the  model  with  a  normal  distribution  as  defined, 
wind  values  were  obtained  from  the  model  and  tested  using 
a  Kolmogorov-Smirnov  Goodness-of-Fit  Test.  Figure  I-l  con¬ 
tains  the  Vogelbach  Computing  Center  Statistical  Package 
for  the  Social  Sciences  (SPSS)  Computer  Program  used  to 
test  the  wind  output  distribution.  The  Kolmogorov-Smirnov 
Test  tested  the  SLAM  wind  distribution  against  a  theoreti¬ 
cal  normal  distribution  with  the  p  and  a  previously  defined 
The  null  hypothesis  (Hq)  used  in  the  test  was  that  "There 
is  no  significant  difference  between  the  observed  data  and 
that  given  by  normal  distribution  with  mean  equal  to  .917 
and  a  standard  deviation  equal  to  .16.  For  each  case  the 
Kolmogorov-Smirnov  critical  value  was  greater  than  the 
Kolmogorov-Smirnov  value  calculated  by  the  SPSS  Program; 
therefore,  Hq  was  not  rejected.  However,  since  the  mean 
and  standard  deviation  are  estimated  parameters  the  Lillie- 
fors  Test  Statistic  should  also  be  used  to  test  Hq .  The 
example  in  Figure  I-l  contains  ten  samples.  The  Kolmogorov 
Smirnov  (K-S)  critical  value  for  ten  samples  and  a  =  .05 


153 


VOGELBACK  COMPUTING  CENTER 
NORTHWESTERN  UNIVERSITY 

SPSS  —  STATISTICAL  PACKAGE  FOR  THE  SOCIAL  SCIENCES 
VERSION  8.0  —  JUNE  18,  1979 


RUN  NAME 
VARIABLE  LIST 
N  OF  CASES 
INPUT  FORMAT 
NPAR  TESTS 
READ  INPUT  DATA 


K-S  SLAM  WIND  VARIANCE  VALIDATION 
X 
10 

FREEFIELD 

K-S (NORMAL, .9, .16)=X 


GIVEN  1  VARIABLES,  INITIAL  CM  ALLOWS  FOR  1736  CASES 

MAXIMUM  CM  ALLOWS  FOR  11336  CASES 


IK-S  SLAM  WIND  VARIANCE  VALIDATION 

FILE  NONAME  (CREATION  DATE  =  02/26/81) 

—  KOLMOGOROV  -  SMIRNOV  GOODNESS  OF  FIT  TEST 

TEST.  DIST.  -  NORMAL  (MEAN  =  .9000  STD.  DEV.  =  .1600) 

CASES  MAX(ABS  DIFF)  MAX(+  DIFF)  MAX(-  DIFF) 

10  .2368  .1599  -.2368 


K-S  Z  2-TAILED  P 

.749  .629 


Fig.  I-l.  SLAM  Wind  Variance  Validation;  SPSS  Program 
and  Results;  Kolraogorov-Smirnov  Test 
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(Ref.  Shannon,  Systems  Simulation ;  The  and  Science , 

Prentice-Hall,  1975)  is  0.410.  The  critical  value  for  ten 
samples,  a  =  .05  from  Lilliefors  (1967)  is  0.258.  Since 
our  largest  deviation  was  .2368,  we  do  net  reject  the  Hq 
and  the  wind  distribution  is  validated  under  with  the  K-S 
and  Lilliefors  criteria. 

The  second  stage  of  verification  consisted  of 
ensuring  the  internal  mechanical  operations  of  the  SLAM 
model  were  functioning  as  designed.  To  accomplish  this  the 
SLAM  Output  Trace  (see  Appendix  G)  was  completely  analyzed 
for  mechanical  verification.  The  items  verified  were  the 
time,  distance,  and  the  fuel  calculations  in  the  user  func¬ 
tions  and  the  network  portion  of  the  model.  The  mechanical 
verification  of  the  ASSIGN  nodes  was  accomplished  for  each 
node.  For  example,  in  ASB  (see  Appendix  K  for  SLAM  descrip¬ 
tion)  (Yokota  Scenario — C-141/KC-135)  the  C-141  cruise  to 
rendezvous  point  fuel  consumption  is  computed  in  USERF  13 
as  19,300  lbs.  and  this  fuel  is  added  to  ATRIB(4),  airlifter 
total  fuel  consumption.  Looking  at  a  SLAM  Trace  (see 
Appendix  G) ,  ATRIB(4) 's  value  prior  to  the  completion  of 
node  ASB  was  11,280  lbs  which  represented  fuel  consumption 
up  to  that  point.  Upon  departing  node  ASB,  ATRIB(4) 
equalled  30,  650  lbs  (11,280  +  1,936)  which  was  determined 
to  be  correct.  Throughout  the  SLAM  model  each  node  was 
verified  in  this  way  to  insure  the  results  were  mathe¬ 
matically  and  functionally  correct. 


^ -  - - -  ....aJi— Aj 
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The  second  stage  of  internal  mecli:  nical  verifica¬ 
tion  was  the  verification  that  the  SLAM  discrete  event 
simulation  was  providing  the  SLAM  network  model  with  the 
proper  input  information.  This  was  easily  verified  by 
assigning  SLAM  input  parameters,  such  as  aircraft  ramp 
fuel  and  takeoff  times,  calculated  by  the  analytical  model, 
to  global  variables  in  the  discrete  part  of  the  SLAM  model, 
node  EVENT  1.  These  input  values  are  then  made  to  print 
out  on  the  SLAM  Trace  as  they  are  used  in  the  user  func¬ 
tions.  The  agreement  between  the  input  values  assigned  to 
global  variables  and  those  used  in  the  user  functions  was 
100  percent.  Therefore,  this  verified  that  the  SLAM  net¬ 
work  model  was  indeed  using  the  proper  input  parameters 
defined  at  the  start  of  the  simulation.  By  verifying  the 
SLAM  network  input  parameters  by  this  method,  validation 
of  the  input  values  themselves  was  not  necessary;  only  the 
verification  that  they  were  transferred  correctly  to  the 
user  functions  as  required  was  necessary.  Since  these 
SLAM  input  parameters  have  been  previously  validated  in 
the  analytical  model,  it  can  be  said  with  confidence  that 
the  input  parameters  for  SLAM  are  not  only  verified  but 
also  validated. 

By  insuring  that  the  internal  mechanical  function 
of  the  model  is  correct,  a  significant  degree  of  confidence 
is  built  up  in  the  model.  If  the  model  could  then  be  shown 
to  predict  flight  plan  fuel  consumptions  correctly,  then 
validation  of  the  SLAM  model  is  completed. 
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To  validate  the  SLAM  model,  the  ir.iut  parameters 
computed  for  the  analytical  model,  which  have  already  been 
validated  against  real  world  data,  were  used  to  manually 
compute  total  fuel  consumed  for  each  aircraft.  This  was 
done  for  each  aircraft  combination  and  each  scenario.  The 
fuel  consumption  figures  using  SLAM  compared  to  manual  com¬ 
putations  are  shown  in  Table  4.  As  can  be  seen,  the  total 
aircraft  fuel  consumption  SLAM  figures  are  within  7  percent 
of  the  manual  computations.  These  reassuring  results 
increase  the  validity  of  the  SLAM  model. 

The  verification  and  validation  resulted  in 
increasing  confidence  in  the  SLAM  model.  When  the  analyti¬ 
cal  model  verification  and  validation  processes  are 
included  the  models  have  been  validated  to  the  point  where 
experimental  data  can  be  obtained  and  used  with  confidence 
over  a  varying  range  of  input  parameters. 

The  final  validation  of  FLTPLN  was  accomplished 
using  the  SLAM  model.  The  FLTPLN  output  (fuel  loads, 
distances,  takeoff  times,  etc.)  was  entered  as  SLAM  input 
parameters  and  SLAM  output  (total  fuel  consumption  per 
aircraft)  .  Figures  3-6  show  the  percent  deviation 
of  the  SLAM  fuel  consumption  mean  compared  to  the  analytical 
model  fuel  consumption  value.  A  negative  value  means  that 
the  SLAM  model  used  less  fuel  than  the  analytical  model  and 
a  positive  value  means  the  SLAM  model  consumed  more  fuel. 

The  specific  SLAM  results  are  discussed  in  con¬ 
junction  with  Tables  L-1  to  L-4  by  aircraft.  Generally, 


the  SLAM  results  for  all  aircraft  and  cajcjo  loads  compare 
favorably  with  the  analytical  results,  which  indicates 
that  the  analytical  data  is  realistic  and  can  be  used  with 
confidence.  Overall,  the  tankers  burned  about  1  percent 
less  fuel  in  the  SLAM  model  than  in  the  analytical  model . 
Conversely,  the  airlifters  burned  about  2.5  percent  more 
fuel  in  the  SLAM  model  than  in  the  analytical  model.  This 
difference  in  airlift  fuel  can  be  explained  from  the  nature 
of  the  two  aircraft  missions.  The  tanker  flies  to  a 
rendezvous  point  and  returns  to  his  departure  base  which 
is  relatively  close  compared  to  the  airlifter  destination; 
the  airlifter,  however  flies  to  the  rendezvous  points  and 
then  continues  on  to  his  destination  which  is  thousands 
of  miles  farther  than  the  tanker  route  of  flight.  Conse¬ 
quently,  the  airlifter  has  more  time  to  be  affected  by 
increases  in  the  wind  factor,  and  these  increases  plus  any 
delay  fuel  consumption  would  cause  the  increased  fuel 
requirement.  There  are  two  reasons  the  tanker  SLAM  results 
are  less  than  the  analytical  consumption.  The  first  is 
that  an  increase  in  wind  factor  has  less  effect  on  the 
tanker  than  on  the  airlifter.  The  headwind  increase  to 
the  rendezvous  point  becomes  a  tailwind  increase  when 
returning  to  destination.  The  second  is  the  difference  in 
the  loiter  fuel  computations.  In  FLTPLN  the  tanker  loiter 
fuel  is  always  based  on  a  fifteen-minute  3,oiter  duration; 
however,  in  SLAM,  the  loiter  time  is  based  on  actual 


aircraft  arrival  s  which  usually  requires  less  than  fifteen 
minutes  loiter  time.  For  the  tankers,  fifteen  minutes 
loiter  fuel  is  more  than  the  difference  between  SLAM  and 
FLTPLN  fuel  values. 

The  final  validation  of  the  model  is  subjective 
and  consisted  of  asking  the  question,  "Does  the  output 
data  make  sense?"  Speaking  as  a  MAC  pilot  and  navigator, 
the  output  fuel  consumptions  do  make  sense.  As  distance 
between  the  departure  and  destination  bases  increases, 
the  total  fuel  consumption  also  increases.  In  the  case  of 
the  tanker,  as  the  distance  from  the  departure  base  to  the 
rendezvous  entry  point  increases,  total  tanker  fuel  con¬ 
sumption  also  increases.  As  the  cargo  decreases  on  the  air- 
lifter,  the  total  fuel  consumption  decreases  for  each  air- 
lifter  in  each  instance.  Also,  for  a  route  segment  with  a 
tail  wind,  fuel  consumption  decreases  compared  to  a  com¬ 
parable  route  segment  with  a  head  wind.  All  of  these 
general  ooservations  from  the  data  make  sense  to  an  opera¬ 
tional  pilot  or  navigator  and  none  of  them  is  counter¬ 
intuitive.  The  result  of  these  general  observations  by 
professional  aircrew  is  that  the  model  validity  is  increased 
and  will  help  reassure  a  future  decision  maker  or  user  of 
the  completeness  of  the  research  and  the  validity  of  the 
results.  By  obtaining  reasonable  results,  user  concerns 
about  the  details  of  the  model  or  the  methodology  are 


minimized. 


This  appendix  explains  the  methodology  used  to 
calculate  the  95  percent  confidence  intervals  portrayed 
in  Tables  L-1  to  L-4,  Appendix  L,  and  provides  a  sample 
calculation. 

A.  The  results  of  the  SLAM  model  replications  for 
a  particular  mission  scenario,  aircraft  combination,  and 
cargo  load  were  combined  to  obtain  a  sample  mean.  How¬ 
ever,  since  they  are  only  estimates,  a  confidence  interval 
surrounding  the  mean  is  a  means  of  indicating  the  accuracy 
of  the  results  compared  to  the  true  population  mean.  Given 
the  sample  mean  (x) ,  the  sample  size  (n) ,  and  the  sample 
standard  deviation  (s) ,  a  confidence  interval  can  be 
created  using  the  t  statistic.  From  the  confidence  inter¬ 
val,  with  an  upper  and  lower  boundary,  it  can  be  said  that 
the  probability  of  the  mean  being  between  the  upper  and 
lower  limits  is  equal  to  probability:  1  -  a.  Alpha  (a)  is 
defined  as  .05,  hence  a  95  percent  confidence  interval 
will  be  created.  If  normality  is  assumed  (see  Part  B  for 
normality  validation  of  fuel  consumption  data)  a  confidence 
interval  can  be  estimated  using: 

P  (LL  <  UL)  =  1  -  a  =  1  -  .05  =  0.95 


where, 

LL  =  X  -  *  s  //  n-1 

UL  =  X  -  t^^.^  *  s/  J  n-1 


=  critical  values  of  t  for  n-l  degrees  of 
freedom;  and 

a  =  .05. 

and  where, 

LL  =  lower  limit; 

UL  =  upper  limit; 

X  =  sample  mean; 

s  =  sample  standard  deviation;  and 
n  =  sample  size. 

From  Hines  and  Montgomery  (1972) : 

tcrit  “  2.26  for  a  sample  size  of  10  with  9 
degrees  of  freedom. 

Therefore,  the  confidence  interval  for  a  KC-135  departing 
McGuire,  refueling  a  C-141B  with  55,000  pounds  of  cargo 
with  10  samples  is: 

SLAM  Mean  (x)  =  90,840  lbs. 

Sample  Standard  Deviation  (s)  =  1,690  lbs. 

=  (2.26)  (1.69)  11000)  „  „ 

X  +  *  s//  n-1  =  90,840  +  1270  =  92,110  lbs 

X  -  *■  s/  /  n-1  =  89,570 

Therefore,  the  confidence  interval  is: 

P{89570  <  M  <  92,110)  =  .95 
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B.  Normality  Verification.  To  compute  a  confi¬ 
dence  interval  using  the  sample  t  statistic  as  in  Part  A 


of  this  appendix,  normality  must  be  assumed.  Using  the  K-S 
Test,  using  Lilliefors'  Statistic,  normality  will  be  vali¬ 
dated  for  the  SLAM  output.  Figure  J-1  is  the  SPSS  Program 
and  results  for  a  SLAM  run  (C-141B  /  KC-135 ,  tanker  base 
is  McGuire  AFB,  cargo  load  is  55,000  lbs).  For  the  KC-135 
the  largest  deviation  was  .2368;  since  this  is  not  larger 
than  the  Lilliefors  critical  value  of  .258  (n  =  10,  a  =  .05) 
the  Hq  that,  "There  is  no  significant  difference  between 
the  data  and  a  normal  distribution,"  cannot  be  rejected. 
Therefore,  normality  can  be  assumed  for  the  SLAM  output 
and  the  .t  statistic  can  be  used  to  create  a  95  percent  con¬ 
fidence  interval. 


VOGELBACK  COMPUTING  CENTER 
NORTHWESTERN  UNIVERSITY 

SPSS  —  STATISTICAL  PACKAGE  FOR  THE  SOCIAL  SCIENCES 
VERSION  8.0  —  JUNE  18,  1979 


RUN  NAME 
VARIABLE  LIST 
N  OF  CASES 
INPUT  FORMAT 
NPAR  TESTS 
READ  INPUT  DATA 


K-S  SLAM  KC-135  NORMALITY  VALIDATION 
X 
10 

FREEFIELD 
K-S (NORMAL) =X 


GIVEN  1  VARIABLES,  INITIAL  CM  ALLOVJS  FOR  1736  CASES 

MAXIMUM  CM  ALLOWS  FOR  11336  CASES 


IK-S  SLAM  KC-135  NORMALITY  VALIDATION 
FILE  NONAME  (CREATION  DATE  =  02/26/81) 

—  KOLMOGOROV  -  SMIRNOV  GOODNESS  OF  FIT  TEST 


TEST  DIST.  -  NORMAL  (MEAN  =  90.8400  STD.  DEV.  =  1.6959) 


MAX(-  DIFF) 
-.2294 

K-S  Z  2-TAILED  P 

.725  .669 


CASES  MAX(ABS  DIFF) 

10  .2294 


MAX(+  DIFF) 
.1402 


Fig.  J-1.  SLAM  Output  Normality  Validation; 
SPSS  Program  and  Results; 
Kolmogorov-Smirnov  Test 
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The  experimental  model  consists  of  two  models, 
an  analytical  model  and  a  SLAM  model.  This  appendix 
provides  a  detailed  description  of  both  models. 

"FLTPLN"  is  a  FORTRAN  model  which  computed  the 
minimum  fuel  required  to  operate  an  airlifter  and  tanker 
aircraft  through  the  optimal  rendezvous  point  for  a  spe¬ 
cific  scenario.  Inputs  to  the  program  are: 

1.  The  aircraft  combination. 

2.  The  airlift  cargo  weight. 

3.  The  average  wind  for  the  scenario. 

4.  The  geographic  locations  of  the  departure, 
destination,  and  abort  bases. 

The  program  outputs  include: 

1.  The  coordinates  of  the  optimal  rendezvous 

point . 

2.  The  total  fuel  and  all  route  segment  fuels. 

3.  The  takeoff  fuel  loads  for  both  aircraft. 

Additional  parameters  which  are  generated  and  used  as 
inputs  to  the  SLAM  simulation  will  be  described  later  in 
this  appendix. 

"FLTPLN"  is  the  executive  program  calling  several 
subroutines  to  accomplish  the  flight  planning  computations. 
The  FLYME  set  of  subroutines  provides  the  fuel  data  base 
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for  all  calculations.  The  FORTRAN  sourer;  r-c  de  for  FLTPLN 
is  listed  in  Appendix  C. 

The  "FLTPLN"  program  begins  by  computing  the  maxi¬ 
mum  allowable  fuel  loads.  This  forms  one  boundary  of  the 
region  of  feasible  rendezvous  points.  A  geographic  point 
on  this  boundary  is  defined  as  the  rendezvous  exit  point. 

A  second  point,  250  NM  closer  to  the  airlifter  departure 
base,  is  defined  as  the  rendezvous  entry  point.  These  two 
points  define  the  rendezvous  track.  If  the  exit  point  is 
within  400  NM  of  the  airlifter  departure  base  the  entry 
point  is  redefined  at  the  distance  flown  by  the  airlifter 
in  initial  climbout.  Rendezvous  is  then  planned  when 
25,000  feet  is  first  reached.  The  exit  point  is  redefined 
250  NM  downtrack.  These  rendezvous  points  fall  along  the 
Great  Circle  course  between  the  airlifter 's  d-^parture  and 
destination  bases.  Subroutine  RHOTHTA  returns  this  course 
and  distance  from  the  equation: 


9  =  COS 


-1  I 


-1  |SIN(LAT2)  -  SIN(LATl)  *  COS(D) 


SIN(D)  *  OOS(LATl) 


where , 


9  =  the  Great  Circle  course; 

LATl  =  the  latitude  of  the  first  point; 

LAT2  =  the  latitude  of  the  second  point;  and 

D  =  the  distance  between  these  points  computed 
from  the  equation; 
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D  =  cos  SINdATl)  *  SIN(LAT2)  +  COS(IATl)  *  COS(LAT2) 


*  COS(IOaG2-LONGl) 


where, 

LONGl  =  the  longitude  of  the  first  point;  and 
L0NG2  =  the  longitude  of  the  second  point. 


The  geographic  coordinates  of  these  points  are  computed  in 
subroutine  LATLONG  by  the  equations: 


LAT3  SIN  SIN(LATl)  *  COS(D)  +  COS(IATl)  *  SIN(D)  *  COS (6) 


WNG3  =  LONGl  +  COS 


(OOS(D)-SIN(LATl)  *  SIN(LAT2) 
00S(IAT1)C0S(LAT2) 


All  angles  used  in  these  equations  are  computed  in  radians 
with  the  following  relationships: 

1.  One  degree  of  latitude  =  50  nautical  miles  (NM) . 

2.  One  radian  =  57.2957  degrees. 


The  fuels  for  the  individual  route  segments  between  the 
coordinates  are  computed  using  the  equation: 


FUEL  =  DI STANCE /GROUNDSPEED  *  FUEL  RATE 


where  the  distance  is  obtained  from  subroutine  RHOTHTA,  the 
fuel  rate  is  obtained  from  the  subroutine  FLYME  data  base, 
and  ground  speed  is  obtained  from  the  statement  function  of 
the  wind  triangle  equation: 
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GROUND  SPEED  =  TAS  *  COS(SIN~^  (W/TAS*Sr';  CTL-OX:) ) ) 

-  WV  *  CX)S(WD*TC) 

where  TAS  is  the  aircraft  true  airspeed  for  a  given  cruise 
mach  number  at  the  cruise  altitude  and  is  computed  in  sub¬ 
routine  TASP,  WD  is  the  wind  direction,  WV  is  the  wind 
velocity,  and  TC  is  the  rumbline  course  between  the  points 
computed  in  subroutine  RHOTHTA.  The  wind  is  held  constant 
throughout  this  study  at  260°  at  55  KNOTS.  The  fuels  for 
eight  route  segments  are  computed  in  the  following  order: 

1.  Airlifter  rendezvous  exit  to  destination. 

2.  Airlifter  departure  to  rendezvous  entry. 

3.  Airlifter  rendezvous  entry  to  departure  base 
for  abort. 

4.  Airlifter  rendezvous  entry  to  alternate  base 
for  abort. 

5.  Airlifter  rendezvous  track. 

6.  Tanker  rendezvous  exit  to  recovery  base. 

7.  Tanker  rendezvous  track. 

8.  Tanker  departure  to  rendezvous  entry. 

The  fuels  are  computed  in  the  reverse  order  from  which  they 
are  flown  so  that  fuel  loads  determined  for  each  route 
segment  are  the  minimum  feasible  to  accomplish  the  remainder 
of  the  flight.  The  fuel  rate  returned  by  subroutine  FLYME 
is  dependent  on  the  cruise  altitude  and  gross  weight  of  the 
aircraft.  A  step  climb  profile  is  approximated  by  com¬ 
puting  the  fuel  rate  at  the  average  gross  weight  and 


altitude  for  the  route  segment.  These  a'7cra.ge  values  are 
obtained  by  estimating  half  of  the  fuel  to  be  used.  The 
fuel  computations  are  repeated  until  the  fuel  rate  repre¬ 
sents  the  average  for  the  route  segment.  Experimentation 
with  this  process  determined  that  the  fuels  converged  to 
within  1000  pounds  at  the  second  iteration,  50  pounds  at 
the  third,  and  2  pounds  by  the  fourth.  Three  iterations 
were  selected  for  the  averaging  process. 

Climb  profiles  are  computed  for  the  initial 
departure  segments  up  to  cruise  altitude  and  the  post¬ 
rendezvous  segments  from  25,000  feet  to  the  cruise  alti¬ 
tude.  The  climb  time,  range,  and  fuel  consumed  are 
returned  from  FLYME  and  added  to  their  corresponding  com¬ 
ponents  in  the  fuel  computations.  The  initial  takeoff 
fuel  load  for  the  airlifter  is  the  maximum  allowable  for 
the  specified  cargo  load.  The  fuel  to  be  transferred 
from  the  tanker  is  the  difference  between  the  airlifter 's 
total  fuel  requirements  and  its  takeoff  fuel.  The  tanker' 
takeoff  fuel  load  is  computed  as  its  total  fuel  require¬ 
ments  plus  the  transfer  fuel.  All  aircraft  maintain  a 
fuel  reserve  to  accomplish  holding,  flight  to  an  alternate 
and  landing  after  the  final  cruise  segment.  This  reserve 
is  30,000  lbs  for  the  C-5A  and  KC-10  and  15,000  lbs  for 
the  C-141  and  KC-135. 

After  each  fuel  computation,  the  fuel  load  is 
checked  for  feasibility  against  the  current  maximum  fuel 
capacity.  Both  aircraft  must  maintain  sufficient  fuel  to 


abort  to  the  departure  base  or,  for  the  i-irlifter,  to  the 
alternate  base,  if  the  alternate  is  closer  than  the  depar¬ 
ture  base.  If  the  required  fuel  load  is  feasible,  the  sum 
of  the  route  segment  fuels,  not  including  the  abort  fuels 
is  stored  as  a  matrix  data  point.  If  infeasible,  the 
value  777.777  is  stored  in  the  matrix. 

Next,  the  coordinates  of  the  rendezvous  exit  point 
are  adjusted  in  five  degree  increments  north  or  south, 
whichever  is  closer  to  the  direction  of  the  tanker  base 
and  the  entire  routine  is  recomputed.  Five  latitude  values 
covering  twenty  degrees  are  used.  The  latitude  is  then 
reset  at  the  initial  value  and  the  longitude  is  adjusted 
in  two  degree  increments  towards  the  airlifter  destination. 
The  five  values  of  latitude  are  repeated  for  each  of  the 
thirteen  values  of  longitude  covering  twenty-four  degrees. 
The  5  X 13  matrix  of  total  fuel  values  is  sorted  and  the 
minimum  value  from  the  sixty-five  points  is  selected  as  the 
output  variable.  The  fuel  loads  for  both  aircraft  for  this 
point  are  increased  by  the  respective  start,  taxi,  and 
takeoff  constant  and  these  values  represent  ramp  fuel 
onloads  and  are  listed  along  with  the  planned  transfer 
fuel . 

At  this  point,  the  airlifter  takeoff  fuel  load  is 
decremented  by  20,000  pounds  generating  matrices  for  ten 
values  of  the  fuel  load.  Only  matrices  which  contain  at 
least  one  feasible  value  are  printed  and  if  no  feasible 
values  are  computed,  the  statement  "NO  FEASIBLE  SOLUTION" 


is  printed.  The  minimum  values  for  each  ■r.atrix  are  sorted 
and  the  overall  minimum  fuel  is  printed  for  this  specific 
scenario,  cargo  load,  and  aircraft  combination. 

The  initial  cargo  weight  for  each  airlifter  is  the 
maximum  allowable  cargo  defined  by  operational  or  perform¬ 
ance  limitations.  Two  additional  cargo  weights  are  cor^- 
sidered.  One  weight  is  the  lowest  that  could  be  expected 
if  the  airlifter  were  filled  with  bulky,  lightweight  cargo, 
and  the  other  is  the  mean  of  the  maxim.um  weight  and  the 
minimum  expected  weight.  For  the  C-5A,  these  weights  are 
100,000,  85,000,  and  70,000  lbs.  For  the  C-141B,  the  cargo 
weights  are  70,000,  55,000,  and  40,000  lbs.  Three  separate 
runs  are  iteratively  made  at  these  cargo  weights  for  a 
specific  airlifter  and  tanker  combination.  Since  data  was 
desired  from  all  four  aircraft  combinations,  for  a  spe¬ 
cific  scenario,  the  combinations  are  changed  within  the 
FLTPLN  program.  The  combinations  are  listed  in  Table  1. 
Thus,  only  the  specific  scenario  must  be  changed  to  obtain 
data  for  all  four  aircraft  combinations  at  the  six  desired 
cargo  weights. 

The  SLAM  model  was  used  to  verify  the  analytical 
solution  with  the  addition  of  two  stochastic  variables, 
winds  and  takeoff  delays.  The  SLAM  model  also  determines 
the  flight  times  for  each  aircraft.  The  model  follows  the 
mission  profile  shown  in  Figure  K-1  and  conforms  to  the 
flight  planning  segments  discussed  in  Chapter  II. 

Figure  K-2  is  the  SLAM  network  description  of  the  model. 
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Fig.  K-1.  SLAM  Structural  Mission  Profile 


Fig.  K-2,  SLAM  Network 
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See  Appendix  D  for  a  description  of  SLAM  syirbols  and  con¬ 
cepts.  The  model  uses  a  combination  discrete  event-network 
orientation  to  simulate  the  aircraft's  flight  and  consists 
of  two  parallel  networks  representing  the  flights  of  the 
tanker  and  airlift  aircraft.  The  output  variable  and  mea¬ 
sure  of  effectiveness  is  fuel  consumed  by  the  tanker  and 
airlifter.  Only  the  airlift  network  is  described  here  and 
it  is  broken  down  into  the  segments  of  Figure  K-1.  However 
differences  in  the  tanker  network  will  be  noted.  An 
in-depth  nodal  description  follows  the  general  description. 

General  Description 

The  SLAM  model  receives  some  of  the  input  param¬ 
eters  shown  in  Figures  K-3  and  K-4  by  calling  the  flight 
plan  analytic  model  (FLTPLN) .  The  global  variables  (XX(i)) 
contain  the  initial  conditions  and  other  pertinent  data. 
Attributes  1-5  are  used  as  shown  in  Table  K-1. 

The  first  segment  of  the  model  represents  the 
start,  taxi,  and  takeoff  phase  (node  ASl  to  AS 7  for  the 
airlifter  and  AS2  to  AS6  for  the  tanker) .  Delays  can  be 
incurred  by  each  aircraft.  A  logic  sequence  (ASS  to  AS7 
for  the  airlifter)  prevents  the  second  aircraft  from  depart 
ing  until  the  first  scheduled  aircraft  has  departed. 

The  second  segment  of  the  model  represents  the 
initial  climb  from  the  base  to  the  cruise  altitude  and  the 
cruise  to  the  rendezvous  entry  point  (AS7  to  ASX  for  the 
airlifter  and  AS6  to  ASY  for  the  tanker) .  Included  in  this 


XX  (1) 

XX  (2) 

XX(3) 

XX(4) 

XX  (5) 

XX(6) 

XX  (7) 

XX(8) 

XX  (9) 

XX(IO) 

XX(ll) 

XX(12) 

XX(13) 


XX(15) 
XX (16) 
XX(19) 
XX(20) 
XX(21) 

XX(22) 
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Airlift  departure  to  rendezvous  entry  point  dis¬ 
tance. 

Airlifter  rendezvous  exit  point  to  destination 
point. 

Airlifter  rendezvous  entry  point  to  abort  base 
distance . 

Tanker  departure  to  rendezvous  entry  point  dis¬ 
tance. 

Tanker  rendezvous  exit  point  to  destination. 

Transfer  fuel. 

Initial  fuel  load — tanker. 

Initial  fuel  load — airlifter. 

Gross  weight  of  the  tanker  at  time  (t)  . 

Gross  weight  of  the  airlifter  at  time  (t) . 

Tanker  takeoff  time. 

Airlifter  takeoff  time. 

Aircraft  code: 

C-5/KC-135  =  1 
C-141B/KC-135  =  2 
C-5/KC-10  =  3 
C-141B/KC-10  =  4 

Airlifter  takeoff  delay. 

Tanker  takeoff  delay. 

Initial  climb  distance  airlifter. 

Initial  climb  distance  tanker. 

Arrival  time  of  airlifter  at  rendezvous  entry 
point  (Ql) . 

Arrival  time  of  tanker  at  rendezvous  entry  point 
(Q2)  . 


Fig.  K-3.  Global  SLAM  Variables 
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XX(23)  =  Airlifter  fuel  consumption  from  reidezvous  entry 
point  to  abort  base. 

XX (24)  =  Tanker  fuel  consumption  from  rendezvous  entry 
point  to  abort  base, 

XX (27)  =  Distance  traveled  in  the  refueling  track. 

XX (31)  =  Airlifter  departure  base  elevation. 

XX (32)  =  Airlifter  initial  cruise  altitude. 

XX (28)  =  Time  required  in  the  refueling  track 

XX (33)  =  Airlifter  cruise  time  from  leveloff  to  rendezvous 
point. 

XX(35)  =  Airlifter  initial  climb  fuel. 

XX (40)  =  Tanker  initial  climb  fuel. 

XX (41)  =  Tanker  departure  base  elevation. 

XX (42)  =  Tanker  initial  cruise  altitude. 

XX (45)  =  Airlifter  destination  fuel  reserves. 

XX (46)  5=  Tanker  destination  fuel  reserves. 

XX (60)  =  Tanker  post-rendezvous  climb  fuel. 

XX(61)  =  Airlifter  post-rendezvous  climb  fuel. 

XX(62)  =  Tanker  rendezvous  exit  point  to  destination  time. 

XX (63)  =  Airlifter  rendezvous  exit  point  to  destination 
time. 

XX(71)  =  Airlifter  cruise  altitude  to  destination. 

XX (73)  =  Airlifter  average  cruise  altitude  to  destination. 
XX(75)  =  Airlifter  ramp  fuel. 

XX (76)  =  Tanker  ramp  fuel. 

XX (81)  =  True  course  for  leg  1. 

XX (82)  =  True  course  for  leg  2. 

Fig.  K-3 — Continued 
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XX (83)  =  True 
XX (84)  =  True 
XX (8 5)  =  True 
XX(94)  =  Zero 
XX (9 5)  =  Zero 


course  for  leg 
course  for  leg 
course  for  leg 
fuel  weight  of 
fuel  weight  of 

Fig.  K-3 


3. 

4. 

5. 

the  airlifter 
the  tanker. 

-Continued 
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TANKER  MODEL 


TABLE  K-1 


SLAM  ATTRIBUTES  SUMMARY 


ATRIB(i) 

Description 

(1) 

Aircraft  takeoff  time 

(2) 

Aircraft  conJjination  code 

(3) 

Cruise  time 

(4) 

Cumulative  fuel  consumed 
for  a  particular  aircraft 

(5) 

Aircraft  identifier 

1  =  Airlifter 

2  =  Tanker 

NOTE;  Attributes  are  numerical  characteristics 
assigned  to  entities  (aircraft)  traveling  through  the  net¬ 
work  simulation  model.  These  attributes  travel  with  a  spe¬ 
cific  aircraft  and  are  unique  to  that  aircraft.  This 
table  shows  the  attributes  used  in  the  SLAM  model. 


segment  are  fuel  consumption  calculations  for  climb, 
cruise,  and  loiter  times.  The  loiter  time  is  the  time  dif¬ 
ference  between  the  arrival  of  the  second  and  first  air¬ 
craft  at  the  rendezvous  point  (Ql  and  Q2) .  The  aircraft 
fuel  supply  must  be  sufficient  to  return  from  the  rendezvous 
point  to  a  suitable  abort  base,  if  the  refueling  is  aborted. 

The  third  segment  (ASX  to  ASE  for  the  airlifter 
and  ASY  to  ASF  for  the  tanker)  represents  the  refueling 
maneuver  and  calculates  the  time  and  fuel  consumed  during 
the  fuel  transfer. 

The  fourth  segment  (ase  to  T3  for  the  airlifter  and 
ASF  to  T4  for  the  tanker)  represents  the  route  segment  from 
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the  rendezvous  exit  point  to  landing.  Inc]uded  in  these 
segments  are  time  and  fuel  calculations  for  post-rendezvous 
climbout,  cruise  to  destination,  and  approach  and  landing. 

At  the  terminal  nodes,  T3  and  T4,  the  total  fuel  consumed 
respectively  by  the  airlifter  and  tanker  aircraft  are 
printed. 

Nodal  Description 

This  section  describes  the  mission  profiles  through 
the  SLAM  model  node  by  node.  The  discrete  flight  planning 
model  provides  takeoff  times  for  the  airlifter  and  tanker 
and  calls  the  aircraft  network  model.  (Appendix  E  is  the 
SLAM  computer  code.)  A  logical  IF  statement  in  subroutine 
FLTPLN  will  call  the  ENTERl  or  ENTER2  node  depending  on 
which  aircraft  takes  off  first.  For  simultaneous  takeoffs, 
only  one  ENTER  node  will  be  called  so  multiple  aircraft 
are  not  generated.  The  ENTER  node  simulates  engine  start 
for  the  first  aircraft.  If  the  airlifter  takes  off  first, 
the  ENTERl  node  is  entered  and  the  next  node  (ASl)  assigns 
the  value  of  one  to  global  variable  II,  indicating  the  air¬ 
lifter  should  take  off  first.  II  is  assigned  a  value  of  two 
in  the  tanker  network  if  the  tanker  is  scheduled  to  take 
off  first. 

The  next  event  consists  of  two  regular  activities, 
one  continues  the  airlifter  to  takeoff  and  one  initiates 
the  tanker  engine  start  sequence.  The  activity  from  ASl 
to  AS3  involves  no  time  and  at  ASS  assigns  ATRIB(5)  equal 
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to  the  aircraft  code  to  facilitate  reading  the  SLAM  Output 
Trace.  The  second  activity  leading  from  ASl  occurs  with 
a  completion  time  from  XX (11)  which  is  the  takeoff  time 
difference  between  aircraft.  In  this  way  each  aircraft 
is  assured  of  completing  the  takeoff  sequence  at  the 
correct  time. 

The  next  section  models  takeoff  delays  encountered. 
The  fifteen  minutes  activity  from  AS3  to  either  GOll  or 
G033  simulates  the  time  for  start,  taxi,  and  takeoff.  This 
is  a  standard  value  used  in  flight  planning.  The  decision 
to  go  to  GOll  or  G033  is  based  upon  the  aircraft  code  in 
ATRIB(2).  If  the  aircraft  is  a  C-5  the  aircraft  will  delay 
in  accordance  with  the  network  following  GOll;  if  the  air¬ 
craft  is  a  C-141B  the  aircraft  will  delay  through  G033. 

If  a  delay  in  the  scheduled  takeoff  time  for  the  first  air¬ 
craft  occurs,  the  second  aircraft  will  remain  grounded 
until  the  first  aircraft  is  repaired  or  the  mission  is  can¬ 
celled.  Stops  to  prevent  the  second  scheduled  aircraft 
taking  off  prior  to  the  first  are  discussed  later.  The 
delay  occurs  probabi'' istically  and  always  for  fifteen 
minutes.  Those  delays  which  do  not  burn  fuel  are  not 
modeled.  Once  the  airlift  delay  is  determined,  it  is 
assigned  to  global  variable  XX (15)  which  will  be  used  in 
future  fuel  computations. 

The  next  segment  computes  the  fuel  used  in  the 
start,  taxi,  and  takeoff  event  including  the  delays. 
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ATRIB(4)  is  used  as  a  fuel  accumulator  for  each  aircraft. 


It  is  initially  set  at  zero  and  all  fuel  consumed  is  added 
to  ATRIB(4)  as  it  is  calculated.  ASS  calculates  the  fuel 
consumed  up  to  that  point  by  calling  user  function  one 
which  assigns  the  start,  taxi,  takeoff  fuel  for  the  air- 
lifter  and  calculates  the  amount  of  fuel  used  in  any  delays. 
See  Appendix  F  for  user  function  descriptions.  ASS  also 
assigns  the  value  of  "one"  to  XX (17)  which  will  be  used  in 
determining  takeoff  sequences. 

To  determine  when  to  actually  launch  the  second  air¬ 
craft,  the  activity  network  emanating  from  ASS  checks  three 
questions  in  this  order: 

1.  Has  the  tanker  taken  off  first? 

2.  Should  the  tanker  take  off  first  and  is  he 

delayed? 

3.  Should  the  airlifter  take  off  first? 

If  the  first  answer  is  yes,  the  airlifter  takes 
off  and  travels  to  node  AS7.  If  the  answer  is  no  and  the 
second  answer  is  yes,  the  tanker  was  delayed  and  the  air¬ 
lifter  must  also  delay  to  provide  the  correct  takeoff 
separation  between  the  two  aircraft.  The  airlifter  will 
stay  grounded  until  the  tanker  is  airborne  plus  the  original 
takeoff  time  difference.  This  is  accomplished  by  adding 
the  original  takeoff  separation  to  the  actual  time  the 
first  scheduled  aircraft  departs.  If  the  second  answer  is 
no,  the  third  answer  is  yes  and  the  airlifter  departs 
first.  The  takeoff  maneuver  is  modeled  as  three  minutes 
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(brake  release  to  landing  gear  retraction) .  The  SLAM  model 
uses  a  conditional  take  first  logic  on  the  three  activi~ 
ties,  thus  allowing  only  the  first  "yes"  activity  to  be 
released.  The  logic  sequence  between  ASZ  and  AS6  deter¬ 
mines  the  takeoff  time  and  sequence  for  the  tanker.  After 
the  takeoff  is  accomplished  AS7  computes,  with  USERF3 , 
the  takeoff  fuel  and  adds  it  to  ATRIB(4).  AS6  performs 
this  function  for  the  tanker. 

The  next  phase  of  the  flight  is  the  initial  climb 
to  altitude.  The  airlift  climb  time  is  computed  in  USERF5. 
AS9  then  calculates  the  climb  fuel  and  climb  distance 
using  USERF7  and  USERF9  respectively.  USERF8  and  10  cal¬ 
culates  the  values  for  the  tanker.  The  climb  fuel  is  added 
to  ATRIB(4)  in  AS9  for  the  airlifter  and  in  ASS  for  the 
tanker. 

The  next  phase  of  flight  is  the  cruise  from  level- 
off  to  the  rendezvous  entry  point.  The  cruise  time  for  the 
airlifter  is  confuted  in  USERFll  and  for  the  tanker, 

USERF12  is  used.  ASB  computes  cruise  fuel  to  the  entry 
point  (Ql) ,  adds  it  to  ATRIB{4)  and  sets  XX{21)  equal  to 
the  entry  point  arrival  time.  XX (21)  will  be  used  to  deter' 
mine  aircraft  loiter  time.  If  the  tanker  arrives  first  it 
will’ hold  at  Q2 . 

The  next  event  is  the  rendezvous  maneuver  and 
refueling.  The  maneuver  begins  with  the  arrival  of  the 
second  aircraft  at  the  rendezvous  point  Ql  or  Q2.  When 
the  first  aircraft  arrives  at  the  queue  node,  it  is 
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prevented  from  proceeding  further  by  match  node,  MATl. 

The  arrival  of  the  second  aircraft  releases  MATl  and  the 
refueling  maneuver  begins.  Once  through  the  match  node, 

ASC  and  USERF15  computes  the  loiter  fuel  for  the  airlifter. 
ASD  and  USERF16  perform  this  function  for  the  tanker. 

After  each  aircraft  has  arrived  at  the  rendezvous 
point,  but  before  actual  refueling  occurs,  two  conditions 
must  be  met: 

1.  Aircraft  must  be  capable  of  safely  returning 
to  a  suitable  abort  base,  if  hookup  is  unsuccessful. 

2.  Aircraft  must  be  able  to  fly  to  the  destina¬ 
tion  base  with  the  fuel  available  after  the  transfer. 

If  these  conditions  are  not  met,  then  unexpected  fuel  con¬ 
sumption  due  to  excessive  loiter  fuel,  delays,  or  unfore¬ 
seen  wind  changes,  have  made  the  refueling  event  infeasible 
due  to  a  lack  of  available  fuel.  Each  aircraft  will  abort 
to  their  respective  bases.  The  abort  computation  is  calcu¬ 
lated  in  USERF17  for  the  airlifter  and  USERF18  for  the 
tanker.  The  SLAM  model  compares  the  difference  between  the 
fuel  available  and  the  fuel  required  to  abort.  A  differ¬ 
ence  of  less  than  zero,  represents  an  infeasible  situation 
and  the  simulation  stops  by  aborting  the  applicable  air¬ 
craft  to  terminal  node  1  or  2. 

If  the  aircraft  meet  the  feasible  criteria,  they 
proceed  with  the  refueling  maneuver.  Refueling  occurs  on 
the  refueling  track,  prearranged  in  length  by  subroutine 
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FLTPLN.  The  time  in  the  maneuver  is  identical  for  each  air¬ 
craft  and  calculated  by  ASX  and  USERF19  for  the  airlifter, 
and  ASY  and  USERF19  for  the  tanker.  Fuel  consumed  in  the 
rendezvous  maneuver  is  calculated  in  ASE  and  ASF  for  the 
airlifter  and  tanker  using  USERF21  and  20  respectively. 

Once  the  rendezvous  is  completed  each  aircraft 
departs  the  exit  point  to  their  respective  destinations. 

The  time  to  destination  is  computed  by  USERF23  and  22  for 
the  airlifter  and  tanker  respectively  and  assigned  to 
ATRIB(3).  The  fuel  to  reach  each  destination  is  computed 
by  USERF25  and  24  for  the  airlifter  and  tanker  respectively 
and  added  to  the  total  fuel  consumed  in  ASG  and  ASH. 

Approach  and  landing  fuels  and  fuel  reserves  are  included 
in  these  fuel  computations.  Fuel  reserves  used  are  from 
the  flight  plan  parameters. 

When  the  airlifter  and  tanker  reach  their  destina¬ 
tions,  they  pass  through  their  terminal  node  T-3  (airlifter) 
and  T-4  (tanker) .  When  both  terminal  nodes  are  released, 
the  simulation  is  stopped  and  the  results  are  displayed  in 
a  SLAM  Nodal  Trace  of  the  simulation  from  takeoff  to  land¬ 
ing.  Total  fuel  consumed  by  each  aircraft  is  obtained  from 
the  column  labeled  attribute  four  at  T3  and  T4 .  Appendix  G 
lists  a  typical  SLAM  nodal  trace. 
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Appendix  L 

FLTPLN  and  SLAM  Summary  of  Results 


i 
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This  appendix  contains  four  tables  of  comparisons 


between  SLAM  and  FLTPLN  fuel  consumptions.  Table  L-1  shows 
these  results  for  the  C-5,  Table  L-2  for  the  C-141B, 

Table  L-3  for  the  KC-10,  and  Table  L-4  for  the  KC-135. 

The  data  in  each  table  represents  the  combined  fuel  con¬ 
sumption  data  for  all  scenarios  (McGuire  AFB  to  Tehran, 

Iran  and  Travis  AFB  to  Yokota  AB,  Japan)  for  a  particular 
aircraft.  The  SLAM  data  was  obtained  using  a  normal  wind 
factor  variant  (y  =  55,  a  =  10  KNOTS) . 

In  Tables  L-1  to  L-4  the  first  column  (Tanker  Base) 
lists  the  tanker  departure  base.  By  specifying  only  the 
tanker  base  the  scenario  is  also  fully  defined  since  only 
the  east  coast  tanker  bases  are  used  in  the  Tehran  scenario 
and  only  the  west  coast  tanker  bases  are  used  in  the  Yokota 
scenario.  The  second  column  (Complementary  Aircraft) 
defines  the  aircraft  combination.  The  third  column  (Air- 
lifter  Cargo)  represents  the  airlifter  cargo  load;  those 
missions  which  were  infeasible  at  a  particular  cargo  weight 
are  marked  "Infeas."  The  fourth  column  (Analytical  Fuel 
Consumption)  represents  the  total  fuel  consumption  for  the 
applicable  aircraft  obtained  from  the  analytical  model 
(FLTPLN) .  The  fifth  column  (SLAM  Fuel  Consumption)  repre¬ 
sents  the  SLAM  fuel  consumption  sample  mean.  The  sixth 
column  (%  Dev)  represents  the  percent  deviation  of  the 
SLAM  mean  compared  to  the  analytical  model  value.  For 
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exanple,  the  KC-10  fuel  consumption,  for  the  Travis-Yokota 
scenario,  refueling  the  C-141  with  a  cargo  load  of  55,000 
lbs,  was  138,000  lbs.  The  SLAM  mean  fuel  consumption  was 
131,200  lbs.  The  percent  deviation  is: 


131.2 

138.2 


-5.5% 


A  negative  value  means  that  the  SLAM  model  used  less  fuel 
than  the  analytical  model  and  a  positive  value  means  the 
SLAM  model  consumed  more  fuel.  The  seventh  column  shows  a 
95  percent  confidence  interval  as  a  total  fuel  consumption 
range  such  that  the  probability  of  the  true  SLAM  mean  being 
in  that  interval  is  95  percent.  Appendix  J  contains  a 
sample  confidence  interval  calculation.  To  obtain  a  con¬ 
fidence  interval,  normality  must  be  assumed.  Kolmogorov- 
Smirnov  Tests  of  the  SLAM  data  were  used  to  verify  normality 
and  Appendix  J  contains  a  normality  verification  sample 
calculation  using  the  Vogelbach  Computing  Center’s  SPSS 
Kolmogorov-Smirnov  Nonparametric  One  Sample  Test.  A  "yes” 
in  column  eight  (SLAM  Feasible)  indicates  that  the  SLAM 
total  mission  fuel  consumption  for  the  specific  aircraft 
using  the  upper  limit  of  the  confidence  interval  as  the  fuel 
required  for  a  successful  mission,  was  feasible  for  the  par¬ 
ticular  scenario.  The  fuel  available  for  the  airlifter  is 
the  ramp  fuel  plus  the  transferred  fuel.  The  fuel  avail¬ 
able  for  the  tanker  aircraft  is  the  ramp  fuel  minus  the 
transferred  fuel.  In  Figures  3-6  the  percent  deviation  for 
each  aircraft  is  plotted. 
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SLAM  VERSUS  ANALYTICAL  FUEL  CONSUMPTION  RESULTS 
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